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I.  Introduction 


This  final  report  covers  research  carried  out  during 
the  period  09/01/79  to  09/09/82  under  U.S.  Army  Contract 
DAAG  29-79-0203.  The  principle  goal  of  the  work  was  the 
development  of  high  power  Cerenkov  sources  in  the  lower 
mm  wavelength  range.  It  was  demonstrated  that  a  Cerenkov 
maser,  which  consists  of  a  mildly  relativistic  electron 
beam  and  a  tubular  dielectric  resonator,  is  capable  of 
producing  power  at  the  multi-hundred  KW  level  with 
efficiencies  in  the  ten  percent  range.  A  detailed  dis¬ 
cussion  of  this  device  is  found  in  Section  III  c) , 
Characteri zation  of  a  Cerenkov  Maser,  by  S.  Von  Laven. 

Also  contained  in  Section  III  are  a  general  review  of  the 
process  of  stimulated  Cerenkov  radiation,  and  a  paper 
which  examines  the  concept  of  submm  FIR-tunable  Cerenkov 
devices.  Published  abstracts  are  reproduced  in  Section  II 
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Abstract  Submitted 
for  the  Washington  Meeting  of  the 
American  Physical  Society 


April  23-26,  1979 


Suggested  title  of 
session  - 

Electron  Beans  and 
free'  Electron  Lasers 


Stimulated  Raman  Scattering  from  Electron  Beams 
in  Dielectric  Resonators .  K.  BUSBY,  K.  FELCH,  R.W. 
LAWMAN,  J.E.  WALSH,  Dartmouth  College. -'--The  addition 
of  a  dielectric  resonator  to  an  electron  beam  Raman 
scattering  experiment  can  dramatically  shorten  the  wave¬ 
length  of  the  output  radiation.  Well  above  the  guide 

cutoff  frequency  the  output  wavelength  A  is. 

s 

eon)/6o 

where  \p  is  the  punp  ripple  length  and  80  is  average 
electron  beam  velocity.  A  combination  Cerenkov- Raman 
scattering  experiment  has  been  designed  to  test  this 
reduction.  We  use  a  tubular  dielectric  resonator  and  a 
rippled  magnetic  field  as  a  pump  field.  The  effective 
intensity  of  the  pixitp  field  is  approximately  20  IW/cm2. 
This  together  with  electron  beam  voltages  and  currents 
of  100-200  KV  and  10  A  respectively  will  lead  to  stimu¬ 
lated  Cerenkov- Raman  growth  lengths  of  .2-. 5  cm-l.  The 
device  performance  will  be  discussed. 

*Vfork  supported  in  part  by  U.S.  Army  Grant  #DAAG39-78- 
C-0032  and  in  part  by  U.S.AP0SR  Grant  #  77-3410 
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Abstract  submitted  for  the 
1979  IEEE  INTERNATIONAL  CONFERENCE  ON  PLASMA  SCIENCE 

June  4-6,  1979 


A  Cerenkov-Raman  Radiation  Source.  K.  BUSBY, 

K.  FELCH,  R.W.  LAYMAN,  J.E.  WALSH,  Dartmouth  College.* 
■—Stimulated  Raman  Scattering  from  an  electron  beam  in 
a  dielectric  medium  shows  a  number  of  interesting  and 
potentially  useful  features In  a  dielectric  the 
particle  velocity  can  exceed  the  speed  of  light  and 
hence  there  are  two  possible  dispersion  curves  along 
which  energy  and  momentum  can  be  conserved.  Further¬ 
more,  in  the  neighborhood  of  the  Cerenkov  velocity 
threshold  very  large  Doppler  shifts  can  be  obtained 
and  thus  this  interaction  can  be  the  basis  of  a  short 
wavelength  radiation  source. 

In  order  to  test  the  potential  of  such  a 
device  an  experiment  which  uses  a  mildly  relativistic 
electron  beam,  a  dielectric  resonator,  and  a  rippled 
magnetic  field  for  a  pump  wave  has  been  constructed. 
The  transverse  component  of  the  magnetic  field  is 
approximately  300  G  which  implies  an  effective  pump 
wave  intensity  of  20  MW/cm^. 

With  beam  voltages  between  100  and  200  KV  and 
beam  currents  in  the  10A  range,  operation  in  the 
lower  and  middle  part  of  the  mm  wave  length  region  of 
the  spectrum  is  anticipated.  The  performance  of  the 
device  will  be  discussed. 

*Work  supported  in  part  by  U.S.  Army  Grant  #DAAG39-78- 
-C0032  and  in  part  by  U.S.AFOSR  Grant  #  77-3410 

1.  J.E.  Walsh,  G.  Crew,  Bull.  Am.  Phys.  Soc.  23,  748, 
(1978) 
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Generation  of  Millimeter  Microwaves 


Usinqthe  Cerenkov  Interaction.  K.  Felch. 
Busby,  jJeI  Walsh ,  and  R.w.  Layman,  Dartmouth 
College.  — Earlier  experiments  involving  the 
Cerenkov  interaction  between  a  relativistic 
electron  beam  and  a  dielectric  loaded  waveguide 
yielded  6  mm  microwaves.  These  experiments 
indicated  that  microwaves  of  significantly 
shorter  wavelengths  might  be  possible  using  the 
same  procedures.  The  present  experiments 
extend  the  previous  attempts  into  the  3  mm 
range.  Methods  by  which  shorter  wavelengths 
are  achieved  using  the  Cerenkov  interaction  and 
the  short  wavelength  limit  of  such  techniques 
will  be  discussed.  Results  of  the  latest 
experimental  investigations  will  be  presented. 

*Wor ^supported  in  part  by  U.S.  Army  Grant 
#DAAG39-78-C-0032  and  in  part  by  U.S.  AFOSR 
Grant  #77-3410. 
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Multiple  Frequency  Output  of  A  Raman 
Maser  With  A  Magnetic  Higgler  PunpV  K.  BUSBY, 

K.  FELCH,  R.W.  LAYMAN ,  D.  KAPILOW,  J.E.  WALSH, 

JL  JL  , 

Dartmouth  C. — Microwave  output  at  two 
distinct  frequencies  from  a  Raman  maser  has 
been  observed.  The  microwave  production  is 
due  to  the  interaction  of  a  40-100  KeV,  10  A 
electron  beam  with  a  0.9  cm  period,  zero 
frequency,  magnetic  wiggler  field.  The  undulations 
are  a  15%  perturbation  to  the  guiding  magnetic 
field.  The  interaction  takes  place  within  a 
waveguide. 

The  electron  velocities  are  too  low  to  inter¬ 
act  with  the  fundamental  wavelength,  A_  =  0.9  cm, 

•J  j-J 

undulations,  but  the  electrons  do  interact  with 
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higher  harmonic  wavelengths.  Only  the  odd  harmonics 
need  be  considered  since  the  magnetic  wiggler 
appears  approximately  as  a  square  wave. 

As  shown  in  the  figure,  synchronism  occurs 
when  the  negative  energy  cyclotron  beam  mode  dis¬ 
persion  (u>b  =  -Q  +  k.VQ>  couples  to  the  magnetic 
wiggler  dispersion  (Wpuap  =  0,  kpurnp  =  2-nApump, 
n  **  1,3,5,...)  such  that  the  waveguide  dispersion 
is  satisfied.  The  two  horizontal  lines  are  the 
third  and  fifth  harmonics  of  the  magnetic  wiggler 
dispersion. 


u> 


synchronism 


waveguide  dispersion 


magnetic  wiggler 


electron  beam  mode 


electron  velocity 
(l-Vo2/c2)_Js 


With  73  KeV  electrons,  Bg  =  6200  gauss,  and  a  TEgi 
cutoff  frequency  of  29.1  GHz,  the  frequencies  at 
synchronism  are  30  and  48  GHz,  the  observed  fre¬ 
quencies  are  29  and  42  GHz.  Convective  instabilities 
at  (+|w| ,  +  |k|)  are  possible,  but  are  at  higher 
frequency. 

Gain  calculations,  in  the  collective  regime, 
for  the  Raman  interaction  between  a  relativis  ic 
electron  bean  and  a  static  magnetic  higgler  .eld 
will  be  presented.  Waveguide  effects  v'1.1  e 

included. 

♦Submitted  by  R.W.  Layman 

♦♦Supported  in  part  by  U .  S .  Amy  Grnnt-DAAf'3°-78-C- 
0032  and  in  part  by  U.S.  AF0SR  Grant  !■  77-3410. 


Subject  category  number: 


•  (  X  )  Prefer  oral  session 

•  (  )  Prefer  poster  session 

•  (  )  No  preference 

•  (  )  Special  requests  for  placement 

of  this  abstract 

If  more  than  one  session 

is  offered,  we  want  to  be 
in  the  first  session. 


•  Submitted.^:  /) 

v  (Signature)  / 

R.W.  layman _ 

(same  name  typewr  :;en) 

Physics  Department 
(full  address) 

Dartmouth  College 


Hanover .  NH  03755 


•  I  am  a  member  of  the  C mm  ‘me  or; 
Plasma  Sc'once  and  App  ...aliens 


(  )  yes 


(  X  '  no 


Abstract  Submitted 

For  the  Twenty-second  Annual  Meetin » 
Division  of  Plasma  Physics 
November  10  to  14,  1980 

Subject  Category  Number  4 . 8 _ ’ _ 


Stimulated  Cerenkov  Radiation  at  Lower  mm 
Wavelengths ,  J.  WALSH,  D.  KAFILCW,  R.  LAYMAM,  D.  SPEER, 
J.  BRANSCUM,  D.  WISE  and  C.  HANNA,  Dartmouth  College — 
Stimulated  Cerenkov  radiation  is  generated  by  a  mildly 
relativistic  electron  beam  dielectric  resonator 
combination.  The  200KV,  5-10A  1  usee  electron  beam 
pulse  is  produced  by  a  transformer  and  a  pulse  forming 
network.  The  radiation  is  contained  in  two  modes,  the 
first  in  the  4-8  mm  region  and  the  second  between  1  and 
2  mm.  The  operating  wavelength  is  determined  by  the 
beam  voltage  and  the  resonator  dimensions.  The  wave¬ 
lengths  are  measured  with  internal  reflection  filters, 
cutoff  filters,  and  a  mechanically  tuned  Fabry-Perot 
resonator.  The  effects  on  device  performance  of  a 
variety  of  quasi-optical  feed  back  and  mode  control 
techniques  will  be  discussed. 
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Stimulated  Cerenkov  Radiation  Experiments. 

J.E.  WALSH,  J.  BRAXSCUM,  J.  GOLUB,  D.  SPEER,  D.  WISE, 
and  R.  LAYMAX,  Dartmouth  College.* — An  electron  beam 
generator  with  a  voltage  range  of  up  to  400  kv  and  an 
output  current  of  0. 1A  to  20A  has  been  used  to  conduct 
stimulated  Cerenkov  radiation  experiments.  A  tubular 
quartz  dielectric  resonator  is  used  as  a  slow  wave 
structure  for  TM()n  waves  which  interact  synchronously 
with  the  electron  beam.  Using  a  series  of  resonator 
designs  coherent  oscillations  have  been  obtained  in  the 
wavelength  range  between  1.5  cm  and  1.5  mm.  Power 
output  at  the  longer  wavelengths  is  in  the  range  of 
tens  of  kws  while  that  at  shorter  wavelengths  it  is 
approximately  one  kw.  At  the  present  time  the  beam  is 
collected  on  a  mirror  set  at  a  45°  angle  to  the  end  of 
the  resonator.  The  mirror  reflects  the  output  through 
a  collimating  lens  to  a  Fabry-Perot  interferometer. 

A  description  of  the  apparatus  and  experiments 
with  various  resonator  parameters  such  as  dielectric 
constant,  radial  dimensions  and  length  will  be 
described.  A  discussion  of  theory  and  the  correlation 
of  theory  and  experiment  will  be  presented. 

*Work  supported  in  part  by  U.S.  Army  Grant  //DAAG  29-79- 
C0203 
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which  la  linearly  dependent  on  the  well  depth  D,  any 
Instability  in  D  affects  the  stability  of  the  v!  reso¬ 
nance*  A  lower  potential  with  greater  stability  should 
allow  us  to  obtain  i  single  or;  tor.  iinvwidth  of  G.  D05  Hz 
limited  by  detection  circuit  darning.  Ccruari.vrn  of 
proton  and  electron  cyclotron  resonances  yield  a  pre¬ 
liminary  miss  ratio  r»/ae  ■  1836. 15 >00(23)  wntre  the  un¬ 
certainty  has  been  increased  by  0.1  spa  to  cover  a  possi¬ 
ble  systematic  error.1  This  value  is  20  tines  more 
accurate  than  our  previously  reported  value. - 
^R.S,  Van  Eyck,  Jr.,  P.5.  Sohwinber,*,  ar.d  5.K.  Salley  in 
Atomic  Masses  In  run  1  i-ental  C^stmt  »  6  (Plenum,  i 0 9. 0 ) . 
^SL.S.  Van  Dyck,  Jr.,  ana  ?.B.  Schwinberg,  to  be  publish¬ 
ed  in  Bull.  Am.  Phys.  Soc.  26  (1981). 


■SESSION  E:  SPECTROSCOPY 
Saturday  morning,  13  April  1931 
Salisbury  Labs,  Room  121  at  9:00  A.M. 

T.  H.  foil,  presiding 

E  l  An  Exact  Solution  to  the  Problem  of  Sal f -Diffusion 

Via  the  Vacancy  Mechanic  in  the  Sudden  Jump  Appro  ■: mat  ion, 
B.  CRALLA  and  L.  JVNIHLR,  Tufts  L*. — An  exact  solution  has 
been  obtained  for  the  problem  of  self-diffusion  bv  the 
vacancy  mechanism  in  the  sudden  lump  dpproximat ion4 .  fhe 
probability  distribution  function  P(n,t)  for  locating  a 
marked  atom  at  a  lattice  posit ior  r  at  a  time  t  can  be  ex¬ 
pressed  as  an  integral.  The  theory  is  relevant  over  time 
scales  on  the  order  of  the  mean  jump  tine  or  greater. 
Therefore,  not  only  can  we  obtain  the  diffusion  constant 
but  we  can  also  calculate  the  effect  of  self-diffusion  of 
a  Mossbauer  atom  on  Che  Mdssbauer  spectrum.  In  the  case 
of  the  linear  chain  we  have  a  solution  in  terms  of  known 
functions.  In  particular  the  mean  square  displacement 
<n2>  of  a  marked  atom  is  proportional  to  t1'-  as  c  app¬ 
roaches  infinity,  which  is  a  well-known  result'.  The 
Bardeen-Herring  correlation  factor  f  for  the  diffusion 
constant  for  a  square  lattice  has  been  found  equal  to 
(it— l)"” ^ ,  in  agreement  with  the  value  obtained  bv  Scnoen 
and  Lowen3  using  direct  combinator ial  techniques.  Work 
is  in  progress  in  calculating  the  correlation  factors 
for  three-dimensional  lattices  and  on  calculating  the  co¬ 
rresponding  Mossbauer  spectra. 

1  M.C .  Dibar  Ire  anJ  ?.A.  r'Unn,  Phys.  Rev.  15B,  1261 
(1977). 

*  D.L.  Huber,  Phys.  Rev.  B!5,  533  (1977). 

*  A.H.  Schoen  and  R.W.  Low«_n.  Bull.  Am.  Phvs.  Soc.  5_, 

280  (I960). 


>  v(12)  >  ...  refer  tc  the  r*  ■  0  l,  n  ■  l  •*  2,  ... 

sitior.s  spaced  b/  150  Hz,  ani  H  ( C 1  j  ,s  f.r  a  r.  - 

0  -  I  Rabi  frequency  of  53  Hz.  \ow,  *  -r  m  -  s- 
m.^rfect  adiaba'-ic  fast  passages  roving  t r.e  *  *?*r  *>r.  frnr 
n  =»  0  -•  1,  1  *»  2,  2  -  1  ...  the  excit.n; 

becomes  tna  synchrocyclotron  rr-cr.ir.i  * .  ->  :  ■ 

of  the  nucicn  folir»z.vj  vtt‘.  Thu  corrcs  r.i,:* ;  *■  ,  :  in¬ 

crease  is  easily  -tsvtvea  vie  a  a-.ift  *  '  rz.  far 
tr.e  electron  initially  in  n  ■*  1,  xr.  tm.e  or c.er  .-.at.:,  in-* 
first  adiabatic  fast  passage  moves  it  ;  rvi  r  -  v-.  -r  i  _r 

it  lerains  as  the  n  *  1  -*■  i  tzar  sit  ur.  ca..r;t  70  r. .  In 

an  alternate  faster  senene  the  n-state  vnerqy  ii  j-.  Ie:- 

tively  multiplied  in  ''•10  )  usee  by  tne  ir.hc- 14*  4rv-»ld 

of  a  microwave  at  0 1 )  -  30  kHz}.  The  ex:  its*. :  :n  fc  ir.  1 
is  carefully  aaju  it  a  from  v(01)  to  f.f21)  -  3'  *nz!. 

Now,  on  the  average,  the  n  *  1  state  energy  ;s  m cl-.i:  lied 
xlOO  while  the  n  *  G  state  energy  tfjtni  t.ne  ii>.  de¬ 
tection  is  again  via  mass  increase.  We  thank  the  National 
Science  Foundation  for  support. 


(SESSION  F:  NUCLEAR/PARTICLE  PHYSICS 
Saturday  morning,  1-5  April  1981 
Olin  Hall.  Room  107  at  9:00  A.M. 

V.  Bluemel,  presiding 

F  I  Achromatic  Proton  Microbe  tm,  F.  W.  YA  Si  IS, 

M I c re -•  copy  A.;  c  a ,  Inc,  and  G01£S£!E, 

Polytechnic  Ir<gt  1 r. ute . --.At.  achromatic  quadrupole  lens 
doublet  is  Scir.g  evaluated  as  a  means  or  forming  *  iv.elv 
focused  ion  beans.  Experimental  measurements  indicate 
that  a  small  image  broadens  less  than  5  microns  w  ea  v:i* 
energy  of  transmitted  protons  is  varied  fron  5'0  zn 
800  keV.  A  chromatic  aberration  coefficient  of  less 
than  0.4  cm  is  derived.  Typical  measurements  0:  beam 
dimension,  made  using  the  Rutherford  backseat  ter  :ny 
technique,  wiLL  be  presented.  Procedures  will  be 
discussed  for  minimizing  parasitic  first-  and  second- 
order  aberrations  by  electrical  rotation  of  the  prin¬ 
cipal  sections  of  a  lens*  and  by  means  of  an  off-center 
aperture  placed  to  rive  orthogonal  measurement  of  the 
undesired  hexanole  field  component.  An  open  house  will 
be  held  at  the  WPI  Van  de  Craaf  Laborator/  (Physics 
building,  basement,  south  end,  Room  004)  on  Saturday, 
April  18,  11-12  noon. 

*F.  W.  Martin,  3rd  International  Conference  on  Ion 
Implantation,  Kingston,  Canada,  Jul>  1930. 


E  2  Optical  Detection  f  "tnrvd  ^7*  Ions  in  an  RF  Trap. 

Vi.  NAGOVRIiEY  a:. a  H .  zi-Z‘.V\.T,  _ r  M-ninvoi.  --  We 

have  observed  t  latruscer.re  from  tne  3“  3*.  tj  j?ri/,  tran- 
8  it  ion  at  27  j  c:  "t'  ;t.~  ;  r.tair.e-i  in  a  srr.jIl.~Rf 
trap  ( 2z  =  .  1  •  *m,  2r  .1-.-).  icr.s  were  ti  .1:  red  in 

a  well  ato:t  13  v.  Its  ty  2-  .icntion  ot  c  ,0  volt 

RF  potential  \t  LI  MHz  tv.-  -.r.-  nn:  ...nd  -nd  caps. 

The  2794JA  s was  a  f re.iu*.*r.cy  cavity  stacil- 

ized  single  frequency  ::r.-.er:ul  dye  laser  producing 
about  1  uW  feejssci  to  1  jT.-  * •: t  diai^it^r .  Viscous 
cooling  -*f  tr.j  10.-.3  in  *  l'  ~  -  vr  ks  buffer  :<m 
ne^essar/  for  observation  jf  tr.e  stror  :vst  sur.als  ,4  x 
10  count,  s).  We  <^ut irate  th_.t  \z>  at  3 G.  lor.s  were 
observed.  At  much  lower  pressures  <■*»  1  C“r"  Torr}  we 
have  observe  1  email  number  a  ^f  icat  (  *  ii.  esttmatv-ii 
for  short  times  ‘an  to  ~.*c  1  .  are  currently  at¬ 
tempting  to  rsliar.vely  :  :cl  the  icr.s  by  down  tuning 
the  laser;  tjr  ultimate  is  to  observe  a  sii.uie  1  * 

ion  cooled  tc  considerably  re low  r ocm.  tenceratuce.  This 
work  is  supported  by  tr.e  JNR. 

E3  Fag  ter.  ;  -y  1  •  r  2  :  r  ♦:  r>  .  "  ?  .  ■.  j  t .  r-yi.  i  n  =  d  .1  j  r.  o.-c- 

r.’-  am.  ii.  iiFMiLr  ar.J  0 .  _ V_.  •  •  t  t :  .  .  ~t .  v. .  -  - 

From  t  3  0  ;o  1  .;»:C  t:  .*  .-ivCir  ..1  11  t  r*  »  «>  rtai*:  :  s 
sufcnectod  to  a  micro  wav?  rj  1  s  ••  »r.  >e  fts.*a;cntv  and  a.i.li- 


F  2  The  Role  of  Excited  State  Polarizabilities  in 
Inelastic  Collisions.*  E.  POLLACK  and  J.  STtVSNS,  •  'n. iv . 
of  Conn.  —  In  A'’"  +  B  A+  +  collisions  the  incii^nc 
potential  at  large  r  is  V*-2ae^/r\  where  x  is  the 
polarizability  of  the  ground  state  of  £.  This  attractive 
potential  crosses  Cculor.b  repulsive  curves  of  final 
exothermic  A+  +  3^  states  and  may  populate  :hym.  I ho 
potential  for  A"^  t  5*  is  given  by  V  ■  -2a'e'*’/ z'* ,  where 
a*  is  the  polarizability  cf  the  excited  state  S*.  At 
infialte  separaticn  V*  lies  above  V  by  the  excitation 
energy  of  B.  Sir..-e  a*  is  much  larger  than  n  (s»e 
ref.  I  as  an  eyerie),  V*  will  cress  V  at  V.ncvn  i  a  tor- 
particle  separation  and  provide  an  intermeoiate  state 
which  crosses  and  may  excite  high  lying  levels  et 
A+  +  B+.  A  similar  model  is  applicable  to  \*^  3j 
collisions.  The  A*  *  R  ■*  A  +  3^  case  may  also  ba  ana¬ 
lyzed  with  a  related  model.  The  effect  of  the  repulsive 
cere  will  be  discu>.>ed .  The  model  will  be  cem.rarei  w:tn 
experimental  results. 

^Supported  by  the  Army  Research  Office  and  the  Univer¬ 


sity  of  Connecticut  Research 
*E.  Pollack,  E.  J.  Robinson, 


Foundation, 
and  B.  Bederson, 


Phys.  Rev.  134,  A1210  (19tw) . 


t i c n  of  Coherent  Mill i'-: 

n-  ' 


r 


relativists  elect. -on  beat  accelerated  ilor.-;  the  axis 
of  v  dielectric  lined  »avc4^;c;  hns  proitieeU  coherent 
Cerer.lo.  r«d.a*iar>  at  ki’.o.att  tower  levels  in  the 
rill-.ncter  rcgijr.  lheor>,  a  Jcscnrt i an  of  me 
experimental  apparatus,  ar.J  sote  experimental  results 
will  be  presor.teu. 

Work  supported  in  part  by  U.S.  Amy  Grant  No.  OAAG- 
29-79-CO  2030. 


F 4  Relativistic  Pear.  Propagation  in  a  Sielectric 
Lined  Navv .:i  ~c  ■  .).  ER-V-dCL"..  J.S.  ■<>  .5:1.  '.IVAN. 

Dartm.i  .~6c’.  U  Har.o.-e-.  N.H.  0^753— Relativistic 
electron  roans  propagating  along  the  axis  of  a 
diel  ectric  lined  cylindrical  waveguide  have  been  shown 
to  proouce  coherent  till  meter  radiation,  i  Surface 
charge  collecting  on  the  dielectric  liner  causes 
problems  of  bean  dynamics  uhxcn  could  Unit  the 
practical  uses  of  this  process.  Experiments  demonst¬ 
rating  the  effects  of  charge  collection  and  some 
possible  solutions  will  be  described. 

1.  It.  Fetch,  K.  Busby,  R.  Luvr.an,  and  J,  Walsh,  , 
Bull.  An.  Phys.  Soc.  24,  1076  (1979). 

Work  Supported  in  part  by  U.S.  Army  Grant  No.  DAAG-  I 
29-79-CO  2030.  _ I 


F5  Cerenhov-Cyclotron  Instability.  J.  GOLUB  and 
J.F.  KALttTT  Uattnouth  >~  1 U  ;e .  Hanover,  \.H.  037S3-- 
It  has  over  shown  that  a  slow  cyclotron  wave  oronaga- 
ting  along  an  election  beam,  both  in  a  dialectric 
medium,  is  unstable. 1  A  slow  space  charge  wave 
propagating  at  the  Cerenkov  angle  in  a  similar  system 
is  also  unstable.  We  have  investigated  the  off-angle 
propagation  of  a  slow  cyclotron  wave.  The  system  is 
shown  to  be  unstable  at  a  wavelength  given  by 

*  *  (Sn-l)/(.7c/2rc) 

where  8  =  v/c  is  the  relative  beam  velocity;  n  is 
the  index  of  refraction  of  the  dielectric;  and  S5C  i? 
the  Telativistic  cyclotron  frequency.  A  small  signal 
temporal  growth  rate  has  been  derived  for  propagation 
at  a  general  angle.  At  the  Cerenkov  angle,  using 
realistic  experimental  parameters  and  1  =  SOOum,  this 
growth  rate  is  approximately  or.e  order  of  magnitude 
greater  than  that  of  the  on-angle  case.  Potential 
application  of  this  instability  to  some  practical 
free  electron  laser  systems  will  be  discussed. 

1.  D.  Kaoilow  and  J.E.  Walsh,  Bull.  Am.  Phys.  Soc., 
25,  6  (19S0) . 

Supported  in  part  bv  The  Office  of  Naval  Research 
Grant  #SC0014-73-C-0760. 


P  SESSION  G.  GRAVITATION/COSMOLOGY 
t  Saturday  morning.  18  April  1931 
k  Olm  Hail,  Room  223  at  9:00  A.M. 
k  H.  W.  Hibinjer,  presiding 


G  1  Mechanism  of  Elcctr'T.aguetic  Radiation , 

H.  HI'.  S ISClir. •  V  rtestar  Fc:vte:-.ric  Institute.— 
Mach's  principle  avserts  that  inertial  effects  are 
caused  r.of  by  acceleration  relative  to  seme  reference 
frame  tut  by  a:r<  deration  relative  tr>  th«  rest  if  the 
mass  of  the  universe;  however,  experimental  confirmation 
of  this  prlnci;le  Is  preiiuded  largcl”  by  the  i  ipossl- 
bi.'.lty  ci  removing  that  ambient  -c3s  or  of  si.leldine  a 
system  from  it;,  effect.  The  correspc.iuir.g  elercrcbiag- 
netlc  principle  would  attribute  r.idlaticr  bv  a  charge 
to  acceleration  relative  to  the  ambient  charge.  An 
experiment  will  be  described  that  can,  in  principle,  ■ 
investigate  the  mechanlss  of  electromagnetic  radlaticr.: 
does  a  cl  it.  rad  btr.aws-  it  has  '..ojutc  *•?..’ er- 
atior.  cr  bar  .:  -.«■  r  >.  ,n  a  •  rlerp.t  ion  relat'.v'  to 
other  charges.  The  possible  results  of  such  an 


experiment  and  their  interpretation  will  also  be 
discussed. 

Submitted  by  L.  R.  rtaa-Mohan 


C  :  A  Me •■•■»)  Tl'-w  f'  ;-■>  Gcl-r  Svstca  - 

Fart  I  ■  nASGli.  f.  clc.l.'f,  '  1 1  r^d .  The 
solar  svstea  has  a  r.ew  muter  2  .  • .  - 62.  with 
dimensions  cf  tas/sec  in  tr.e  plan*  cf  the 
ecllptls.  The  Newtonian  relation  for  orbital 
bodies  is  a  limiting  fort  cf  the  expression 
derived .  Juuiter  .  as  nr.  influential  role.  A 
nor.-llnear  relation  is  quantized  where 
p%  15  a  cci7»tant.  With  V  I?.  as  the 

r  i.  »-<  constant,  where  K  -as  the  value 
unity  and  taxes  care  of  the  Users  lens,  has 
the  value  of  0,01?26£  a.  V-  The  number  h 
taxes  Integral  values  while  P  Is  the  period' 
in  sidereal  tropical  years.  The  value  of  wi 
for  Venus  and  Earth  are  72  and  57  v 

respectfully.  Values  of  W  for  such  diverse 
bodies  as  "Object  Kowal",  Halley's  coset  and 
Ceres  are  1367*  17*9  ar.i  26?  .  The  number 

gives  the  average  sesl-sajor  axes. 

1  *  The  secular  variations  for  the 

outer  planets  are  discussed. 


G  3  Quarks  and  Particles.  HAH  OLD  P. 

SCHVEDE,  Retired .  Kinematic  considerations 
of  the  sasses  of  particles  show  Intriguing 
results.  Phase  pathe  for  Plons  and  Xaons 
are  indicated.  A  derived  number  0.611004091 
is  considered  as  the  mass  of  the  electron.ln 
Mev/  c2  and  a  number  of  about  100  ey/  oz  la 
obtained  which  Is  perhaps  related  to  the 
electron  neutrino.  Relations  are  given  for 
the  proton  and  the  neutron.  The  twin  prime* 
centered  at  1?8  arid  1020  play  a  role. 

Rational  fractions  same  as  the  alleged 
fractionally  charged  quarks  appear  naturally. 
Elliptic  functions,  particularly  the 
kfelers trass  relations,  are  revealing.  A 
fundamental  relation  is  easily  eultable  tt> 
cover  nany  types  of  foroes. 


SESSION  H:  GENERAL 
Saturday  morning.  18  April  1981 
Oiin  Hall,  Room  12G  at  9  . 00  A.M. 
P.  Gl3nz,  presiding 


H 1  Association  of  Hemoglobin*  ns  Studied  bv  Intensi¬ 
ty  Fluctuation  .' s?c*.rof>' orv.  t.  J.  L-T-TT.'T* .  "air.  cf 
vcnneclitul  —  \  '  »  :or»al  acult  he'.vitoin  tK'oA)  'exists 
in  vivo  an  a  teiramer-j  two  chain?  cf  whi*h  »me  of  theaC 
form  anti  two  chain*  cf  which.  »rt /9  .  The  extent  of  the 
aeeociation  of  these  four  chame  jr.to  tetremer,  in  vitro, 
wae  mearurei  bv  mtenritv  fluctuation  eoectroecopr  (IPS) 
and  a  pH  dependent  readies  equilibrium  conatant  deduc¬ 
ed.  Ferultr  are  compared  with  valuer  obtained  bv  other 
methcls.  (2)  formation  of  multi-teirameric  aggreiratea  of 
HbA  la  found  to  occur,  in  vitro,  at  low  ieuic  strength. 
The  averrre  eixe  of  an  a.-rrrgate  end  the  diepereion  lo 
eir.ee  war  determined  bv  ITT  .  Ir.iicatiena  of  lie- tel  ag¬ 
gregation  at  verv  him.  ionic  atrer.cn.-.  wae  e.lfo  ebrervei. 
(})  Apolicaticn  of  I.-f  to  etudirr  of  so-called  eitkHng 
hemof.obinr  (Kbf)  will  be  diecuared.  T.-.e  formation  of 
the  r-l  rv«re  of  deoTv-ubf  mav  be  preceded  bv  etatee  of 
limited  segregation.  Tr.ese  etatee  ehruld  be  amenable  to 
etudv  bv  IF  3. 

*Y.  J.  LaGattuta,  etal.,  Bicpnve.  J. 


HI  HarS-He1;'  Ca '  *ji  ato't  <n  ?J>e  Aoa’.'tI1;  cf 

SjtcOt  .'.''T-M'.-iV  •  T-v-i;'  r-~«tc 

)r ji. --cy-.ent  jd.  .-.c;  ‘r  sfica:e  "sir..;?.;  -ivi  tne 
concept  of  prcje.ticr  ratrices,  le:d  to  the  deve1c;»ient 
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Abstract  Submitted 
For  the  Twenty-third  Annual  Meeting 
Division  of  Plasma  Physics 
October  12  to  16,  1981 

Category  Number  and  Subject  Microwave  Generation  -  4.8 
CD  Theory  0  Experiment 

Power  CharacterizaCion  of  a  Cerenkov  Maser*. 
S.  VON  LAVEN,  J.  BRAN SCUM,  J.  GOLUB,  R.  LAYMAN , 

D.  SPEER,  D.  WISE,  Dartmouth  College  — Stimulated 
Cerenkov  radiation  is  generated  by  a  mildly  relativ¬ 
istic  electron  beam  in  a  waveguide  which  is  partially 
filled  with  dielectric.  Substantial  levels  of  output 
power  have  been  generated  at  wavelengths  varying 
between  1.5  and  10  mm.  The  beam  energies  range  from 
100  to  250  keV,  and  the  currents  from  15  to  20A. 

Pulse  widths  are  of  the  order  of  lp  sec.  Beam  power 
to  microwave  power  conversion  efficiencies  of  1  to 
5%  have  been  achieved  to  date.  At  the  present  time 
the  system  runs  as  a  self  oscillator  with  feedbacK 
provided  by  the  relatively  small  reflection  from  the 
end  of  the  tubular  quartz  guide.  The  overall  cavity 
Q  in  this  case  is  low  and  oscillation  commences  only 
at  relatively  high  gain  levels.  In  order  to  better 
control  the  feedback,  a  number  of  quasi-optical 
cavity  and  output  coupling  schemes  have  been  tested. 
The  effects  of  these  on  output  power  and  wavelength 
will  be  presented  and  discussed. 

*Work  supported  by  ARO  Contract  number  DAAG29-79-C- 
0203. 
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Cerenkov  Masers:  A  Possible  Plasma  Heating  Source 

J.  Walsh,  J.  Branscum,  J.  Golub,  R.W.  Layman, 

D.  Speer,  S.  Von  Laven 

Dartmouth  College 
Hanover ,  N .  H . 

A  Serenkov  maser  consists  of  an  electron  beam, 
a  dielectric  resonator  and  output  coupling  optics.  The 
beam  velocity  can  exceed  the  phase  velocity  of  the  wave 
in  this  system,  and  when  it  does,  a  coherent  instability 
leads  to  beam  bunching  and  a  transfer  of  energy  to  the 
wave.  The  field  in  the  beam  channel  is  also  evanescent. 

The  decay  rate,  however,  scales  as  k/y  where  k  is  the 
axial  wavenumber  of  the  wave  and  y  is  the  ratio  of  the 
energy  of  the  electron  and  its  rest  mass.  Hence  by 
using  mildly  relativistic  electron  beams  (y  =  1. 1-1.6) 
good  beam-to-wave  coupling  can  be  obtained  in  the  lower 
mm  part  of  the  spectrum.  Depending  upon  their  complexity 
and  ultimate  performance  characteristics,  devices  of  this 
kind  may  have  a  number  of  applications  in  plasma 
diagnostics  and  heating. 

In  order  to  test  the  basic  ideas  underlying  such 
devices,  a  high-voltage  (400  Kv  max.)  pulse  transformer- 
based  e-beam  generator  has  been  used  to  drive  tubular 
quartz  resonators.  At  the  present  time,  coherent  output 
has  been  obtained  over  the  range  10mm-1.5mm.  A  summary 
of  theoretical  expectations  and  recent  experimental  results 
will  be  presented. 


Work  supported  in  part  by:  AFOSR  Grant  r77-3410B,  ARO 
Grant  *DAAG-29-79-C-0203  and  OMP.  Grant  =M00-14-79-C-0760 . 


Jj.  C.  Phillips,  Covalent  Bonding  In  Crystals.  Molecules, 
and  Polymers  (U,  of  Chicago  Press,  1969);  H.  Watanabe  as 
cited  in  ref.  1. 
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D 3  Nonlinear  screening  of  negative  point  charges  In 
diamond,  silicon,  and  germanium.  P.  CSAv i T.SZKY  and  K.R. 
BSSCRSTOnT  t'niv.  of  Mair7e.--Cornolti  and  Resta1  (CR )  have 
recently  forniuTated  a  Thonas-Fermi  (Tf)  theory  of  nonlinea 
Impurity  screening  in  semiconductors.  CR  have  obtained  th 
spatially-variable  dielectric  function  for  charges  i  le0, 
t  4e0  (eo  the  proton  charge)  embedded  in  pure  diamond. 
Si  and  Ge.  The  nonlinear  results  differ  importantly  from 
the  results  of  the  linearized  TF  theory2.  tie  have  pre¬ 
viously  solved  the  nonlinear  TF  equation  of  CR  by  an  equiv 
alent  variational  principle3.  We  have  used  a  two-paranete 
trial  solution  and  considered  the  cases  of  +le0,  +2e0,  +  3e 
and  +4e«  in  pure  diamond  Si  and  Ge.  We  have  now  extended 
our  variational  approach  to  the  negative  charges  -le0,  -2e 
-3eo,  and  -4eo  in  the  above  semiconductors'*.  Our  analytic 
results,  using  again  a  two-parameter  trial  solution,  agree 
remarkably  well  with  the  numerical  results  of  CR. 

*F.  Comolti  and  R.  Resta,  Phys.  Rev.  B  17,  3239  (1978). 

2R.  Resta,  Phys.  Rev.  B  16,  2717  (1977). 

3P.  Csavinszky  and  K.R.  Brownsteln,  Phys.  Rev.  B  (to  be 
published). 

*P.  Csavinszky  and  K.R.  Brownstein,  Phys.  Rev.  8  (to  be  l 
published) . 


D 4  Phonon  Conduction  lit  Elastically  Anisotropic 
Cubic  Crystals.  A.K.  McCurdy,  Worcester  Polytechnic 
Institute. —  Striking  differences  In  the  boundary-scat¬ 
tered  phonon  conductivity  are  predicted  along  the  prin¬ 
cipal  axes  of  cubic  crystals.  The  results  are  due  to 
phonon  focusing  arising  from  elastic  anisotropy.  Nor¬ 
malized  curves  of  phor.on  conduc-'vlty  have  been  calcu¬ 
lated  for  samples  of  square  cr  .  action  as  a  function 
of  the  elastic  anisotropy  A  -  iC.44/  (Cn  -Cj2)  and  the 
elastic  ratio  Anisotropies  of  more  than  50X 

are  possible  for  different  rod  axes.  Silicon  and  calci¬ 
um  fluoride,  materials  in  which  this  anisotropy  was 
first  reported,  are  shown  to  be  very  favorable  materi¬ 
als  to  demonstrate  this  anisotropy.  For  silicon  and 
calcium  fluoride  samples  of  rectangular  crcss-sectlou 
the  thermal  conduction  Is  shown  to  depend  upon  the  cry¬ 
stallographic  orientation  and  width  ratio  of  the  side 
faces  for  samples  with  the  same  <110>  rod  axis.  Re¬ 
sults  are  expressed  In  a  convenient  form  for  predicting 
the  phonon  conductivity  of  elastically  anisotropic  cry¬ 
stals  given  the  linear  dimensions,  the  density  and  the 
elastic  constants. 


D S  Equilibrium  Configuration  of  an  Ethylene  Monolayer 
Adsorbed  on  Graphite  by  Eric  Ehrhardt,  tarry  Pratt,  Howard 
Patterson3  (University  of  Maine,  Orono,  Maine  04469)  and 
Larry  Passell  (Brookhaven  National  Laboratory,  Upton,  N.Y. 
11973). 

In  this  talk  we  will  describe  computer  calculations  • 
which  have  been  carried  oat  to  determine  the  equilibrium 
configuration  of  an  ethylene  moniayer  phys-adsorbed  on 
graphite.  A  four  sublattice  structure  was  assumed  from  thf 
results  of  elastic  neutron  scattering  studies.  In  the  cal¬ 
culation  we  have  included  ethylene-ethylene  interaction  as 
well  as  the  ethylene-graphite  substrate  interactions.  Our 
results  for  ethylene  adsorbed  on  grapnite  are  very  similar 
to  those  of  Fusilier,  Gillis,  and  Raich^  for  nitrogen  odsoi 
on  graphite.  That  is,  the  ethylene  molecules  show  a  herrii 
bone  pattern  with  the  ethylene  C=C  axis  almost  perpendlcul; 
to  the  graphite  basal  plane. 

1.  Research  partially  supported  by  NSF,  Department  of 
Materials  Research  DMR  77-07140. 

2.  C.  R.  Fuselier,  N. '$.  Gillis  and  J.  0.  Raich,  Solid 
State  Communications  25,  747  (1978). 


E2  Cerenkov  Radiation  as  a  Source  of 

Millimeter  Radiation.  J.  b  RANSOM,  J.  WALSH* * - - 
An  electron  moving  with  a  velocity  greater  than  tho 
phase  velocity  of  an  electromagnetic  wave  produces 
radiation.  This  situation  can  be  achieved  physically 
in  a  number  of  ways,  one  of  which  is  to  have  the 
electron  move  through  (or  near)  a  dielectric.  The 
resultant  radiation  is  known  as  Cerenkov  radiation. 

This  phenomenon  may  be  used  as  the  basis  of  a  practical 
millimeter-infrared  radiation  source.  Three  problems 
must  be  considered:  1)  method  for  making  o/ck  <  1. 

2)  dependence  of  beam-wave  coupling  on  k,  y,  and 
position.  3)  method  of  insuring  that  energy  moves 
into  the  wave  from  the  beam.  The  first  problem  is 
solved  by  using  a  dielectric  resonator  with  internal 
reflection  types  of  inodes.  The  second  is  approached 
by  examining  the  phase-velocity  dependence  of  the 
strength  of  the  electric  fielJ  which  is  undergoing 
total  internal  reflection.  The  third  problem  is 
approached  by  constructing  a  dielectric  Fabry-Perot 
resonator.  Each  component  will  be  discussed  briefly. 

‘Supported  in  part  by  Office  of  Naval  Research  Grant 
t  MDO0-16-79-C-0760. 


E  3  Millimeter  Wave  Generation  with  a  Relativistic 
Electron  Beam.  R.W.  LAYMAN,  J.  BRANS CUM,  and  J.  WALSH* 
— Production  of  electromagnetic  radiation  in  the  3S  to 
ISO  GHz  range  by  a  mildly  relativistic  electron  beam 
accelerated  along  the  axis  of  a  dielectrically  lined 
cylindrical  waveguide  has  been  reported  elsehwere3 . 

This  process  shows  potential  as  a  tunable  source  of 
high  power  millimeter  radiation.  Results  of  experi¬ 
mental  work  to  demonstrate  this  possibility  will  be 
described. 

1.  K.  Felch,  K.  Busby,  R.  Layman,  and  3.  Walsh,  Bull. 
Am.  Phys.  Soc.  24,  1076  (1981) . 

•Work  supported  in  part  by  U.S.  Army  Research  Office 
Grant  •  DAAG-29-79-C-0203. 


E 4  Evidence  of  Bifurcation  Universality  for  First 
Sound  Subbarmonic  Generation  in  Superfluid  He*.* 

C.W.  SMITH,  D.A.  FARRIS,  and  M.J.  TEJn'ANI ,  University  of 
Maine,  Orono. --Measurement  of  the  subharmonTc  responses 
of  superfluid  helium-4  to  ultrasound  at  3  MHz  is  report¬ 
ed.  A  matched  pair  of  P2T4  thickness  mode  transducers 
are  positioned  parallel  on  a  common  axis  in  an  open 
geometry.  One  transducer  is  employed  as  a  first  sound 
source  and  the  other  as  a  receiver.  The  received  signal 
Is  Fourier  analyzed.  Several  subharmon*c  frequencies 
(f0/n,  where  n  *  2,3,4,..)  of  the  applied  frequency  f0 
are  observed  above  specific  sound  thresholds.  Prelim¬ 
inary  results  for  n  1  2,4,8  have  been  analyzed  in  terns 
of  the  bifurcation  theory  for  a  nonlinear  system  in 
transition  to  chaos'.  To  within  experimental  error  the 
thresholds  for  the  onset  of  the  subharmpnics  aqree  with 
the  theoretically  predicted  value  of  the  universal 
geometric  convergence  constant,  4.7.  Comparison  of  the 
observed  decrease  In  the  amplitude  of  successive  sub- 
harmonics  with  theoretical  prediction,  the  sequence 
n  =  3,6,12  and  apparent  pnase-locking  behavior  are 
currently  under  investigation. 

♦Supported  by  NSF,  DMR800S358  and  AF0SR,  NP  80-151. 

'M.J.  Fergenbaum,  J.  Stat.  Phys.  19,  2b  (1978). 


E  5  Photoacoustic  Studio  of  todlm  Vapor.  A.  BRAN,  M. 
SARIASCUE  and  B.  Dl  TA-.V-. O,  Boston  College.  Ve  con¬ 
structed  a  photoacoustic  apparatus,  using  an  acoustical 
cylindrical  cavity  operating  in  a  longitudinal  rude  and 
used  molecular  iodine  vapor  as  specimen  and  Argon  as 
buffer  gas  to  study  the  phntuacoustic  characteristics  or 
the  system.  The  ledlr.e  rolcculos,  excited  periodically 
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A  Two-Stage  Cerenkov  Oscillator,  J.E.  WALSH, 

S.  VON  LAVEN,  R.W.  LAYMAN,  E.  GARATE  and  R.  COOK, 
Dartmouth  College* — Recently^-,  100  KW  levels  of  output 
power  have  been  obtained  in  Cerenkov  oscillator  experi¬ 
ments.  The  device,  which  has  been  described  in  more 
detail  elsewhere,  consists  of  a  pulse  transformer¬ 
generated,  100-200  KV,  10-20  A,  1  usee  beam  pulse,  a 
dielectric-loaded  guide,  and  a  variety  of  output 
coupling  structures.  The  output  power  levels  indicate 
that  the  beam  is  nearly  100%  modulated  during  its 
passage  through  a  10  cm-long  resonator.  This  in  turn 
suggests  that  enhancement  of  the  output  power  could 
be  obtained  if  the  now  bunched  beam  is  passed  through 
a  second  resonator  which  is  weakly  coupled  to  the  first, 
and  whose  phase  velocity  is  lower  than  that  of  the 
first  resonator.  Such  a  device  is  similar  to  a  multi¬ 
cavity  klystron  or  a  velocity-tapered  traveling  wave 
tube.  A  theoretical  analysis  and  results  of  experi¬ 
mental  investigation  of  this  concept  will  be  presented. 

1.  A  High  Power  Cerenkov  Maser  Oscillator,  S.  Von 
Laven,  J.  Branscum,  J.  Golub,  R.  Layman  and  J.  Walsh, 
Appl.  Phys.  Lett.,  Sept.  1,  1982,  to  be  published. 

Supported  in  part  by  AFOSR  Contract  ?>  82-0168  and  ARO 
Contract  #  DAAG-29-79-C-0203-A01. 
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V58 ,  Edited  by  C.  Marton  (Academic  Press,  1982) 
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INTRODUCTION 


i  a.  Cerenkov  radiation 

The  electromagnetic  wave  produced  by  a  charged 

particle  moving  with  greater  than  light  velocity  in  a 

dielectric  medium  is  known  universally  as  Cerenkov1 

radiation.  Cerenkov's  experiments,  which  were  performed 

independently  during  the  1930 's,  and  the  subsequent  analysis 

2 

of  the  phenomena  by  Frank  and  Tamm  did,  however,  have  some 
precursors. 

3 

Heaviside,  m  1889,  analysed  the  problem  of  the 
radiation  produced  by  a  charged  particle  when  it  moved  with 
uniform  velocity.  This  work  was  done  prior  to  the  develop¬ 
ment  of  the  special  theory  of  relativity,  and  Heaviside 
assumed  that  it  was  possible  for  a  particle  to  move  with  a 
velocity  greater  than  that  of  light  in  a  vacuum.  When  it 
was  so  assumed,  radiation  was  produced.  In  a  formal  sense, 
his  results  were  similar  to  those  of  Frank  and  Tamm. 

4 

Sommerfeld,,  in  1904,  without  apparent  knowledge  of  Heavi¬ 
side's  results,  performed  a  similar  analysis.  There  were 
also  some  experimental  precursors  to  Cerenkov's  work.  M. 

5 

Curie,  m  1911,  observed  that  radiation  produced  in  the 
walls  of  glass  containers  holding  radioactive  material  was 
probably  due  in  part  to  the  penetration  of  the  glass  by  fast 
charged  particles.  Some  experiments  performed  by  Mallet^  in 


in  1926  were,  in  part,  observations  of  Cerenkov  radiation. 

None  of  this  early  work,  however,  lessens  the  importance  of 
pioneering  experiments  of  P.A.  Cerenkov. 

Following  the  initial  experiments  of  Cerenkov  and 

theory  of  Frank  and  Tamm,  an  extremely  large  number  of  both 

theoretical  and  experimental  contributions  have  appeared. 

General  discussions,  with  hundreds  of  additional  references, 

7  8 

may  be  found  in  Jelley,  in  Zrelov  and  in  the  review  article 

.  q 

by  Bolotovski .  The  interest  of  many  contributors  has  been 
the  potential  use  of  the  Cerenkov  process  as  a  practical 
radiation  source.  Notable  among  these  contributions  were  the 
papers  of  Ginzburg,10  in  which  he  considered  a  number  of  ways 
in  which  electrons  could  be  coupled  to  dielectrics  and  be  made 
to  produce  radiation  in  the  millimeter  and  submillimeter 
regions  of  the  electromagnetic  spectrum. 

Much  of  the  early  work  dealt  with  the  radiation 
produced  by  single  electrons.  As  we  shall  see,  however,  this 
spontaneous  radiation  is  a  relatively  weak  process  for  all 
wavelengths  longer  than  that  of  the  blue  ultra-violet  regions 
of  the  spectrum.  Hence,  in  order  to  produce  useful  amounts 
of  radiation,  it  was  natural  to  consider  the  radiation 
produced  by  a  bunched  electron  beam.  At  wavelengths  long 
compared  to  the  length  of  the  bunch,  the  radiated  power  is 
proportional  to  the  square  of  the  number  of  electrons  involved 
and  hence  the  power  emitted  rises  dramatically.  A  number  of 
experiments  were  designed  to  explore  the  properties  of  the 
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Cerenkov  radiation  produced  by  prebunched  electron  beams 

moving  in  close  proximity  to  a  dielectric  surface. 

Important  contributions  were  made  by  Coleman,11  by  Danos,12 
13  14 

by  Lashinsky  ,  and  by  Ulrich.  In  these  experiments, 
no  provision  was  made  for  feeding  back  the  emitted  radia¬ 
tion  on  subsequent  bunches  and  hence  they  could  be  cate¬ 
gorized  as  observation  of  enhanced  spontaneous  emission. 

Suggestions  have  also  been  made  that  Cerenkov 
radiation  could  be  used  as  the  basis  of  a  microwave 
tube  . 15,16,17  In  these,  a  dielectric  tube  was  used  as  a 
slow  wave  structure.  The  general  configuration  suggested 
was  similar  to  that  used  in  traveling  wave  tubes.  When 
electron  beams  in  the  energy  and  current  range  found  in 
conventional  microwave  tubes  are  used,  however,  the 
resulting  devices  are  unsatisfactory  for  several  reaons. 

We  will  develop  this  line  of  argument  carefully  in  subse¬ 
quent  sections,  since  these  difficulties  must  be  surmounted 
in  constructing  a  useful  Cerenkov  source. 

A  major  difficulty  in  constructing  a  Cerenkov  source 
that  is  capable  of  producing  useful  amounts  of  radiated 
power  is  the  coupling  of  the  electron  beam  to  the  dielectric. 
In  elementary  discussions,  it  is  usually  assumed  that  the 
electron  is  passing  right  through  the  dielectric.  This  can 
actually  be  done  for  the  limiting  case  of  very  high  energy 
particles  and  gaseous  or  liquid  dielectrics.  In  this 
regime,  Cerenkov  radiation  actually  finds  wide  practical 
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application  as  a  diagnostic  tool.  There  have  also  been 

18  19  20 

serious  attempts  '  '  to  observe  stimulated  Cerenkov 

radiation  in  the  visible  and  ultra-violet  region  from  a 

high-energy  beam/gaseous  dielectric  combination.  In  these 

latter  experiments,  momentum  modulation  '  by  an  applied 

electromagnetic  has  been  observed,  but  as  yet  there  is  no 

clear-cut  evidence  of  true  stimulated  emission.  An 

alternative  to  passing  an  electron  beam  directly  through  a 

dielectric  is  to  let  a  beam  propagate  along  a  channel. 

21  22  23 

Recent  experiments  '  '  in  which  millimeter-wavelength 

stimulated  Cerenkov  radiation  has  been  observed  have  been 
of  this  type. 

A  primary  purpose  of  the  present  paper  is  to  explore 
the  potential  of  the  latter  option.  We  will  establish 
criteria  necessary  for  producing  usable  levels  of  stimulated 
Cerenkov  radiation  at  wavelengths  which  are  short  compared 
to  the  characteristic  scale  length  of  both  the  transverse 
and  longitudinal  dimensions  of  a  dielectric  resonator. 


i  b.  Cerenkov  masers 


The  goal  of  the  general  area  of  research  pertaining 
to  devices  now  often  called  free-electron  lasers  iS  to 
produce  coherent,  tunable,  moderate,  and  high-power 
radiation  in  parts  of  the  electromagnetic  spectrum  where 
such  a  source  is  not  now  available.  All  of  the  devices 
suggested  to  date  have  much  in  common  with  microwave  tubes, 
and  hence  the  designation  "maser"  or  "laser"  could  be  the 
subject  of  debate.  It  is  possible,  but  not  necessary,  to 
formulate  the  equations  of  motion  quantum-mechanically . 

The  electron  transitions  are  between  continuum  states.  The 
recoil  due  to  single-photon  emission  is  negligible,  and  thus 
Planck's  constant  does  not  appear  in  any  final  working 
formula.  A  classical  analysis  based  on  either  fluid  or 
kinetic  equations  will  lead  to  the  same  expressions.. 
Therefore,  much  of  what  is  know  about  microwave  tubes  will 
apply  also  to  free-electron  lasers.  Microwave  tubes, 
however,  operate  at  wavelengths  comparable  to  or  greater 
than  the  device,  while  the  opposite  will  be  the  case  for  any 
free-electron  laser  or  maser.  This  difference,  although 
minor  from  some  viewpoints,  accounts  for  many  of  the  diffi¬ 
culties  encountered  in  attempting  to  build  short  wavelength, 
beam-driven  radiation  sources. 
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II.  THEORY 


A  series  of  calculations  aimed  at  establishing  the 
beam  energy,  current,  and  velocity  spread,  which  are 
required  in  order  to  obtain  growth  of  stimlated  Cerenkov 
radiation,  will  be  presented  in  this  section.  The  analysis 
will  proceed  along  classical  lines  similar  to  those  used  in 
traveling  wave  tube  and  beam  plasma  theory.  In  Sections  II  A 
and  B  we  will  examine  the  exponential  gain  of  stimulated 
Cerenkov  radiation  obtained  when  it  is  assumed  that  either 
a  strongly  magnetized  or  a  completely  unmagnetized  mono- 
energetic  electron  beam  passes  directly  through  a  dielectric 
medium.  The  limit  implied  by  the  assumption  that  the  beam 
is  monoenergetic  will  be  examined  in  Section  II  C,  and 
modified  gain  formulas  will  be  derived.  Section  II  D  will 
then  be  devoted  to  some  resonator  configurations  which  are 
more  practical  for  the  present  application.  Emphasis  will 
be  on  the  slab  geometry,  since  in  this  case  it  is  possible 
to  present  a  reasonably  compact  analytic  result.  The 
results  obtained  from  other  geometries  will  be  similar.  A 
few  brief  comments  and  calculations  related  to  nonlinear 
effects  will  be  outlined  in  Section  II  E. 


ii  a.  Cerenkov  gain  on  a  strongly  magnetized  beam 

We  consider  first  the  case  of  a  plane  wave  propagating 
at  an  angle  to  strongly  magnetized  electron  beam.  The 
geometry  is  shown  in  Fig.  3. 

II  A.  1.  Current  Modulation 

When  the  beam  is  strongly  magnetized,  the  beam  density 

and  modulation  are  one-dimensional  and  lie  along  the  beam 

and  magnetic  axis.  In  this  limit,  the  linearized  equation 

for  the  velocity  modulation  has  only  one  component,  v  , 

z 

where : 


dv 
_ z 

dt 


(1) 


The  solution  of  this  equation  for  the  assumed  Ez  is  readily 
found : 


v 


z 
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U)-k  vQ 


E 
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(2) 


This  result  together  with  a  linearized  equation  of  continuity 
gives  for  the  density  modulation  n  , 
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[in_e/ (my3)  ]  kEz/ (u)-kvQ)  2 

(3b) 

Thus  the  current  produced  by  the  wave  is  given  by: 


J  =  -  n  ev  -  nev 
20  o 


2  2 

where  w  =  4nn  e  /m  is  the  beam  plasma  frequency, 
p  o 


(4a) 

(4b) 


II  A.  2.  The  Wave  Equation 


In  writing  the  second  of  these,  we  assume  that  the  wave  and 
the  beam  are  in  a  dielectric  medium  where: 

D  =  €E  (6) 

Taking  the  time  derivative  of  the  second  of  Maxwell's 
equations  and  substituting  the  first  gives  a  single-wave 
equation: 


V  x  V  x  E  +  ~ 
c 


3£ 
3 1 


(7) 


There  is  no  current  component  in  the  direction  perpendi¬ 
cular  to  the  beam  and,  hence,  the  perpendicular  component 

of  eq.  (7)  may  be  used  to  express  E  in  the  terms  of  E  . 

x  z 

Doing  this  and  substituting  into  the  longitudinal  component 
of  eq.  (7),  and  making  use  of  the  assumed  time  and  z-depen- 
dence,  we  obtain  a  single  wave  equation  for  E_: 

Z 


Since  we  have  also  assumed  a  plane  wave  dependence  in  the 
perpendicular  as  well  as  longitudinal  direction,  we  also 
obtain  immediately: 


co  2 


(“4  .  k2)i 


(co-kvQ) 


0 


(9) 


where  p  is  the  perpendicular  component  of  the  wave  number 


II  A.  3.  The  Dispersion  Relation 

We  are  obviously  interested  in  the  case  Ez  ?  0  and, 
hence,  the  coefficient  of  eq.  (9)  is  the  dispersion 
relation: 


for  a  plane  wave  propagating  through  a  dielectric  medium 
at  an  angle  to  a  strongly  magnetized  electron  beam. 

Equation  (10)  is  a  quartic  in  both  w  and  k  and 

hence  it  has  four  roots.  When  v  <  c/Ve  ,  all  four  roots 

o 

are  real,  while  if  vq  >  c/Ve  it  has  two  real  roots  and  a 
complex  conjugate  pair.  One  of  the  real  roots  is  related 
to  a  wave  propagating  in  the  direction  opposite  to  that 
of  the  beam  (in  the  negative,  z-direction) .  The  other 
three  result  from  the  coupling  of  an  electromagnetic 
wave  propagating  in  the  positive  z-direction  and  two  beam 
space  charge  waves.  The  latter  two,  fast  and  slow  space 
charge  waves,  would  be  normal  modes  of  the  free  beam.  In 
the  presence  of  the  dielectric,  however,  they  become 
coupled  to  the  electromagnetic  wave.  When  the  velocity 
threshold  vQ/c  =  l/VT  is  exceeded,  the  beam-wave  dielectric 
system  becomes  unstable. 


II  A.  4.  Cerenkov  Gain 


The  presence  of  the  beam  is  obviously  felt  most 
strongly  for  waves  near  "synchronism",  i.e.  when 

u  r  kv  (11a) 

o 

and 

kv  ~  w  ,  (lib) 

o  ok 

2 

where  we  define  w  ,  as 

ok 

wk2  »  c2 (k2+p2) / e  ,  (11c) 

the  dispersion  relation  of  the  electromagnetic  waves  in 
the  absence  of  the  electron  beam. 

In  the  region  where  eqs.  (11a)  and  (lib)  are  valid, 
the  dispersion  relation,  eq.  (10) ,  becomes  an  approximate 
cubic : 

2 

(uj-kv  ) 3  -  u)(l-l/e2e)  =  0  (12) 

•  °  2yJe 

Equation  (12)  follows  from  (10)  when  kvQ  is  set  equal 

to  a)  in  those  terms  where  the  substitution  does  not  give 

2 

zero.  This  is  a  valid  assumption  provided  u)p  is  small 
in  a  sense  which  we  will  define  shortly. 
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The  root  corresponding  to  eq.  (13b)  is  an  exponentially 
growing  wave  in  either  time,  Imu  ^  0,  or  space,  Imk  ^  0. 
The  choice  between  these  will  be  determined  by  initial 
and  boundary  conditions. 

We  will,  for  the  moment,  assume  that  the  spatial 
growth  is  of  interest  and  we  will  let  Imk  =  a  ,  then: 

a<, 

Examination  of  eq.  (14)  shows  that  the  spatial  gain 
increases  with  the  two-thirds  power  of  the  beam  density 
and  the  one-third  power  of  the  frequency.  It  vanishes 
as  the  beam  energy  approaches  the  Cerenkov  threshold 
and  decreases  as  e  and  y  become  large. 

Shown  in  Fig.  4  are  sketches  of  free  wave  dispersion 
curves  for  two  different  perpendicular  wave  numbers,  p^ 
and  P2  .  The  curves  leave  the  k  =  0  axis  at  the  point 
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w/c  =  p//T  ,  cross  the  speed  of  light  at  oj/c  =  p/Ve-1  and 
then  asymptotically  approach  a  wave  prpagating  in  the 
z -direction.  Along  this  curve,  the  angle  of  propagation 
is  varying  from  0  =  tt/2  to  0  =  0.  Also  shown  in  Fig.  4 
is  a  beam  "velocity"  line,  w  =  ckg  .  The  points  at  which 
this  line  intercepts  the  dispersion  curves  are  points  at 
which  the  beam  velocity  and  the  phase  velocity  of  the  free 
waves  are  the  same;  they  are  in  "synchronism". 

Consideration  of  eqs.  (11a)  and  (lib)  shows  that, 
at  this  point,  the  angle  of  propagation  is  the  same  as  the 
Cerenkov  angle,  0  =  cos  ^ (l/Bv'e  ) .  At  this  point,  the 

dispersion  is  modified  by  the  beam  and  the  wave  will  grow 
at  a  rate  given  by  eq.  (14).  If  e,  and  the  beam  density 

are  left  unchanged,  the  rate  of  growth  at  the  synchronous 
point  on  the  p.^  curve  will  be  greater  than  that  on  the  p2 
curve  by  an  amount  equal  to  the  frequency  ratio  to  the 
one- third  power.  Thus  the  stimulated  Cerenkov  process  is 
a  potential  short  wavelength  radiation  source. 

Growth  will  also  occur  at  angles  other  than  the 
Cerenkov  angle.  Shown  in  Fig.  5  is  a  numerical  solution 
of  the  complete  dispersion  relation  (eq.  (10)).  We  see 
that  there  are  three  solutions  in  the  positive  ,  positive 
k  quadrant  of  the  w  -  k  plane.  One  is  purely  real,  while 
the  other  two  are  a  complex  conjugate  pair  in  the  region 
below  and  near  synchronism  and  real  above  this  point.  The 
gain  peaks  just  below  synchronism  (the  shift  is  equal  to 
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-1 


the  Re(w-kv0)  given  in  eq.  (13))  and  goes  identically  to 


zero  at  the  point  u>  =  kvQ.  On  the  small  k  side  the 

Im(<i>/c)  goes  to  zero  more  slowly.  The  exact  shape  of  this 

curve  will  depend  upon  y,  e  and  the  beam  density. 

We  have  now  established  that  by  controlling  the 

angle  of  propagation,  e,  and  the  beam  energy,  the  frequency 

1/3 

at  this  maximum  growth  occurs  increases  as  w  .  It  will  be 
instructive  to  consider  the  magnitude  of  the  gain  as  these 
parameters  are  manipulated.  In  order  to  do  this,  we 
rewrite  again  eq.  (14) ,  now  in  this  form: 


/Bw  2w\1/3 

•-w 


g(yt)  F(y,yT) 


(15) 


where 


Y  2  _  e 

T  £-1 


is  the  threshold  energy  (BT  =  1/e), 


U-1AT2>  1/3 


g(yt)  -  - 573 

yt 


and 


F(Y,YT)  = 


d-Yr 


1/3 


(1-1/Y  ) 


(16a) 


(16b) 


(16c) 


so  that 


2 

One  power  of  6  has  been  inserted  in  front  of  w 

P 

we  may  subsequently  express  it  in  terms  of  the  beam  current 

a  form  which  we  will  find  convenient  in  our  numerical 

evaluation  of  the  gain.  Before  we  do  this  evaluation,  we 

will  examine  the  functions  G,  F,  and  yt/y  • 

The  function  G  depends  only  upon  the  dielectric 

constant  of  the  material.  A  sketch  is  shown  in  Fig.  6a. 

It  shows  a  vertical  rise  at  yt  =  1,  the  point  where  the 

dielectric  constant  of  the  material  approaches  infinity, 

2 

reaches  a  maximum  at  y  T  =  7/5  (e  =  7/2) ,  and  finally 
-5/3 

decreases  as  yt  as  YT  becomes  large  (e  1)  .  Thus,  in 
considering  a  practical  Cerenkov  source,  one  cannot  move 
profitably  in  the  direction  of  low  beam  energy,  optically- 
dense  materials  (y  -*•  1/  yt  if  £-►<*>)  since  the  gain 
vanishes  rapidly  in  this  limit.  As  a  practical  matter, 
one  could  not  propagate  a  beam  in  this  type  of  material 
in  any  event.  In  the  opposite  limit,  we  would  have  gasses 
(e-*-l).  In  this  region,  the  gain  will  also  decrease,  but 
conclusions  as  to  the  usefulness  of  this  limit  must  also 
include  consideration  of  the  term.  It  is  interesting 

and  perhaps  important,  for  practical  mm-submm  wavelength 
devices  that  G  peaks  in  the  region  of  the  dielectric 
constant  of  quartz. 

The  function  F  depends  both  upon  the  threshold 
energy,  yt»  and  the  beam  energy  of  y.  It  rises  vertically 
from  y  =  Ym  and  asymptotically  approaches  unity  from  below 
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Sketches  of  F,  yt/y  and  their  product  are  shown  in  Fig.  6b. 
There  is  obviously  a  local  maximum  in  the  growth  rate.  The 
value  of  the  product  at  this  maximum  is  about  .5  . 

Before  we  consider  some  actual  numerical  values  for 
the  growth  rate,  it  will  be  useful  to  consider  one  further 
scaling,  which  will  be  that  of  the  beam  density.  We  assume 
for  the  present  that  the  beam  is  now  a  rectangular  slab  of 
thickness  and  that  the  variation  of  E  in  the  x-direction 
is  still  given  by  exp  (ipx) .  The  term 
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can  be  re-expressed  as 


(17b) 


where 


r 


o 


2.  2 

e  /me 


(18a) 


I  =  ec/r  (18b) 

o  o 

and  I  is  the  electron  beam  current.  When  I  is  measured 
in  amperes,  IQ  has  the  value  17  kA.  Hence,  the  factors 
preceding  the  energy  and  material  form  factors  in  the 
expression  for  gain,  eq.  (15) ,  are  given  by: 


1 


(19) 


When  I  is  approximately  3A,  the  first  of  the  three  factors 
is  approximately  equal  to  0.1,  while  if  it  3ma  it  becomes 
0.01.  The  second  factor  may,  in  principle,  vary  from  zero 
to  a  moderately  large  number,  and  the  characteristic  scale 
length  a  may  be  anything  from  .01  to  1  centimeter. 

Hence,  substantial  gain  is  possible  in  principle.  A 
discussion  of  ways  in  which  this  may  be  achieved  in 
practical  cases  will  be  deferred  until  after  we  have  made 
some  mention  of  wave— guiding  structures. 


II  B.  GAIN  FROM  AN  UNMAGNETIZED  BEAM 


The  preceding  analysis  presumes  that  the  current 
density  modulation  occurs  only  in  the  z-direction.  As  we 
will  see  in  later  discussion,  one  class  of  Cerenkov  device 
will  make  use  of  mildly  relativistic  electron  beams  and 
will  somewhat  resemble  microwave  tubes.  The  beams  in 
these  devices  will  almost  certainly  propagate  along  a 
strong  axial  guide  field,  and,  in  this  limit,  the 
assumptions  made  in  the  last  section  be  at  least  approx¬ 
imately  valid. 

Another  class  of  device,  however,  might  make  use  of 
a  more  relativistic  beam  such  as  that  used  in  the  injector 
of  a  linear  accelerator,  a  linear  accelerator  itself,  or 
perhaps  some  other  type  of  accelerator.  The  beam  in  this 
case  may  very  well  not  be  magnetized.  It  will  then  have 
rapidly  varying  components  in  the  transverse  as  well  as 
the  longitudinal  direction,  and  the  gain  formulas  will  be 
modified. 

When  the  beam  is  unmagnetized,  the  linearized 
equation  for  the  perpendicular  motion  is: 


while  the  longitudinal  motion  is  still  governed  by  eq.  (1) . 
Assuming  the  same  geometry  given  in  Fig.  1,  the  one  non¬ 
vanishing  component  of  this  equation  will  lie  in  the 
x-direction. 
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Equation  (21),  with  the  aid  of  Faraday's  law,  may  be 
restated  in  the  form: 
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The  v  x  B  term  gives  rise  to  an  E  as  well  as  E  dependence 
for  v2. 

The  solution  to  eq.  (21) ,  together  with  the  linearized 
equation  of  continuity,  may  be  used  to  construct  expressions 
for  the  current  density.  These  are: 
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The  current  terms  can  now  be  substituted  in  eq.  (7) .  When 
this  is  done,  we  have  as  our  new  wave  equation: 


*p2vo2^ 

-kvo>7 


=  0  (24) 
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The  determinational  equation  for  eq.  (24)  is  now  the 
dispersion  relation  for  the  unmagnetized  beam-dielectric 
combination.  It  is: 
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(25) 


Equation  (25) ,  which  appears  quite  cumbersome  in  comparison 
with  eq.  (9)  ,  is  still  a  quartic  in  either  <d  ,  k  or  both. 

All  qualitative  comments  made  about  the  strongly  magnetized 
case  apply  here  as  well.  However,  the  results  are 
quantitatively  somewhat  different.  Again,  the  strongest 
coupling  region  of  the  beam  to  the  wave  is  in  the  velocity 
synchronism  (u>/ck  =  B). 

2 

If  terms  proportional  to  l/(id-kvQ)  are  collected 
separately,  we  obtain  for  the  dispersion  relation: 


Near  synchronism,  this  reduces  to: 


td2(c-l)  (B2e-1) 
((d-kvQ)  2 


(27a) 
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Once  again,  the  dispersion  relation  is  cubic  and  the 

frequency  and  the  dependence  upon  the  size  of  the  Cerenkov 
-1  2 

angle  (0  =  sin  (1-1/3  e ))  are  the  same.  However,  the  beam 
energy  and  e  dependence  are  different.  If  we  use  the 
functions  defined  earlier,  we  have  for  the  spatial  growth 
rate: 


a 


G(Y  )  f(y,yt) 


1/3 


(28) 


-1/3  •  • 

The  energy  dependence  is  now  y  in  the  high  energy  limit, 

as  opposed  to  the  more  constrictive  y  1  dependence  in  the 
strongly  magnetized  limit.  If  all  other  factors  are  the 
same,  the  gain  in  the  unmagnetized  limit  will  be  greater 
than  that  for  the  strongly  magnetized  beam.  This  is 
because  the  electrons  in  the  beam  can  now  do  work  on  the 
wave  in  both  the  transverse  and  longitudinal  direction. 
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II  C.  BOUNDED  STRUCTURES 


Excepting  the  possibly  interesting  limit  of  extremely 
relativistic  beams  and  gaseous  dielectrics,  it  is  not 
practical  to  have  the  beam  penetrating  the  dielectric. 
Hence,  in  assessing  the  practicality  of  Cerenkov  sources, 
it  is  important  to  consider  dielectric  wave  guides  and 
resonators  which  have  channels  for  the  beam  propagation. 
This  complicates  the  analysis.  Thus,  before  we  take  up 
the  cases  quantitatively,  it  will  be  useful  to  consider, 
at  this  point,  the  regime  where  the  results  of  the 
preceding  section  are  qualitatively  useful. 

First,  we  note  that  with  minor  changes,  the  results 
of  the  last  section  will  apply  exactly  to  a  metal-bounded, 
cylindrical,  dielectric  waveguide  through  which  an 
electron  beam  propagates.  The  perpendicular  wave  number, 
p,  is  now  a  root  of  zero  order  Bessel  function  and  is  no 
longer  completely  free.  The  only  other  change  is  that  the 
factor  7r  in  the  current  term  no  longer  appears,  because 
the  beam  is  now  also  cylindrical.  The  field  symmetry  is 
now  transversely  magnetic. 
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II  C.  1.  Cylindrical  Guide  with  a  Beam  Channel 

When  the  beam  propagates  in  a  hole  in  the  dielectric, 
we  have  a  situation  such  as  that  sketched  in  Fig.  7. 

If  the  diameter  of  the  hold  is  sufficiently  small,  a 
concept  which,  shortly,  will  be  quantitatively  specified, 
the  results  of  the  preceding  section  might  be  expected  to 
apply  more  or  less  exactly. 

It  is  obviously  the  relative  size  of  the  hole  which 
is  the  fundamental  difference.  Fortunately,  it  is  possible 
to  attain  considerable  insight  into  its  effect  with  little 
analysis.  We  consider,  for  the  moment,  a  metal-lined 
guide  partially  filled  with  dielectric.  The  dispersion 
curves  sketched  in  Fig.  7  are  similar  to  those  shown  in 
Fig.  4.  The  main  difference  is  the  shape  near  the  light 
line,  ui  =  ck.  The  point  where  the  curve  crosses  this  line 
is  now  controlled  by  the  relative  filling  factor,  d/b,  as 
well  as  the  dielectric  constant  of  the  material.  As  d/b 
and  e  become  small,  the  point  where  the  partially-filled 
guide  becomes  a  slow  wave  structure,  w/ck  <  1,  can  thus 
still  be  made  to  occur  at  an  arbitrarily  high  frequency. 

When  u/ck  >  1,  outside  the  light  line,  the  field  in 
the  hole  is  proportional  to  JQ(pr),  an  ordinary  Bessel 
function.  In  this  regime  it  peaks  in  the  center  of  the 
hole.  However,  we  must  operate  in  the  regime  u/ck  <  1, 
and  in  this  case,  the  radial  dependence  is  proportional 
to  a  modified  Bessel  function  ,  IQ(qr)  .  The  field  is  now 
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a  minimum  at  r  =  0,  and  the  beam  wave  coupling  is 
obviously  decreased. 

A  sketch  of  the  field  dependence  in  the  two  regimes 
is  shown  in  Fig.  8.  The  wave  number  in  the  dielectric,  p, 
is  still  given  by: 


(29) 


while  the  wave  number  in  the  hole,  when  w/ck  <  1,  is  now 
given  by 


q 


2 


(30) 


The  latter  is  obviously  one  measure  of  the  field 
depression  in  the  hole.  Since  we  operate  near  synchronism, 
a)  =  ck8  ,  we  have  for  q: 

q  =  k/y  Ola) 

q  =  u)/c6y  (31b) 

or 

q  =  2tt/16y  (31c) 

Hence,  when  non-relativistic  beams  are  used  fiyiav/c,  the 
field  drops  off  away  from  the  dielectric  in  a  distance 
small  compared  to  a  wavelength.  If,  however,  the  beam  is 
at  least  mildly  relativistic.  By  £  1,  the  opposite  limit 


applies  and  we  can  operate  with  wavelengths  that  are  small 
compared  to  the  hole. 

The  latter  considerations  actually  apply  to  any 
structure  supporting  a  wave  for  which  u>/ck  <  1.  One 

might  then  ask  about  the  relative  advantages  of  a 
dielectric  tube  since,  if  8y  1,  then  coupling  would  be 
improved  at  short  wavelengths  for  only  slow  wave  structure. 

The  advantages  of  the  dielectric  tube  also  lie  in  the 
short  wavelength  range.  In  a  conventional  slow  wave 
structure,  the  periodicity  must  also  be  comparable  to  the 
wavelength.  Structures  of  reasonable  length  are,  therefore, 
a  great  many  wavelengths  long  and  they  become  very  difficult 
to  fabricate  at  relatively  long  wavelengths  (a  few  milli¬ 
meters)  .  It  is  possible,  but  not  easy,  to  build  conventional 
structures  with  a  fundamental  period  smaller  than  a  few 
millimeters.  The  dielectric  is,  however,  a  smooth  structure 
and  easy  to  fabricate.  When  the  beam  is  relativistic,  the 
coupling  impedance  becomes  comparable  to  that  of  other 
structures.  Modifications  of  this  basic  structure,  such  as 
a  dielectric  tube  with  no  metal  boundary  and  multiple 
coupled  tubes,  may  also  be  of  practical  use.  Another  basic 
structure,  a  dielectric  slab  bounded  on  one  side  by  a 
conductor,  also  shows  promise  for  application  in  the 
shorter  wavelength  region.  This  follows  from  the  fact  that 
a  greater  mode  separation  at  small  wavelengths  can  be 
obtained  from  this  more  open  structure.  Hence,  it  may  well 


be  easier  to  make  single-mode  devices  with  this  type  of 
structure,  and  for  this  reason  we  will  analyze  it  in  some 
detail . 

The  basic  geometry  is  shown  in  Fig.  9.  Assuming, 
for  the  moment,  that  no  beam  is  present,  we  have  for  the 
TM  modes  of  the  guide: 


E  =  <0,Ey,Ez) 


(32a) 


where 


[$  •&->■)] 


k  )  I  E_  =  0 


and 


ik  «, 
«.2e-k2  3x 


(32b) 


(32c) 


In  the  region  0  —  y  sS  d  ,  the  dielectric  constant  e  appears 
while  in  the  region  y  *£  d,*»«"  is  set  equal  to  unity. 

Anticipating  the  fact  that  we  are  concerned  only  with 
slow  waves  bound  to  the  surface  guide,  we  have  for  the 


electric  fields: 


s 


1 1 

!  r 


b 
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where 


(33b) 


in  the  region  0  £  y  ^  d  .  Outside  the  dielectric,  the  field 
is: 


Ez  =  Be“qy  (34a) 

2  2  2  2 

q  =  k  u>  /c  (34b) 

Matching  the  tangential  electric  and  magnetic  fields  can  be 
used  to  eliminate  the  constants  A  and  B.  Thus  we  have  for 
the  dispersion  relation  of  the  dielectric  slab  wave  guide: 

eq  cot  pd  =  p  (35) 

A  plot  of  the  roots  of  this  function  is  given  in  Fig.  10. 
The  lowest  order  mode  has  no  cutoff.  It  comes  up  along  the 
light  line,  u>  =  ck,  until  pd  gets  somewhat  closer  to  the 
neighborhood  of  tt/2.  Thereafter,  as  w  becomes  larger,  it 
asymptotically  approaches  the  speed  of  light  in  the 
dielectric.  In  the  region  tt/2  <s  pd  3s  tt  there  are  no 
solutions  to  q.  (45),  while,  when  tt  i  pd  <  3ir/2,  a  second 
mode  which  has  a  finite  w  cutoff  frequency  can  also  propa¬ 
gate.  At  successively  higher  frequencies,  more  of  these 
inodes  appear.  Several  are  shown  in  Fig.  10. 


II  C.  2.  Coupling  of  Beam  to  Bounded  Resonator 


Also  shown  in  Fig.  10  is  a  beam  speed  line,  w  =  ckg  . 
It  is  obvious  that  if  the  beam  velocity  satisfies  the 
Cerenkov  conditions, 6  >  1/Ve  ,  phase  synchronization 
between  an  electron  beam  and  a  wave  can  be  maintained. 

When  the  beam  is  added,  the  wave  equation  in  the  vacuum 
region  becomes: 


e,i  ~  1  “ 


2,  3 
MP  /Y  ^ 

(u>~kvo) 4 


(36a) 


(36b) 


In  arriving  at  eqs.  (36a)  and  (36b) ,  it  has  been  presumed 
that  beam  density  modulation  occurs  only  in  the  z-direction, 
that  the  left  edge  of  the  beam  is  close  to  the  dielectric, 
and  that  the  beam  extends  indefinitely  in  the  region  y  >  d. 
The  size  of  the  actual  gap  between  the  beam  and  the 
dielectric  will  be  an  important  parameter  in  a  short  wave¬ 
length  device  and  its  role  will  be  discussed  separately. 


II  C.  3.  The  Beam-Guide  Dispersion  for  a  Bounded  Structure 

When  a  bounded  structure  is  used  to  support  the  wave, 
as  it  must  in  almost  any  practical  source,  the  dispersion 
relation  becomes  a  transcendental  as  opposed  to  an  algebraic 
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function.  It  will  be  more  or  less  straightforward  to 
obtain  values  for  the  roots  by  numerical  means,  but  it  is 
not  immediately  obvious  how  to  obtain  a  good  qualitative 
understanding  of  the  roots. 

One  method  which  is  appropriate  for  relatively  weak 
beams  is  the  following.  Assume  a  relation  of  the  form: 


D(u),  k,  o>  )  =  0 


where  the  presence  of  in  eq.  (37)  indicates  the  presence 
of  the  beam.  If  the  beam  is  weak,  we  can  write: 


D  (a),  k,  a)  2)  =  D (0)  (a),  k)  +  a)  2 

**  ^  3<d 

P 


(38a) 


where 


(w,  k)  =  D  (k,  a),  0) 


(38b) 


is  the  dispersion  relation  for  the  waves  supported  by  the 
structure  when  no  beam  is  present.  This  function  can,  in 
a  region  near  to  the  solution  D ^  (id,  k)  =  0  be  written 


D(0)  <(d,  k)  =  (id— u>k)  3D°/8w  (39) 


where  are  the  roots  of  eq.  (35) . 

The  second  term  in  eq.  (38a)  can  also  be  further 
reduced.  The  dependence  of  the  dipsersion  relation  upon 
(Dp  always  enters  through  e  and  hence  the  second  term  of 


Equation  (38a)  will  also  have  the  form: 


2 

0) 

P 


3D 


3  co 

P 


to  2  C  (to  ,k) 

3  .  ,  .  3 

Y  (to-kvo) 


(40) 


where  C(io,  k)  is  a  function  which  will  depend  upon  the 
detail,  of  the  structure.  It  may,  for  example,  have  zeros, 
but  it  will  not  have  any  poles  near  either  to  =  ky,  or 
to  =  to^. 

Thus,  near  synchronism,  to.  =  kv  ,  and  for  beams 

JC  o 

which  are  not  too  strong,  I/Iq  <<  1»  we  again  have  a  cubic 
dispersion  relation: 


or 


0  =  (to— cok) 


+  3.“p___  ,i>  c(“k’k' 


Y  (<o-kv0) 


(41a) 


(co-kvQ)  3 


to  2  C(to.  ,k) 

P  - - - 

Y3  3D°(toK,k)/3to 


(41b) 


Thus,  the  qualitative  nature  of  the  roots  is  rather 
independent  of  the  exact  geometry  of  the  wave-supporting 
structure. 

When  the  wave-supporting  structure  is  a  dielectric 
slab  and  the  assumptions  made  earlier  apply,  the  dispersion 


relation  becomes: 


The  expansion  procedure  outlined  in  the  preceding  sub¬ 
section  then  gives  for  eq.  (41b) : 


(•-?)  *  -(¥-) 


(1-1/Be) 

04  3 
Y 


.  2  , 
sin  pd 

i 


(43) 


kd  Y"  .2 

—  + - j  sin  Pd 

'  eYt 


The  first  two  groups  of  factors  on  the  right-hand  side  of  eq. 
(43)  are  identical  to  the  results  obtained  when  it  was  assumed 
that  the  beam  propagated  in  the  dielectric,  and  much  of  the 
discussion  presented  at  that  point  applies  here  as  well.  The 
last  group  of  factors  contains  the  dependence  on  the  geometry. 
It  can  be  seen  that,  in  addition  to  the  Cerenkov  threshold 
dependence,  the  coupling  also  goes  to  zero  as  the  thick¬ 
ness  of  the  slab  goes  to  zero,  and  is  a  result  that  could 
easily  be  anticipated. 

The  other  trends  in  the  gain  can  be  understood  as 

follows:  On  the  fundamental  mode,  the  value  of  pd  varies 

from  0  at  (i)  =  0  up  to  7r/2  as  to,  k  °°.  On  the  higher 

branches,  it  varies  from  nrc  at  cutoff  (to  =  ck)  to 

(2n  +  1)it/2  as  the  curve  asymptotically  approaches  the 

2 

speed  of  light  in  the  dielectric.  The  value  of  sin  pd 

thus  varies  monotonically  from  zero  to  one.  Assuming  that 

the  velocity  synchronism  is  maintained  along  the  dispersion 

curve,  the  gain  will  vanish  at  to  =  ck,  because  in  this 

limit,  y  and  it  becomes  increasingly  difficult  to 

modulate  the  beam.  Furthermore,  as  0  •+  1  VT  ,  the  gain 

2 

also  vanishes  due  to  the  factor  (1-1/0  e)  in  the  coefficient. 
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The  gain  thus  vanishes  at  both  ends  of  the  dispersion  curve 
and  peaks  in  between.  A  sketch  of  the  general  behavior  is 
shown  in  Fig.  11. 

The  maximum  value  that  the  gain  can  achieve  is  similar 

to  that  of  the  filled  guide  case.  Some  typical  results  for 

a  thin,  quartz  slab  waveguide  are  shown  in  Figs.  12a  and 

12b.  In  these  plots,  the  factor  (3u)  2/c2  1//3)  has  been 

P 

omitted  for  convenience.  The  remaining  factors  contain  all 

the  relevant  frequency  and  energy  dependencies.  Maximum 

values  somewhat  greater  than  unity  are  obtain  for  this 

2  2  1/3 

particular  set  of  parameters.  The  omitted  term  (3u>  /c  ) 

r 

2 

is  actually  the  beam  current  density  in  A/cm  divided  by 
I  (  =  17  kA)  all  to  the  one-third  power.  It  is  relatively 
easy  to  obtain  values  of  0.1  for  this  number/  hence  the 
plots  shown  in  Figs.  12a  and  12b  demonstrate  that  with  a 
quartz  slab  guide,  it  is  possible  in  principle  to  have 
relatively  large  gain  (a  =  .233)  gives  (ldb/cm)  well  into 
the  submillimeter  part  of  the  spectrum. 

The  gain  plot  in  Fig.  12  also  indicates  that  the  gain 
is  a  bit  higher  on  the  higher  order  modes.  This  trend  is 
a  reflection  of  the  factor  in  the  gain.  It  is  real,  but 

it  depends  upon  two  assumptions  whose  validity  are  also 
frequency  dependent.  These  are:  first,  that  the  beam  is 
infinitesimally  close  to  the  dielectric  and  second,  that  the 
beam  is  monoenergetic .  The  first  of  these  will  be  discussed 
now  and  the  second  point  will  be  covered  in  a  later  section. 


-35- 


i 

r 

K3 

l 


II  C.  4.  Finite  Gap  Between  Beam  and  Resonator 

If  we  assume  that  there  is  a  small  gap  between  the 
beam  and  the  dielectric  surface,  we  would  have  a  situation 
such  as  that  shown  in  Fig.  13.  The  analysis  procedes  as 
before,  but  the  resulting  dispersion  relation 

(44) 

is,  at  first  sight,  much  more  complicated.  However,  if  we 
again  assume  that  the  roots  at  synchronism  lie  along  the 
dispersion  curve  for  the  free  modes,  the  situation  simpli¬ 
fies  considerably  and  the  end  result  is  that  the  gain  is 
modified  by  an  exponential  factor  which  depends  upon  the 
size  of  the  gap: 

a  =  ot(d)  e  2 /Y  (45) 

As  long  as  kd 2 /y  is  small,  the  gain  on  the  higher  order 
modes  will  be  comparable  to  or  greater  than  the  gain  on 
the  fundamental  mode.  Values  of  62  of  about  1  millimeter 
would  be  conservative  and  fractions  of  this  are  easily 
obtained.  Hence,  provided  that  ones  uses  0y  st,  1,  the 
quart  guide  system  discussed  above  will  still  be  viable 
well  into  the  submillimeter  region  of  the  spectrum. 

We  have  been  assuming  that  the  beam  extends  indefini¬ 
tely  in  the  positive  y-direction.  As  long  as  the  beam  is 


at  least  a  few  e-foldings  thick,  this  assumption  does  not 
affect  the  gain.  Since  we  are  primarily  interested  in  high 
frequencies,  this  assumption  will  normally  be  valid. 

The  fall-off  of  the  field  in  the  transverse  direction 
may  also  be  useful  in  obtaining  some  mode  selection.  If  a 
relatively  thin  beam  is  used,  the  fields  for  the  lower  order 
modes  may  penetrate  through  to  the  other  side.  If  a  lossy 
material  is  placed  above  the  beam,  it  may  be  possible  to 
further  reduce  the  gain  on  the  lower  order  mode. 

In  the  fall-off  of  the  electric  field,  operation  at 
arbitrarily  short  wavelengths  could  be  obtained  if  y  is 
allowed  to  become  large,  i.e.  kd2/y  will  remain  small. 

This  will  involve  a  penalty  in  the  maximum  value  of  gain 
obtainable,  but  since  it  is  relatively  large  to  begin  with, 
the  resulting  system  will  still  be  potentially  useful.  In 
this  way,  with  more  relativistic  e- beams,  it  might  be 
possible  to  operate  well  into  the  infrared  portion  of  the 
spectrum.  This  will  be  discussed  further  in  a  later  section. 


II  D.  THE  EFFECT  OF  BEAM  VELOCITY  SPREAD 


Prior  to  this  point  in  our  discussion,  we  have  assumed 
that  the  electron  beam  was  perfectly  monoenergetic.  It  is 
intuitively  plausible  that  this  is  a  wavelength-dependent 
assumption,  and  we  will  now  examine  its  consequences.  The 
discussion  will  be  divided  into  three  parts.  First,  we  will 
determine  wavelength  limit  for  a  simple  beam  space  charge 
wave.  Then,  this  result  will  be  compared  with  a  similar 
criterion  for  a  Cerenkov  instability.  Finally,  having  set 
the  limiting  wavelength  for  treating  the  beam  as  mono- 
energetic  we  will  derive  gain  expressions  valid  in  the 
region  where  the  assumption  is  violated. 

II  D.  1.  Beam  Space  Charge  Waves 

The  linearized  equation  of  motion  for  a  strongly- 
magnetized  electron  beam  is  given  by  eq.  (1) .  If  this  is 
taken  along  with  the  equation. of  continuity  eq. -(3a),  Poisson's 
equation,  and  assumptions  similar  to  those  of  that  section,  the 
dispersion  relation  for  space  charge  waves 

to  =  kvQ  ±  Wp/y3^2  (46) 

may  be  easily  derived. 

The  upper  (lower)  sign  in  eq.  (46)  corresponds  to  a 
fast  (slow)  space  charge  wave.  See  Fig.  14a. 
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We  will  concern  ourselves  with  a  slow  space  charge 
wave.  Shown  in  Fig.  14b  is  a  sketch  which  illustrates  the 
meaning  of  the  statement:  "the  wave  is  resolved  from  the 
beam" .  The  wave  is  clearly  resolved  when  the  beam  may  be 
regarded  as  a  delta  function  in  frequency  space,  (the  arrows 
located  at  u>  and  kvQ)  .  If ,  on  the  other  hand,  the  velocity 
spread  of  the  beam,  Av,  is  such  that  the  self-consistent 
frequency  separation,  Acj  =  oi  -  kvQ  (derived  under  the 
assumption  that  the  beam  was  monoenergetic)  is  less  that 
kAv,  the  assumption  is  violated.  A  quantitative  criterion 
for  this  critical  k  is: 


kcAV 


(47) 


I 


Equation  (47)  may  be  re-expressed  in  terms  of  physically 
more  intuitive  variables  if  we  write:  kc  =  u>/cB  =  2tt/AcB, 
Av  in  terms  of  Ay  and  0 ,  and  Up  in  terms  of  the  beam 
current  density  J^.  Then  we  have: 


X 


c 


TT 


(48) 


where  I  is  still  ec/r  =  17  kA.  If  By  is  of  order  unity, 
o  o 

Ay/y  is  of  order  10~2  and  Jb  is  a  reasonable  fraction  of  an 
ampere/millimeter2  than  Ac  is  a  fraction  of  a  millimeter. 
These  are  relatively  modest  requirements,  and  thus  we 
predict  that  it  should  be  possible  to  make  effectively 
cold  beams  well  into  the  submillimeter  part  of  the  spectrum. 


i 
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The  critical  wavelengths  given  by  eqs.  (47)  and  (48) 

are  dependent  upon  the  assumption  of  a  simple  space  charge 

wave.  When  we  are  considering  a  Cerenkov  instability, 

3/2 

however,  Au®  w  -  kv  is  actually  larger  than  w  / y  '  ,  and 

\J  p 

hence  the  beam  can  be  effectively  colder  at  a  given  wave¬ 
length.  The  criterion  for  resolution  is: 


WI 

kcAv  =  T 


(49) 


u 


where  the  right-hand  side  of  eq.  (49)  is  the  real  part  of 
the  detuning  (eq.  13b) .  Substitution  of  the  expressions 
for  cl> j  can  be  made  for  the  appropriate  case. 

When  the  beam  propagates  through  the  dielectric, 
eq.  (13b)  applies  directly  and  we  have: 


V ^  (frf2  (i^) 


^  1 


’77 


777- 


(50) 


The  current  density  dependence  is  similar  to  that  of  eq.  (48) , 

but  provided  the  beam  is  at  least  mildly  relativistic 

(By  £  1),  the  energy  dependence  is  more  favorable.  Overall, 

presuming  that  yT  and  y/yT  are  not  excessively  large,  the 

value  A  ,  given  by  eq.  (50) ,  will  be  at  least  as  small  as 
c 

that  given  by  eq.  (48) .  The  addition  of  the  form  factor 
associated  with  a  more  practical  resonator  will  not  alter 
this  essential  conclusion. 
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II  D.  2.  Gain  in  the  Warm  Beam  Limit 

When  the  criteria  given  in  the  preceding  paragraphs  are 
violated,  the  beam  is  to  be  regarded  as  "warm"  at  the  wave¬ 
length  in  question.  The  gain  does  not  vanish  in  this  limit, 

but  it  does  begin  to  drop  as  w  ^  ,  as  opposed  to  the  general 
1/3 

a)  '  trend  in  the  cold  beam  limit.  This  trend  means  that 
oscillators  can  probably  be  built  in  the  warm  beam  limit, 
but  amplifiers  will  be  impractical. 

In  calculating  the  gain,  we  will  use  the  Vlasov  equation 
as  the  basic  equation  of  motion  and  we  will  retain  the 
assumption  that  the  beam  is  strongly  magnetized.  In  this 
case,  the  Vlasov  equation  is: 


3f 

at 


+ 


3f  .  r, 
li  +  Ps 


(51) 


If  this  is  linearized,  f  =  f  +  if,  and  is  Fourier-transformed, 
we  have  for  the  perturbed  component  distribution: 


6  f 


ie  E. 


3f  /  3p 
o  z 
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(52) 


The  current  is  now  given  by: 


*  f 


v  5  f  3p 
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(53a) 


_  f  <*«/*.>  a 

le  Ez  vz  M-kv_  dvz 


(53b) 


Substitution  of  this  into  eq.  (8)  will  then  lead  to  a 
dispersion  relation.  If  the  beam  distribution  is  a  delta 
function,  then  the  integral  can  be  performed  immediately 
and  the  current  given  by  eq.  (4b)  is  recovered.  However, 
we  are  now  interested  in  the  limit  where  the  beam  velocity 
spread  is  finite. 

An  exact  solution  of  a  dispersion  relation  containing 
a  integral  kernel,  such  as  that  of  eq.  (54) ,  can  be  found 
using  numerical  techniques.  The  results  of  such  a  procedure 
will  be  discussed  below.  However,  some  insight  into  the 
general  behavior  can  be  obtained  in  the  limit  where  k 
times  the  width  of  the  beam  distribution  is  broad  in 
comparison  with  the  gain  which  would  be  obtained  from  a 
calculation  in  which  it  were  assumed  that  the  beam  were 
cold  (monoenergetic) . 

The  dispersion  relation  obtained  from  the  above 
procedure  is: 


In  its  present  form,  the  integral  which  appears  in  eq.  (54) 
is  to  be  performed  along  the  real  momentum  line,  and  hence 
the  procedure  for  handling  the  singularity  at  synchronism 
is  not  yet  defined.  Borrowing  from  the  theory  of  plasma 
physics,  we  handle  it  by  formally  extending  the  integral 
into  the  complex  plane.  First,  we  re-express  eq.  (54)  as  a 
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velocity  integral.  Then: 


(g>2e/c2-k2) 


id ‘e  ,2  2  .  2  u  e/c 
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we  now  let 


1 

w-kv 


where  P  stands  for  the  principal  part  of  the  integral. 

2 

Multiplying  through  by  c  /e  and  using  the  condition  v  =  oo/kioe, 
we  find  for  the  imaginary  part  of  the  dispersion  relation: 


D"  =  - 


tto)  2B2(1-1/B2c) 


8F  (co/k) 
3(w/k) 


It  is  sufficient  for  the  purpose  of  the  present  discussion 
to  ignore  the  small  correction  to  the  real  part  of  the 
dispersion  represented  by  the  principle  part  of  the  integral. 
In  the  limit,  the  real  part  of  the  dispersion  is: 


D 1  =  iti  -  co. 


(58a) 


and  providing  the  growth  is  small,  the  imaginary  part  of  the 
frequency  is  adequately  represented  by: 


D" 

OD'/aw 


(58b) 
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(58c) 


Obviously,  there  will  be  wave  growth  (inverse  Landau  damping) 

in  the  region  of  velocity  space  where  a>k/kcvo)  .  Sketches 

of  the  unperturbed  dispersion  relation,  oj  =  oj^  and  are 

shown  in  Fig.  15.  The  region  of  positive  w  lies  on  the 

larger  k  side  of  the  synchronous  wave  number,  k  =  w/v  , 

s  o 

and  peaks  at  a  velocity  which  is  below  vq  by  an  amount 
approximately  equal  to  the  half  width  of  the  velocity 
distribution. 

Thus,  the  wave  growth  in  the  warm  beam  limit  has  a 
shape  which  is  complimentary  to  that  obtained  in  the  cold 
beam  limit.  This  result  may  seem  to  imply  that  the  growth 
due  to  inverse  Landau  damping  is  fundamentally  different 
from  the  growth  obtained  in  the  warm  beam  limit.  However, 
this  is  not  the  case.  If  the  roots  of  eq.  (55)  are  followed 
as  the  beam  width  is  varied  from  a  value  of  zero  up  to 
kAVujj  (cold)  ,  we  find  that  one  regime  passes  smoothly  into 
the  other.  The  peak  absorption  shown  in  Fig.  15  is  always 
somewhat  less  than  the  gain.  It  is  a  composite  of  the  long 
wavelength,  cold  beam  gain  and  Landau  damping  caused  by  that 
part  of  the  beam  which  has  3F/3v  <  0.  The  peak  gain 
obtainable  in  the  warm  beam  limit  will  always  be  less  that 
Uj  (cold) .  This  occurs  because,  as  the  beam  distribution 
becomes  arbitrarily  sharp,  the  self-consistent  beam  density 


dependent  frequency  shift  will  move  the  phase  velocity  of 
the  wave  downward  relative  to  the  position  of  maximum 
3F/9V. 

In  spite  of  this,  the  warm  beam  limit  may  very  well  be 

of  interest.  An  assessment  of  this  requires  that  we  estimate 

eq.  (58c)  which,  in  general,  depends  upon  the  detailed  shape 

of  the  beam  distribution.  We  can,  however,  proceed  without 

undue  complication  if  we  recognize  that  if  ***  1/Av,  and 

2 

thus  3F/yv  at  its  maximum  is  v*  1/Av  .  Then,  in  terms  of 
Ay,  an  approximate  expression  for  becomes: 


U)  « 
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(59a) 


or 


(59b) 


It  is  interesting  to  evaluate  the  possible  gain  in  the  10y 

-3 

range.  In  this  case,  A  =  10  centimeters.  If  we  can 

-3  -2 

achieve  J^/IQ  ^  10  and  A3/y  >10  ,  it  appears  that 

a  v*  .1  cm.  1  are  within  the  realm  of  possibility.  Thus,  it 
might  be  possible  to  construct  Cerenkov  lasers  down  to  wave¬ 
lengths  comparable  to  those  achieved  in  stimulated  Compton 


devices. 


It  is  also  interesting  to  evaluate  eq.  (59b) ,  when  A 


is  equal  to  A .  Upon  substitution,  we  obtain: 

■  =^0)'  (60 

Examination  of  eq.  (60)  supports  the  conjecture  that 
a  ^.  =  .1cm.  1  are  attainable  in  cold  to  warm  crossover 
region  (A  sc  A  )  . 
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II  E.  COMMENTS  ON  NONLINEAR  BEHAVIOR 


It  would  be  possible  at  this  point  to  develop  a  reasonably 
complete  nonlinear  theory  for  single-mode  Cerenkov  devices. 

This  would,  in  part,  follow  lines  of  argument  originally 
established  to  explain  microwave  tubes,  beam  plasma  inter¬ 
actions,  and  more  recently,  f ree-electron  lasers.  in  the 
wavelength  regime  which  is  of  primary  interest,  however,  the 
ratio  L/X  is  large  and  there  will  be,  at  the  very  least,  a 
few  axial  modes  within  the  half -width  of  the  gain  curve. 

There  may  also  be  a  mixture  of  transverse  modes,  although  if 
the  device  can  be  made  to  operate  on  a  single  transverse 
mode,  it  will  be  advantageous  to  do  so.  Relatively  less  is 
known  about  electron  beam  devices  in  the  multimode  region, 

and  the  development  of  a  complete  nonlinear  theory,  analogous 
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to  that  developed  by  Lamb  for  the  gas  laser,  should  take 
account  of  multimode  operations.  This  would  be  a  substantial 
undertaking.  It  could  be  productive,  however,  since  we 
generally  expect  that  Cerenkov  devices  will  exhibit  many 
phenomena  intrinsic  to  all  multimode  oscillators,  and  that 
some  of  these  may  be  useful  in  applications  (e.g.,  mode 
locking) .  Hence,  because  the  development  would  at  this 
point  omit  some  of  the  most  interesting  parts  of  the  problem 
(on  account  of  its  length  and  because  the  motivation  for 
experimental  development  rests  primarily  on  the  prediction 
of  the  linear  gain  which  would  be  expected  in  specific  wave¬ 
length  regions)  we  will  restrict  discussion  of  nonlinear 
problems  to  a  few  simple  scaling  arguments. 


II  E.  1.  Nonlinear  Scaling  Arguments 


In  the  operating  regime  where  the  beam  velocity 
distribution  can  be  regarded  as  cold  and  the  motion  one¬ 
dimensional,  the  relative  density  modulation  is  given  by: 


6n 

n 

o 


kv 

w-kv 

o 


(61) 


Now  the  change  from  electron  orbits,  which  move  progressively 
forward  compared  to  the  phase  of  the  wave  (untrapped)  to 
orbits  which  are  winding  up  (trapped) ,  occurs  in  the 
vicinity  of  |6n/nQ|  1.  Thus  it  occurs  where 


kfiv  iss  |  in  -  kvQ| 
or 


(62a) 
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(62b) 


This  can  be  converted  to  a  prediction  of  the  magnitude  of 
the  axial  component  of  the  electric  field  at  which  saturation 
of  the  linear  growth  is  in  progress.  We  find  from  eq.  (1) 
and  eqs.  (61a)  and  (61b) : 


E 
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(63) 


The  Poynting's  flux, and  hence  also  the  total  power  carried 

o 

to  the  wave,  will  be  proportional  to  E  and  thus  up  to 

z 

The  latter  will  in  turn  be  proportional  to  (I/IQ)4/^. 
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Hence,  up  to  form  factors  (which  can  vary  between  small 

numbers  and  unity)  the  overall  power  at  the  separatrix 

crossing  value  of  E  (eq.  (62))  will  be: 

z 


p 


wave 


beam 


(64) 


This  is  a  conversion  factor  which  can  be  expected  in  any 
traveling  wave  device.  Enhanced  conversion  could  be  obtained 
by  tapering  the  phase  velocity  and  thus  "deepening"  the 
trapping  well.  Still  more  energy  could  be  recovered  if 
the  beam  were  collected  at  high  potential  (depressed 
collector  operation) .  An  estimate  of  the  overall  efficiency 
obtainable  from  a  ferenkov  device  is  thus  a  subtle  and 
complex  matter.  In  the  final  analysis,  however,  tube-like 
efficiencies  of  perhaps  fifty  percent  could  be  attainable. 


III.  A.  THE  ELECTRON  BEAM 


The  single  most  important  component  of  an  electron 
beam-driven  radiation  source  is  the  beam  itself.  Its 
parameters  will,  in  large  part,  determine  the  performance 
of  the  system.  In  order  to  examine  the  potential  of 
Cerenkov  sources  in  the  millimeter  and  submillimeter 
parts  of  the  spectrum,  it  will  be  useful  to  examine  the 
parameters  of  some  typical  electron  beam  generators  which 
may  be  used  in  this  application. 

A  few  general  types  of  electron  beam  generators  and 
their  parameters  are  listed  in  Table  1.  All  of  the  beams 
are  at  least  mildly  relativistic;  a  choice  is  dictated  by 
coupling  considerations  developed  in  earlier  sections. 
However,  the  values  of  the  beam  current  and  the  modes  of 
operation  vary  widely.  We  will  begin  our  discussion  of 
the  entries  in  the  table  by  considering  the  role  of  the 
beam  energy. 

The  beam  energy  helps  to  determine  the  operating  wave¬ 
length  in  several  interrelated  ways:  first,  by  synchronism; 
second,  by  the  magnitude  of  the  gap  between  the  beam  and 
the  resonator,  which  must  be  present  in  any  real  system;  and 
third,  by  its  entry  into  the  equations  which  determine  the 
beam  modulation.  The  first  of  these  alone  does  not  place 
any  stringent  limit  on  the  attainable  wavelength.  This  is 
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because  the  design  of  dielectric  resonators  which 
will  support  a  wave  of  any  reasonable  phase  velocity  does 
not  present  a  problem.  The  second  and  third  aspects  of  the 
energy  dependence  will  therefore  be  more  important  in 
setting  the  short  wavelength  limit  to  a  device.  The  rate 
at  which  the  electric  field  decreases  as  the  distance  from 
the  resonator  increases  is  given  by  q  =  2tt/13y.  The  gain, 
however,  increases  with  frequency  and  hence  we  would  expect 
it  to  peak  somewhere  in  the  vicinity  of  qaistl,  where  a 
is  a  characteristic  distance  between  the  beam  and  resonator. 

A  modest  value  for  a  would  be  approximately  one  milli¬ 
meter,  and  a  more  difficult  but  attainable  value  would  be 
about  one-tenth  of  this.  Taking  this  range  and  assuming 
qa  m  1,  we  can  determine  the  limiting  wavelength  for  good 
beam-to-resonator  coupling.  The  range  of  this  wavelength 
is  shown  in  the  table.  Examination  of  the  table  and  the 
figure  will  show  that  relatively  compact  machines  could 
ultimately  be  expected  to  work  well  into  the  submillimeter 
region  of  the  spectrum,  while  if  one  extends  the  range  of 
beam  generator  complexity,  operation  in  the  infrared  could 
be  possible. 

While  the  beam-to-resonator  coupling  decreases  away  from 
the  resonator  due  to  the  fall-off  of  the  electric  field  with 
distance,  an  arbitrary  increase  in  beam  energy  will  lower 
the  minimum  wavelength  by  a  corresponding  amount.  However,  the 
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beam  energy  also  enters  into  consideration  through  the 
equations  of  motion  for  the  electrons,  and  increasing  y 
furthers  the  difficulty  of  attaining  good  beam  modulation. 
Hence,  these  two  effects  must  be  traded  against  each  other. 
When  the  beam  is  strongly  magnetized,  the  energy -dependence 
of  the  growth  is  approximately  a  y-1/  and  when  it  is  not 
magnetized,  a ^  y-1/3.  Since  the  criterion  for  strong 
magnetization  becomes  more  difficult  to  meet  as  X  decreases, 
short  wavelength  operation  will  probably  require  unmagne¬ 
tized,  or  at  least  weakly  magnetized,  electron  beams. 

The  current  available  from  the  variety  of  generators 

listed  in  Table  1  also  covers  a  wide  range.  When  the  beam 

1/3 

is  cold,  the  gain  will  scale  as  ( ,  and  thus  a  beam 
2 

with  17q  a/ cm  will  make  this  parameter  0.1.  We  have  seen 
earlier  that  if  this  is  used  with  reasonable  values  of  the 
other  parameters,  the  gain  will  be  in  the  .1  -  .5  cm.  1 
range.  The  first  and  third  entries  in  the  table  can 
probably  achieve  values  in  this  vicinity,  while  the  second 
and  fourth  entries  could  do  so  without  question.  The 
greater  current  available  from  the  second  and  fourth  type 
of  generator  might  also  make  up  for  dif f iciencies  in  another 
parameter . 

Field  emission  diode  generators  produce  very  large 
currents.  Hence,  they  are  in  principle  capable  of  producing 
a  lot  of  gain.  It  was  partly  for  this  reason  that  a 
generator  of  this  type  was  used  in  early  experiments 
designed  to  demonstrate  the  utility  of  stimulated  Cerenkov 
radiation.  They  also  have  the  ability  to  produce  beams 
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whose  energy  is  sufficient  to  couple  well  into  the  sub¬ 
millimeter  region.  Their  drawback  for  short  wavlength 
operation  may,  ironically,  be  the  fact  that  the  current  is 
large.  This  is  because  the  self-fields,  which  have  been 
neglected  in  our  analysis  of  gain,  may  lead  to  larger 
Ay/y,  and  hence  to  a  limit  of  the  usable  range  of  wavelength. 

The  accelerators  listed  in  the  table  will  also  be 
capable  of  operating  at  current  densities  which  give  usable 
gain.  The  peak  current  will  be  low,  but  the  focussing  could 
be  better.  The  current  in  a  linear  accelerator  will  also 
have  a  complicated  time  structure  (the  typical  value  of  I 
is  the  peak  in  the  micropulse) ,  and  this  will  complicate 
the  gain  calculations.  If  used,  however,  such  a  generator 
will  be  designed  to  work  at  short  wavelengths  and  the 
microstructure  may  be  a  minor  feature.  Both  this  compli¬ 
cation  and  the  role  of  self-fields  are  worth  further 
analysis. 

The  pulse  length  and  peak  power  entries  in  the  table 
are  largely  self-explanatory,  although  one  consequence  of 
the  pulse  versus  continuous  operation  is  worthy  of  comment. 

We  will  see  below  that  when  the  Q  of  a  cavity  is  reasonable, 
the  current  density  required  to  initiate  oscillation  will  be 
quite  modest.  Thus,  systems  with  relatively  low  gain, 
a  -  0.01  to  0.1,  may  be  very  usable  as  oscillators  well  into 
the  infrared,  whereas  the  same  beam  generator  would  not  be  a 
suitable  source  for  an  amplifier.  Overall,  it  is  to  be 
expected  that  the  pulsed  beam  generators,  due  primarily  to 
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the  larger  current  densities  available,  could  be  used  as 
both  oscillators  and  amplifiers  while  the  steady-state 
generators  would  be  largely  restricted  to  application  as 
oscillators. 

In  addition  to  its  relation  to  the  gain  of  the  device, 
the  current  density  plays  a  role  in  determining  the  wave¬ 
length  of  which  the  beam  may  no  longer  be  regarded  as  "cold”. 
Clearly,  if  all  parameters  other  than  current  remain  fixed, 
increasing  the  current  decreases  the  wavelength  for  which 
the  beam  must  be  regarded  as  warm.  The  current  density, 
the  energy  and  the  energy  spread  are  not  truly  independent, 
but  we  will,  for  the  purpose  of  discussion,  treat  them  as 
though  they  were. 

Restating  the  expression  kAv  =  uij/2  in  terms  of  the 
spatial  gain,  a (cm.  ,  and  a  beam  energy  spread  Ay,  we 
have  for  the  wavelength  at  which  the  cold/warm  transition 
occurs: 


It  is  clear  that  the  beam  generator  in  the  lower  energy 
end  of  the  range  considered  must  achieve  relatively  better 
energy  collimation  if  it  is  to  operate  in  the  "cold"  regime 
at  any  given  wavelength.  The  lower  energy,  pulsed  electro¬ 
static  devices  with  their  larger  currents  can  be  expected  to 
operate  in  this  mode  for  wavelengths  in  the  upper  submilli¬ 
meter-  to  one-millimeter  range.  As  the  beam  energy  rises, 
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Ill  B.  A  MILLIMETER  WAVELENGTH  EXPERIMENT 

Shown  in  Fig.  16a  is  the  schematic  diagram  of  an 
apparatus  which  has  been  constructed  in  order  to  test  a 
number  of  the  concepts  discussed  in  Section  II.  It 
consists  of  a  high-voltage  pulse  transformer,  a  thermionic 
cathode,  a  single-gap  accelerating  stage,  a  static  magnetic 
field  to  guide  the  beam,  and  a  tubular  quartz  -  dielectric 
resonator.  The  beam  is  guided  along  the  axis  of  the 
resonator  and  is  collected  on  the  surface  of  a  plane  mirror 
placed  at  a  45-degree  angle  to  the  end  of  the  copper  tube 
which  supports  the  quartz.  The  output  is  then  focussed 
into  a  plane-parallel  Fabry-Perot  resonator. 

The  transformer  has  a  peak  capacity  of  400kv  at  200A 
for  10y  seconds  at  a  300hz  repetition  rate.  However, 
typical  operating  conditions  for  the  present  work  are 
100-300kv  at  0.5-20A  for  approximately  ly  second.  At  the 
present  time,  the  pulse  repetition  frequency  is  a  few  hz. 
The  primary  of  the  transformer  is  driven  by  a  hydrogen- 
thyratron  triggered,  single-stage,  pulse-forming  network. 

A  characteristic  high-voltage  trace  and  an  accompanying 
millimeter-wavelength  radiation  pulse  is  shown  in  Fig.  16b. 
The  voltage  range  of  the  generator,  which  provides 
adequate  power  for  stimulated  Cerenkov  radiation  experi¬ 
ments,  is  such  that  millimeter-submillimeter  operation  can 
be  expected. 
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Coherent  output  radiation  has  been  obtained  at  wave¬ 
lengths  extending  from  about  1  centimeter  to  below  1.5 
millimeters.  The  wavelength  of  the  radiation  depends  upon 
the  guide  radius,  the  relative  amount  of  dielectric  and 
its  dielectric  constant,  and  the  beam  voltage.  Configu¬ 
rations  which  should  work  in  the  fundamental  mode  over 
the  1  centimeter  -  3  millimeter  range  have  been  studied, 
and  reasonable  agreement  with  the  theory  of  Section  II  is 
found.  The  diameter  of  the  copper  guide  which  supports  the 
quartz  tube  is  approximately  1.5  centimeter,  and  tubes  with 
1-3  millimeter  wall  thicknesses  have  been  used  most 
recently.  Thus  output  wavelengths  which  are  less  than 
the  transverse  dimension  of  the  waveguide  have  already  been 
obtained. 

At  the  longer  wavelengths,  the  frequency  has  been 
determined  with  the  Fabry-Perot  interferometer,  and  some 
typical  data  illustrating  the  behavior  is  shown  in  Fig.  17. 
In  Fig.  17a,  a  calibration  trace  made  with  a  35  hz  Gunn 
diode  source  is  displayed,  and  in  Fig.  17b,  experimental 
data  with  approximately  the  same  wavelength  is  shown. 

The  trace  shown  in  Fig.  17b  consists  of  many  repetitions 
of  the  e-beam  pulser,  and  it  shows  that  the  average 
spectral  width, which  is  itself  quite  narrow  (.1-1  percent), 
is  primarily  due  to  shot-to-shot  reproducibility.  The 
output  of  a  single  pulse  is  apparently  very  coherent.  An 
interferometer  for  shorter  wavelengths  is  under  construc¬ 
tion*  The  shorter  wavelengths  have  been  measured  with 


cutoff  filters. 


The  output  is  monitored  by  ordinary  microwave  diodes, 
either  IN23's,  IN26's,  or  lN53's,  with  the  latter  being 
used  down  to  wavelengths  below  1.5  millimeters.  Attenuation 
levels  of  30-60  db  are  required  in  order  to  insure  that  the 
output  levels  are  below  the  burnout  levels  of  the  diodes. 

The  absolute  output  power  levels  are  not  yet  precisely 
determined  except  within  an  experimental  uncertainty  which 
becomes  greater  at  shorter  wavelengths  and  varies  between 
.1  and  5  percent  of  the  beam  power.  All  of  the  factors 
which  control  this  conversion  efficiency  are  not  yet  well 
understood.  The  loaded  Q  of  the  resonator  is  modest,  and 
the  system  is  probably  operating  as  a  super- radiant  oscil¬ 
lator.  If  this  is  so  the  previously-stated  conversion 
efficiencies  are  plausible,  but  as  is  the  case  with  the  theory, 
nonlinear  behavior  of  the  experimental  device  is  practically 
unexplored . 

The  experiment  described  above  is  in  a  relatively  early 
stage  of  development.  It  does  appear,  though,  that  milli- 
meter-submillimeter  Cerenkov  devices  are  a  realistic 
possibility. 
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iii  c.  Cerenkov  devices  in  the  short  wavelength  limit 

The  bulk  of  the  analysis  presented  in  the  preceding 

sections  was  devoted  to  systems  which  would  work  in  a  way 

similar  to  microwave  devices  of  the  traveling  wave  type. 

We  emphasized  the  region  which  was  collective  where  there 
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was  an  interaction  of  fields  with  the  beam.  Recently  , 
however,  radiation  with  a  wavelength  of  about  3p  was 
obtair.od  from  an  electron  beam  produced  by  a  linear 
accelerator  and  a  coupler,  which  was  a  helical,  static, 
magnetic  field.  This  device  operated  in  the  so-called 
single- particle  limit: 

WpL/cSY3/2  <  <  1  (66) 

In  this  expression,  L  is  the  length  of  the  interaction  region, 
and  the  other  quantities  are  defined  in  earlier  sections. 

The  interaction  in  this  free-electron  laser  depends 
upon  the  coupling  between  the  tranverse  velocity  modulation 
imparted  to  the  beam  by  the  helical  pump,  and  the  magnetic 
field  component  of  the  wave  in  question.  The  beat  wave 
produced,  which  is  sometimes  called  a  ponderomotive  wave, 
travels  at  a  phase  velocity.  Hence  it,  too,  is  a  slow  wave 
interaction. 

There  is  no  reason  why  the  slow  wave  supported  by  a 
thin  film  dielectric  wave  guide  could  not  also  be  used  in 
the  ly  range.  Furthermore,  the  beam  quality  requirements 
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are  the  same  if 


kb/y  ~  1  (67) 

where  b  is  the  beam  thickness.  The  electric  field  asso¬ 
ciated  with  the  slow  wave  structure  can  be  greater  that  that 
of  the  ponderomotive  wave.  A  detailed  comparison  of  the  two 
modes  of  interaction  (dielectric  versus  static  B  pump)  would 
be  too  lengthy  to  present  here,  but  in  view  of  the  above 
comment,  it  follows  that  there  are  beam-resonator  combina¬ 
tions  for  which  the  Cerenkov  device  gain  will  exceed  the 
gain  of  the  ponderomotive  wave-based  free-electron  laser 
(FEL) .  In  this  case,  the  behavior  of  a  Cerenkov  device  will 
be  more  "forgiving"  than  the  FEL.  An  example  of  this  is  the 
relaxation  of  the  energy  collimation  required  if  a  beam  in  the 
single-particle  limit  is  to  be  considered  monoenergetic.  Simple 
consideration  of  the  phase  shift  accumulated  in  a  length  L 
leads  to  the  requirement: 

ck  AB  L/v  <  tt  (68a) 


or 


A 

L 


> 


(68b) 


Higher  gain  would  allow  a  shorter  length,  and  hence  a 
system  which  was  in  the  "cold"  region  at  a  shorter  wave¬ 
length.  Due  to  this  and  to  other  considerations,  it  is  of 
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interest  to  make  a  detailed  comparison  of  these  two  systems. 
We  will  conclude  this  section  with  a  brief  discussion  of 
two  possible  configurations  which  might  be  used  in  the  lu 
range.  The  first  is  shown  in  Fig.  18. 

Shown  in  the  figure  is  a  thin  dielectric  guide  with 
tapered  ends.  An  electron  beam  passes  over  the  guide  and 
couples  to  the  evanescent  field.  When  the  beam  has  a  high 
enough  energy,  the  coupling  will  be  good.  The  tapered  ends 
face  mirrors  which,  together  with  the  guide,  form  an  optical 
resonator. 

It  might  also  be  useful  to  taper  the  thickness  of  the 
guide  near  its  ends.  By  doing  this,  the  field  energy  in 
the  guide  can  be  increased  at  the  expense  of  that  stored 
above  the  guide.  If  the  taper  is  adiabatic  on  the  scale  of 
k  1 ,  the  dispersion  relation  will  vary  smoothly  as  will  the 
field  distribution.  Shown  in  Fig.  19  are  the  results  of 
one  possible  experimental  configuration  including  a  taper. 

As  the  guide  thickens,  the  field  distribution  is  pulled 
down  into  the  guide  and  formed  into  a  half-sinusoid.  The 
latter  form,  which  is  closed  to  the  normal  mode  distribution 
of  a  guide  which  is  covered  top  and  bottom,  should  give 
better  control  of  the  input-output  coupling. 

Another  possible  version  of  a  short  wavelength  system 
is  shown  in  Fig.  20.  In  this  system,  the  output  optics  are 
placed  below  the  beam.  Incident  radiation  comes  through 
the  end  face  of  the  guide  and  normal  to  the  face,  but  the 
angle  formed  with  the  top  surface  is  greater  than  that 
required  for  total  internal  reflection.  This  guide  is  also 


tapered.  As  it  thins  out,  the  field  is  pushed  out  to  increase 
the  coupling.  It  is  then  pulled  in  again  and  passed  to  a 
second  mirror. 

The  two  configurations  just  discussed  are  highly 
schematic.  The  beam-to-field  coupling  will  be  a  straight¬ 
forward  matter,  but  the  overall  optical  system  may  be  quite 
complex.  Self-reproducing  patterns  of  the  general  type 
discussed  should  exist,  however,  and  well- collimated 
relativistic  beams,  together  with  the  resonator,  could  form 
the  basis  of  far-infrared  Cerenkov  devices. 


IV.  CONCLUSION 


The  theme  of  the  preceding  three  sections  is  the  use 
of  the  combination  of  the  electron  beam  and  the  dielectric 
resonator  to  generate  coherent  radiation  in  the  millimeter 
through  infrared  portions  of  the  electromagnetic  spectrum. 

At  the  present  time,  there  is  experimental  confirmation  of 
some  of  the  basic  theoretical  assumptions  in  the  10-millimeter 
to  1-millimeter  region.  The  theoretical  advantages  of  smooth¬ 
bore,  slow  wave  structures  at  shorter  wavelengths  indicate 
that  submillimeter-infrared  operation  is  also  a  realistic 
possibility.  Beams  with  the  monoenergicity  and  current 
sufficient  to  achieve  oscillation  for  moderate  Q  resonators 
((a^Q/co)  >  1)  currently  exist.  The  long-term  interest  in 
these  sources  will  not  depend  upon  whether  they  work  at  all, 
but  how  they  might  compare  to  other  possible  sources. 
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Cerenkov  maser  components. 

Alternate  configurations  of  the  basic 
device,  a)  slab  waveguide;  b)  possible 
Cerenkov  gas  laser. 

Geometry  for  the  plane-wave  -  e-beam 
interaction  in  a  dielectric  medium. 

The  relations  w  vs  k  for  two  different 
perpendicular  components  of  wavenumber 
(Pi  >  P2>  . 

The  complete  dispersion  relation  for  the 
e-beam  dielectric  system.  The  region  of 
instability. 

The  components  of  growth  rate,  a)  the 
dielectric  contribution;  b)  the  beam  thresh- 
hold  function  F. 

Dispersion  relation  for  the  free  waves  of  a 
partially-filled  guide. 

The  radial  dependent  of  the  axial  electric 
field,  a)  the  fast  wave  region;  b)  the  slow 
wave  region. 
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The  basic  slab  guide  geometry  and  axial 
field  dependence. 

The  dispersion  relation  for  the  slab  guide. 

Gain  curve  shape  vs.  beam  velocity  for  the 
beam-slab  guide  system. 

Numerical  values  of  dispersion  (a)  and 
spatial  gain  (b)  vs.  beam  velocity  for  a 
.025  cm.  thick  quartz  guide. 

Geometry  of  slab  guide/beam  system  when 
there  is  a  finite  gap  between  the  beam  and 
the  guide. 

The  relative  positions  of  the  phase  velocity 
and  the  beam  velocity  distribution  function 
(F8 )  when  the  beam  is  "cold"  and  when  it  is 
"warm" . 

Qualitative  shape  of  the  dispersion  and  gain 
in  the  "warm"  beam  limit. 

The  configuration  of  an  experimental  device 
designed  to  prodce  mm-wavelength  stimulated 
Cerenkov  radiation,  a)  the  device; 
b)  typical  voltage  and  radiation  pulse. 

Fabry-Perot  interferometer  output, 
a)  calibration  (35  Ghz  source);  b)  stimulated 
Cerenkov  radiation. 


Figure  18:  Slab  guide  resonator  configuration  for  a 
possible  far-infrared  device. 

Figure  19:  The  effect  of  tapering  the  slab  guide 
thickness,  a)  the  guide;  b)  the  field 
distribution. 

Figure  20:  Another  possible  configuration  for  short 
wavelength  operation. 
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TUNABLE  CERENKOV  LASERS* 


J.E.  Walsh  and  J.B.  Murphy 

Department  of  Physics  and  Astronomy 
Dartmouth  College 
Hanover,  N.H.  03755 


Abstract 

The  concept  of  a  tunable  Cerenkov  laser  is 
introduced  and  expressions  for  the  operating  wavelength 
and  stimulated  emission  rate  are  developed.  Emphasis 
is  on  the  development  of  criteria  which  must  be  met  if 
tunable,  far-infrared  device  operation  is  to  be  achieved. 


I.  INTRODUCTION 


A  Cerenkov  laser  consist  of  a  mildly  to  moderately 
relativistic  electron  beam,  a  dielectric  resonator,  and 
output  coupling  components.  One  possible  configuration  of 
such  a  device,  which  is  based  in  part  on  concepts  [1]  which 
have  been  experimentally  demonstrated  [2] ,  [3] ,  [4]  in  the 
middle-  to  lower-mm  portion  of  the  spectrum,  is  shown  in 
Fig.  1.  In  Fig.  1,  an  electron  beam  passes  over  a  thin 
dielectri ~  guide  and  interacts  (primarily)  with  the  axial 
component  of  the  TM  modes  of  the  structure.  A  bunching  of 
the  beam  will  occur  and  hence  stimulated  emission  will  be 
produced  in  those  modes  which  have  a  phase  velocity  slightly 
less  than  the  velocity  of  the  beam.  In  the  example  chosen, 
waves  are  launched  from  the  slab  waveguide  by  a  taper.  The 
cavity  is  then  completed  by  mirrors.  Other  configurations 
would  also  be  potentially  useful. 

A  general  result  of  the  analysis  given  in  subsequent 

sections  is  that  a  typical  transverse  scale  length  d  ,  wave- 

2 

length  A  and  relative  beam  energy,  y  =  e/mc  ,  are  related  by 
Ay  ^2Trd.  When  the  scale  length  is  some  fraction  of  a  mm  , 
y  in  the  2  to  10  range  implies  device  operation  in  the  far- 
infrared  portion  of  the  spectrum.  Other  devices  [5]-[8], 
particularly  those  which  use  a  periodic  magnetic  field  pump 
to  couple  an  electron  beam  with  a  radiation  field  are  possible 


tunable  infra-red  sources.  In  a  certain  fundamental  sense, 
all  electron  beam-driven  sources  are  similar,  and  some  very 
brief  comparison  of  Cerenkov  and  other  coupling  schemes  will 
be  Uncluded  in  the  last  section. 

Finally,  before  proceeding  to  the  detailed  analysis, 
it  should  be  pointed  out  that  the  notion  that  Cerenkov 
radiation  could  be  used  as  the  basis  of  a  practical  radiation 
source  is  not  new  [9].  Several  groups  [10] ,  [11] ,  [12]  have 
studied  spontaneous  Cerenkov  emission  and  some  have  used 
prebunched  beams  to  enhance  the  otherwise  weak  spontaneous 
emission  level.  In  addition,  the  general  idea  that  stimulated 
Cerenkov  radiation  could  be  used  to  obtain  laser  action  has 
been  mentioned  [13] ,  [14] ,  [15]  as  a  theoretical  possibility. 

The  first  clear-cut  demonstration  of  the  stimulated  Cerenkov 
radiation  prcess  does  seem  to  be  the  work  reported  in 
Refs.  [2] , [3] ,  and  [4], 

The  remainder  of  this  paper  is  divided  into  four 
sections.  These  deal  with:  the  development  of  an  expression 
for  the  linearized  gain,  a  summary  of  those  features  of  the 
dielectric  slab  guide  behavior  which  would  be  important  in 
Cerenkov  laser  operation,  a  discussion  of  the  implications  of 
the  results  of  the  preceding  two  sections  for  far-infrared 
device  performance,  and  finally  a  few  brief  concluding  remarks. 


II.  THE  STIMULATED  CERENKOV  EMISSION  RATE 


In  this  section,  an  expression  for  the  gain  of  a 
Cerenkov  Maser  will  be  deduced.  The  calculation  is 
facilitated  by  making  a  number  of  simplifying  assumptions: 

1)  that  the  beam  is  monoenergetic;  2)  that  the  beam  density 
is  low  enough  so  that  space  charge  waves  may  be  neglected; 
and  3)  that  the  gain  in  a  single  pass  of  the  resonator  is 
small.  Numerical  criteria  for  these  will  be  established  and 
the  way  in  which  they  limit  possible  device  performance  will 
be  discussed  in  Section  IV. 

The  energy  gained  (d£/dt)  by  the  fields  as  the 
electron  beam  traverses  the  resonator  may  be  expressed  as: 


d 1 

dt  " 


(1) 


J  beam 

where  j  is  the  modulated  current  produced  on  the  beam  by 
the  fields  and  the  integration  covers  the  volume  of  the  beam. 
It  is  most  convenient  to  compute  the  current  with  the  aid  of 
the  linearized  Vlasov  equation.  This  may  be  written  in  the 
form: 


d6f 

dt 


V  f  (P) 
p  o' 


(2) 


Appearing  in  Eq.  (2)  are  5f  and  fQ  ,  the  perturbed  and 
average  components  of  the  beam  momentum  distribution  function. 
The  term  p  is  the  momentum  of  an  electron.  The  expression 


d/dt  on  the  lUS  of  Eq.  (2)  is  the  total  derivative  with 
respect  to  time,  along  a  zero  order  orbit  (one  unperturbed 
by  the  fields) .  If  we  now  assume  further  that  the  coupling 
is  in  the  axial  direction  and  with  a  traveling  wave  moving 
in  the  positive  z  direction 
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we  have 
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This  form  of  the  solution  of 
condition  Sf(z,p,t)  z=Q  * 
resonator  region  unmodulated. 


Eq.  (2)  satisfies  the  boundary 
0  ,  i.e.r  the  beam  enters  the 

The  current  is  given  by: 


j 


5f  dp 


(5) 


where  n  is  the  beam  density. 

The  results  of  Eqs.  (4)  and  (5)  may  be  substituted 
into  Eq.  (1)  and  an  expression  for  the  rate  at  which  energy 
is  added  to  the  resonator  obtained.  It  is  also  convenient 
at  this  point  to  perform  the  axial  portion  of  the  spatial 
integral.  When  this  is  done,  Eq.  (1)  becomes: 
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where 


0  =  (kv  -  cd)  L/2v  (7) 

is  the  relative  phase  change  "seen"  by  a  beam  electron 
traversing  a  resonator  of  length  L  .  Equation  (6)  may  be 
further  reduced  by  re-writing  the  momentum  integral  with 
the  aid  of  the  transformation: 


3  =  d  0  _3_ 

3PZ  dpz  &T  30 


Introducing  the  beam  distribution  in  phase  coordinates  F(0) 
and  performing  an  integration  by  parts.  Eq.  (6)  becomes: 
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z ' 

beam 


d9F  ( 0 ) 
Y3(9)02(0) 


2 

3  sin  0 

se  -pr 


(9) 


We  can  now  make  use  of  our  assumption  that  the  beam  distribution 
is  monoenergetic  and  perform  the  phase  integral.  Before  making 
this  straightforward  substitution,  however,  it  will  also  be 
useful  to  introduce  some  further  definitions. 

Let 


n 


dA  n  |E  I 
1  z  1 


2 


(10) 


where  is  the  beam  area  and  n  is  an  effective,  wavelength- 

2 

dependent,  density.  The  normalizing  field  strength  |e  | 
is  chosen  as  some  convenient  point  (e.g.,  the  surface  of  the 
dielectric  in  the  present  case) . 


G'(9)  =  Jq  ■ 2~-  (12b) 

is  the  derivative  of  the  diffraction  lineshape.  in  writing 

Eq.  (12)  we  have  also  identified  the  total  beam  current, 

I.  =  -n  ecB  A,  ,  and  an  effective  current,  I  =  ec/r 
boon  o  o 

2  2 

(rQ  =  e  /me  )  which  appears  naturally  in  the  problem. 

A  sketch  of  Qb  1  vs  0  is  shown  in  Fig.  2.  Since 
Qb  ^  <  0  when  0  <  6  <  -n,  the  system  will  have  gain  in  this 
region.  When  |Q^|  becomes  less  than  the  Q  of  the  resonator, 
oscillation  threshold  is  reached.  Also  shown  in  Fig.  2  is  the 
electron  beam  velocity  distribution  which  has  been  plotted  as 
a  function  of  the  relative  phase.  A  beam  with  relative  phase 
spread 


(A0)b  <  TT 


(13) 


may  be  regarded  as  monoenergetic.  All  the  electrons  in  such  a 


beam  can  contribute  to  gain.  This  is  the  condition  illustrated 
in  Fig.  2.  An  alternative  and  more  useful  expression  for  this 
condition  may  be  obtained  when  Eq.  (7)  is  used  together  with 
Eq.  (13)  to  relate  the  beam  energy  spread  6y  to  the  wavelength 
for  which  inequality  (13)  is  just  satisfied.  This  is: 

X  >  L  -6y  a  (14) 

(Byr 

Discussion  of  the  magnitude  of  ^  and  the  limiting  wavelength 
will  be  given  after  a  particular  resonator  is  introduced. 


III.  THE  SLAB  WAVEGUIDE  RESONATOR. 


The  dielectric  slab  waveguide  [16]  is  one  possible 
choice  of  the  wave- supporting  structure  for  a  Cerenkov  laser. 
This  choice  may  not  offer  the  highest  possible  gain,  but  it 
would  allow  relatively  straightforward  independent  adjustment 
of  the  beam  and  resonator  axes.  Thus  it  is  an  experimentally 
interesting  choice,  and  we  will  consider  it  in  some  detail. 

The  basic  geometry,  shown  in  Fig.  3,  consists  of  a  beam 
of  thickness  b  ,  at  a  distance  d  above  the  surface  of  the 
guide.  We  assume  that  propagation  occurs  in  the  z  direction 
and  that  all  quantities  are  uniform  in  y  ,  the  direction 
perpendicular  to  both  the  normal-to-the-guide  surface  and  the 
direction  of  propagation.  Later,  in  order  to  facilitate 
numerical  evaluation,  we  will  assume  that  the  beam  and  mode 
widths  in  the  y  direction  are  approximately  b  and  am 
respectively. 

The  electric  fields  associated  with  the  TM  modes  of  the 
structure  shown  in  Fig.  3  are: 
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E  e-qx  eilfcz-wt) 

zo 


ik  5E* 

7  ^ 


in  the  region  x  >  a  .  In  the  region  x  <  a 


(15a) 

(15b) 


Ez  =  Ezl  sinpx  (16a) 

and 

Ev  -  A*  (16b) 

7  tx- 

- * —  ^  ^ . . | ^ 


The  terms : 


and 


2 

w  e 


(17a  ) 


(17b) 


are  the  squares  of  the  transverse  wavenumber  in  the  dielectric 
and  in  the  vacuum  regions  respectively.  The  relative 
dielectric  constant  of  the  guide  materials  is  denoted  by  e. 

The  fields,  together  with  the  boundary  conditions,  can 
be  used  to  obtain  a  dispersion  function: 


D(k,to)  =  p  Ta^Pa.  -  q  (18) 

whose  roots  D  =  0  determine  the  relation  between  w  and  k 
for  the  TM  modes  of  the  structure.  Typical  results,  plotted 
in  units  of  a/X  =  toa/2irc  vs  a/X^  =  ka/2iT  are  shown  in 
Fig.  4.  The  lowest  mode,  TMq^,  has  no  lower  frequency 
cutoff  while  all  higher  order  modes  have  a  long  wavelength 
cutoff  given  by: 


The  trajectories  a/ A  vs  a/Ag  lie  in  the  region  of  the 
plane  bounded  by  the  light  line  (slope  1)  and  the  light  line 
in  the  dielectric  material,  the  line  of  slope  1/e  ( e  =  n)  on 
Fig.  4.  Also  shown  on  Fig.  4  is  a  line  whose  slope  is  equal 
to  8,  the  velocity  of  a  beam  electron  relative  to  the  speed  of 
light.  When  8  >  1/n  the  beam  satisfies  the  Cerenkov  threshold 
condition; 


Y  >  yT  (20a) 

where 

YT2  =  e/(l-e)  (20b) 

=  1/(1  -  8t2)  (20c) 

and 

6T2  =  1/e  ( 20d) 

Velocity  synchronism  is  defined  as  that  point  where  the  phase 
velocity  of  one  of  the  guide  modes  is  equal  to  8  • 

Examination  of  Eq.  (18)  indicates  that,  for  a  given  e, 
the  location  of  the  dispersion  curves  is  a  function  of  the 
guide  thickness  alone.  Thus  we  arrive  at  a  basic  scaling  of 
a  Cerenkov  device  of  this  type;  A  a  .  As  the  guide  thickness 
is  decreased,  the  wavelength  at  which  phase  velocity  synchronism 
occurs  decreases  as  well. 

We  can  now  also  obtain  an  expression  for  n/nQ  ,  the 
effective  beam  density  defined  in  Section  II.  If  we  imagine 


that  we  vary  the  beam  velocity  from  B  =  1  down  to  $  =  we 
sweep  the  synchronous  point  along  the  dispersion  curve.  Thus 
the  output  is  voltage  tunable.  However,  at  the  synchronous 
point  the  transverse  wavenumber  in  the  vacuum  becomes: 

q  =  k/y  (21a) 

or,  since  to  =  ckB  at  synchronism, 

q  =  2  it/ X  By  (21b) 


The  field  decreases  more  rapidly  in  the  x  direction  as  the 
beam  velocity  decreases.  The  coupling  becomes  weaker.  The 
effective  beam  density  for  the  flat  slab  resonator  is  thus 
given  by: 


n 


e-2qd 


'22) 


We  have  assumed  in  deriving  Eq.  (22)  that  the  actual  density 
distribution  of  the  beam  is  approximately  rectangular  in  the 
x  direction.  The  first  factor  in  Eq.  (22)  is  a  decoupling, 
and  the  second  factor  represents  the  fraction  of  the  beam 
seen  by  the  field. 

The  fields  given  by  Eqs.  (15)  and  (16)  may,  together 
with  the  usual  definition  of  energy,  be  used  to  obtain  an 
expression  for  the  energy  stored, £ ,  and  hence  also  for 
|ezq| 2  in  terms  of  £ .  Straightforward  calculation  leads  to: 


IV.  DISCUSSION 


The  results  of  the  two  preceding  sections  may  be  used 
to  infer  some  general  features  of  the  gain  cf  a  Cerenkov 
laser  and  hence  also  to  make  estimates  of  the  parameters 
required  to  achieve  oscillation  threshold  at  a  given  wave¬ 
length.  The  full  expression  for  the  relative  gain  of  a  slab 
guide  Cerenkov  device  is: 


Qb 


1 

(By)5 


(l-YT2/Y2)ka/y 

1  +  ^(-e+1) 


G'  (6) 


(24) 


Shown  in  Fig.  5  is  a  sketch  of  the  TMQ2  dispersion  curve 
together  with  a  sketch  of  the  general  behavior  of  Qb_1(V2). 
The  latter  is  plotted  vs  Y_1  and  is  displayed  in  a  way  that 
emphasizes  how  the  gain  behaves  as  the  beam  velocity  is 
varied  between  B  =  1  and  3  =  BT.  Near  the  light  line  the 
gain  is  small  and  the  behavior  is  dominated  by  the  (By) 
term.  As  B  moves  away  from  the  light  line  the  gain  increases 
until  it  reaches  a  maximum  determined  by  the  exponential 
factors  in  n/nQ.  Further  decreases  in  B  results  in  a  rapid 
decreases  of  Qb  as  3  BT»  Because  of  the  (By)  5  term  the 
gain  will  peak  when  the  arguments  of  the  exponentials  which 
appear  in  n/nQ  are  somewhat  greater  than  unity. 

It  is  also  convenient  at  this  point  to  mention  the 
effect  of  finite  energy  spread  in  the  electron  beam.  The 
gain  expression  displayed  in  Eq.  (24)  holds  when  the  beam 
energy  spread  is  small  enough  to  allow  the  bulk  of  the  beam 
distribution  to  fall  within  one  loop  of  G'(0).  When  this  is 


not  the  case  some  beam  electrons  hinder  wave  growth.  In  order 
to  understand  the  effect  of  this  on  the  gain  it  is  easiest  to 
consider  the  gain  expression  when  the  integration  by  parts  on 
P  (or  6)  is  not  performed,  and  G(0)  can  be  regarded  as  narrow 
compared  to  F(6).  The  gain  is  now  proportional  to  F'(6)  (or 
F '  (P ) ) .  Furthermore,  the  Jacobian  (dP/d6)  will  no  longer 
appear  and  hence  we  lose  a  factor  of  to.  Thus  in  the  warm 
beam  limit: 


Q, 


-1 


to 


-1 


(25) 


and  gain  will  decrease  with  decreasing  wavelength. 

Holding  all  factors  save  By  constant,  a  relative 
maximum  in  Eq.  (24)  occurs  when 


A  Ci 


2nd 

33y 


(26) 


We  anticipate  the  use  of  modest  By  and  thus  we  expect  the 
cold-warm  transition  wavelength  defined  by  Eq.  (14)  will  be 
longer  than  A  determined  by  Eq.  (26) .  If  this  is  the  case, 
the  gain  will  reach  a  maximum  at  a  wavelength  longer  than 
that  at  the  maximum  for  a  cold  beam.  When  this  condition  applies 
the  general  behavior  is  shown  by  the  wavy  line  on  1  in  Fig.  5 
In  order  to  make  a  quantitative  estimate  of  the  magnitude 
of  Qb  we  make  further  specific  choices.  We  assume  that  the 
warm-cold  transition  wavelength  is  about  twice  the  wavelength 
given  by  Eq.  (26)  and  that  we  will  examine  the  possibility  of 
operating  at  the  still  longer  wavelength: 


I 


This  last  choice  is  equivalent  to  setting  2qb  =  2  and  hence 
to  assuming  that  a  transverse  scale  length  adequate  for  the 
maintenance  of  good  coupling  is  assured.  An  expression  for 
the  size  of  1  at  any  point  then  takes  the  following  much 
simpler  form: 


\2vb) 


5  L2 


(l-e~2)  -d/b 
2 


The  factor,  .54,  arises  from  setting  G'(e)  =G'(e)  .  In 

max 

addition  to  the  desired  operating  wavelength,  variables 
remaining  in  Eq.  (28)  are  those  pertaining  to  the  production, 
focussing  and  guiding  of  the  electron  beam.  Before  addressing 
these,  however,  it  will  be  convenient  to  briefly  examine 
cavity  and  output  coupling  losses.  We  will  thereby  establish 
the  approximate  magnitude  of  1  needed  for  threshold. 

We  assume  that  two  losses  predominate;  the  first  is 
absorption  in  the  guide  material,  and  the  second  is  the  output 
coupling  loss.  The  latter,  expresses  as  a  loading  QL,  is 
given  approximately  by: 


(1-p) 

P 


where  p  is  a  reflection  coefficient.  If  we  assume  that  we 
have  no  interest  in  cavities  with  L  less  than  10  cm,  nor  in 


wavelengths  greater  than  1  mm,  relatively  modest  p  will 
place  (Ql  1)  in  the  10  4  range. 

The  general  expression  for  the  energy  lost  in  the 
dielectric  takes  the  form: 


d*D 

dt 


8n 


|E|2dv 

diel. 


(30) 


The  integral  can  be  performed  with  the  aid  of  expressions 
developed  in  Section  II.  Completing  this  and  defining  a 
dielectric  quality  factor  QD  we  have: 


_1_  2ob  /  A  \3  2 

Qd  VT  \27Tby  yT 


(31) 


In  expressing  Qq  1  we  have  also  substituted  for  By  where 
appropriate,  and  denoted  the  dielectric  absorption  coeffi¬ 
cient  as  a  (cm  3)  . 

Since  losses  in  the  submm  and  FIR  regions  are  not 
necessarily  small,  we  will  use  Eq.  (31)  to  establish  a 
condition  for  threshold: 


.  54Ve~  f  A L  \  ^  _  _I_  (1-e  3)  e~d/b  >  ^ 

aa  y  2  l  2Tib2 )  Io  2 
m'  T  \  ’ 


(32) 


2  2 

The  factors  involving  (AL/irb  )  and  I/I  will  show  the 


largest  variation  and  hence  they  will  be  singled  out  for 
scrutiny.  We  asssume,  somewhat  arbitrarily,  that  we  will  be  most 
concerned  with  devices  which  are  based  on  beam  sources  with 


an  available  current  between  1.75  ma  and  17.5  A,  and  in  the 

wavelength  range  lOOy  to  1  mm.  In  the  range  of  the  former 

-7  —  3 

variable  the  ratio  I/IQ  varies  from  10  to  10  . 

2  2 

The  magnitude  of  the  term  (AL/2fTb  )  can  be  easily 
inferred  if  we  re-express  it  in  terms  of  the  value  of  By 
required  by  coupling.  It  is: 


Electron  beam  focussing  such  that  b  ^  1-.5  mm  is  typical 

for  beams  associated  with  microwave  tubes,  and  hence  values 

2  3 

of  L/Byb  in  the  10  -10  range  are  quite  realistic.  This  is 

in  itself  sufficient  to  achieve  threshold  within  our  chosen 

current  range.  Some  further  help  in  attaining  threshold  will 

be  contributed  by  the  terms  in  Eq.  (32)  which  have  a  smaller 

variation.  The  value  of  am  will  be  generally  larger  than, 

but  comparable  to,  b  and  a  &  1  at  a  wavelength  of  lOOy 

can  be  obtained  in  several  relatively  low  index  materials  [17] 

At  longer  wavelengths  (A  1  ram)  values  of  a  as  low  as  2  x  10~ 

are  realistic  and  the  overall  contribution  to  threshold  by  the 

2 

material-related  terms  in  Eq.  (32)  could  be  as  large  as  10  . 

The  beam  energy  collimation  required  if  the  beam  is  to 
by  regarded  as  cold  was  given  first  in  Section  II.  Returning 

now  to  Eq.  (14)  we  see  that  A  =  100vi  and  L,  =  10  cm  would 

-3 

require  6y/By  —  10  if  By  —  1-  If  By  rises  to  about  three, 

-2 

the  requirement  drops  to  6y/By  !£  10  .  The  first  is  an 

excellent  but  achievable  value  and  the  second  is  quite 


straightforward.  Beam  velocity  spread  will  not  be  a  major 
obstacle  to  operation  in  the  100m  -  1  mm  wavelength  range. 

Finally,  it  is  possible  to  make  a  rough  estimate  of 
the  efficiency.  Since  a  change  of  the  beam  energy,  which 
is  sufficient  to  move  the  centroid  out  of  the  gain  region  in 
Fig.  2,  represents  an  approximate  non-linear  limit,  we  have 

6y  (3y) 3  1  (33) 

Y-l  Y“1  L 

or  (Y  +  1)  (3y  ) 2  £ 


as  the  approximate  efficiency. 


V.  CONCLUSIONS 


The  general  conclusion  which  can  be  drawn  from  the 
preceding  analysis  is  that  is  should  be  possible  to  construct 
tunable  Cerenkov  lasers  in  the  far-infrared  portion  of  the 
spectrum.  These  devices  are  beam-voltage-tunable  and  they 
could  in  principle  have  a  reasonable  efficiency.  Moderate- 
to  high-power  output  can  probably  be  obtained.  Like  other 
free-electron  lasers,  decrease  of  the  operating  wavelength 
is  achieved  by  increasing  beam  energy.  However,  the  energy 
wavelength  scaling  A8y  const,  is  such  that  the  energy 
required  to  obtain  far-infrared  operation  will  be  less  than 
that  required  in  a  magnetic  wiggler-coupled  device.  Hence 
Cerenkov  devices  may  be  relatively  compact. 
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FIGURE  CAPTIONS 

B 

1.  Possible  Cerenkov  laser  configuration,  showing  a 
dielectric  slab  guide,  a  beam,  and  mirrors. 

m  2.  The  relation  between  an  approximately  monoenergetic 

beam  distribution  function,  measured  in  relative 
phase  coordinates,  F(9),  and  the  derivative  of  the 
diffraction  line  shape,  G'  (8)  . 

3.  Geometric  details  of  a)  the  beam,  the  guide  and  the 
gap; ,  b)  the  relation  between  the  x  dependence  of 

""  the  field  and  the  preceding  quantities. 

4.  The  free  space  wavelength  vs  the  guide  wavelength, 
both  measured  in  units  of  the  guide  thickness. 

1  5 .  The  relation  between  dispersion-related  and  beam- 

bunching  quantities . 
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ABSTRACT 


The  interaction  of  a  magnetized  electron  beam  with  a  slow 
guided  electromagnetic  wave  has  been  studied  experimentally. 

An  electron  beam  driven,  dielectric-lined,  cylindrical 
waveguide  has  produced  30-100  kilowatts  of  coherent  radiation 
over  an  octave  band  on  the  TMQ1  mode  of  the  waveguide. 
Operation  on  the  TM^  mode  has  been  realized  as  well. 

Impedance  mismatches  at  the  ends  of  the  liner  section  provide 
a  reflected  signal,  which  undergoes  amplification  during 
successive  passes. 

The  growth  rate  of  the  interaction  is  calculated,  taking 
into  account  the  cylindrical  geometry  of  the  experiment,  for 
the  collective  regime  and  for  the  single  particle  regime.  An 
estimate  of  the  saturating  field  amplitude  leads  to  predicted 
powers  for  comparison  with  experiments. 
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CHAPTER  I 


1.1  INTRODUCTION 

The  search  for  new  sources  of  electromagnetic  radiation 
has  led  to  the  construction  of  many  novel  devices.  In  the 
visible  spectrum,  dye  lasers  have  achieved  tunability  while 
multi-stage  CO^  and  neodymium  lasers  have  achieved  high  power. 
In  attempts  to  combine  the  virtues  of  tunability  and  high 
power,  excited  states  have  been  stimulated  to  emit  one  photon 

at  a  widely  tunable  wavelength,  while  emission  in  one  case  of 

1  2 
another  photon,  and  in  another  case  a  phonon,  completes  the 

transition  between  states. 

At  millimeter  and  submillimeter  wavelengths  a  similar  but 
more  hopeful  situation  exists.  Tunable  devices  are  typically 
low  power.  However,  the  free  electron  laser  possesses  both 
tunability  and  high  power.  It  has  so  far  been  an  experimental 
device,  but  a  12  kilowatt  user  facility  is  now  under 

•3  il 

development  at  UCSB.  ’ 

The  Cerenkov  maser,  considered  here,  is  another  tunable, 
high  power,  millimeter  source,  which  may  have  advantages  over 
the  free  electron  laser  in  certain  situations.  The  tunability 
and  power  capability  have  been  tested  experimentally  and  are 


found  to  confirm  predictions.  Significant  power  at  three 
millimeters  has  been  observed. 


rc 

« 

* 

:  • 
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The  Cerenkov  maser  is  one  of  a  class  of  schemes  for 
extracting  electromagnetic  energy  from  electron  beams.  In  the 
Cerenkov  scheme,  a  dielectric,  combined  with  a  guiding 
structure,  is  used  to  support  a  slow  electromagnetic  wave  in 
the  vicinity  of  a  beam.  The  longitudinal  bunching  of  the  beam 
and  the  amplitude  of  the  electromagnetic  wave  drive  each  other 
to  large  amplitude. 

A  variety  of  other  schemes  exist  in  this  class.  The  free 
electron  laser  employs  a  spatially  modulated  transverse 
magnetic  field,  which  couples  the  beam  to  a  radiation  field. 

The  gyrotron  uses  a  carefully  designed  electron  gun  and 
guiding  magnetic  field  configuration  to  define  helical  orbits 
for  the  electrons  in  the  beam.  Bunching  of  the  beam  occurs 
along  the  helix  rather  than  along  the  longitudinal  axis, 
driving  the  amplitude  of  the  beam's  own  cyclotron  radiation. 

An  annular  beam  with  relatively  large  transverse  energy  is 
required . 

The  klystron,  the  magnetron,  and  the  traveling  wave  tube 
depend  on  resonant  cavity  structures  in  different  config¬ 
urations  to  bunch  the  beam.  Variations  of  these  devices  have 

been  scaled  to  millimeter  wavelengths  at  the  expense  of  power 
5 

and  efficiency  and  scaled  to  high  power  at  the  expense  of 
frequency. ^  The  frequency  is  ultimately  limited  by  the  ability 
to  construct  small  resonant  cavity  structures. 
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One  gyrotron  has  produced  1 . 1  megawatts  at  three 

7 

millimeters  in  a  100  microsecond  pulse.  Another  gyrotron, 
which  will  generate  five  millimeter,  high  power,  cw  radiation, 

O 

is  currently  under  development  by  Varian  Associates.  Gyrotron 
wavelengths  are  ultimately  limited  by  the  attainable  magnetic 
field,  since  output  is  at  the  electron  cyclotron  frequency. 

The  field,  B,  required  to  operate  at  wavelength,  x,  is  given 
by 

B  =  ( ymc^/e) /x ,  (1.1) 

thus  scaling  inversely  as  the  wavelength.  The  quantities  m, 

e,  c,  and  y  are  the  mass  and  charge  of  an  electron,  the  speed 

2  -1/2 

of  light,  and  (1-(v/c)  )  ,  respectively,  and  v  is  the  beam 

velocity . 

One  free  electron  laser  has  operated  at  a  wavelength  of 

Q 

3.4  microns  by  employing  a  highly  relativistic  beam.  The 
output  wavelength  is  given  by 

*  =  xp/(B(i+e)Y2) ,  ( 1 .2) 

B  -  v/c, 

where  Xp  is  the  period  of  the  transverse  magnetic  field 
modulation.^  The  pump  wavelengt’  ,  xp  was  3*2  cm  for  the  above 
experiment;  thus  y“90  was  required.  The  lower  limit  on  xp  is 
approximately  one  centimeter.  A  beam  energy  of  79  keV  is 

-  3  - 


required  to  give  x=x  ,  and  175  keV  gives  x=x  /2 .  Thus 

•  r 

relatively  modest  beam  energies  can  be  used  to  operate  a  free 

electron  laser  at  millimeter  wavelengths.  More  complete 

discussions  of  the  operation  and  history  of  millimeter  devices 

11-17 

are  available  and  contain  further  references. 

The  Cerenkov  scheme  imposes  less  severe  requirements  on 

structure  complexity,  magnetic  field  intensity,  and  electron 

beam  energy  than  the  above  schemes.  The  original  suggestion 

that  coherent  Cerenkov  radiation  might  be  obtained  from  a 

1  ft 

bunched  beam  came  from  Ginzburg.  °  Other  early  contributions 

1  q  20  p i 

came  from  Lashinsky  ’  and  Pierce.  Bolotovskii  has 
considered  the  spectrum  of  incoherent  Cerenkov  radiation  in 

pp 

guiding  structures. 

The  availability  of  relativistic  electron  beams  with 

respectable  current  densities  made  possible  the  recent 

experimental  effort.  At  Columbia  a  .5  MeV,  10  kA  beam  has 

been  used  to  drive  a  lucite  lined  waveguide. ^  At  Dartmouth  a 

100-250  keV,  1-40  amp,  pulsed  beam  has  driven  a  quartz  lined 

waveguide.  The  dependence  of  output  frequency  on  beam  energy 

was  measured  in  the  quartz  experiment  and  was  found  to  be 

12  24 

roughly  as  predicted.  ’ 

Further  experiments  at  Dartmouth,  presented  here,  have 
confirmed  the  predicted  frequency  dependence  on  liner 
dimensions,  dielectric  constant,  waveguide  mode  selection,  and 
beam  energy.  Power  measurements  have  been  made  as  well.^ 


These  results  are  presented  in  this  paper  along  with 
supporting  calculations. 

In  Chapter  II  the  electromagnetic  field  components  which 
couple  to  the  beam,  and  the  dispersion  relations  governing 
their  propagation,  are  discussed.  Three  possible  geometries, 
including  the  experimental  geometry,  are  considered  in  order 
to  illustrate  their  similarities  and  differences.  The 
predicted  coupling  frequencies  are  compared  to  measurements. 

In  Chapter  III  the  results  of  Chapter  II  are  used  to 
write  expressions  for  the  growth  rate  of  the  Cerenkov 
instability  in  both  the  collective  and  single  particle 
regimes.  The  criteria  for  choosing  a  regime  as  well  as  for 
overall  validity  of  either  analysis  are  explained  and 
evaluated  for  the  experimental  parameters. 

In  Chapter  IV  the  role  of  feedback  is  discussed.  An 
estimate  of  the  possible  output  power  is  compared  to  measure¬ 


ments  . 


CHAPTER  II 


2.1  FREQUENCY  DETERMINATION 

Several  considerations  restrict  the  set  of  frequencies  to 
which  the  beam  may  couple  in  a  dielectric-lined  waveguide. 

Both  beam  and  guide  are  assumed  to  be  uniform  in  the  longi¬ 
tudinal  ( z)  direction  and  in  one  transverse  (y  or  e) 
direction.  The  electric  and  magnetic  fields  must  satisfy 
Maxwell's  equations  and  the  boundary  conditions  imposed  by  the 
transverse  geometry.  These  define  all  the  possible  modes  of 
the  passive  guide. 

The  longitudinal  component  of  electric  field,  E  ,  must  be 
as  large  as  possible  in  the  volume  of  space  occupied  by  the 
beam  in  order  to  drive  the  axial  bunching.  Thus  TE  modes  are 
eliminated.  Uniformity  in  one  transverse  direction  assigns 
the  value  of  zero  to  the  corresponding  mode  index.  Thus  only 
TMon  modes  need  be  considered. 

Coupling  is  so  far  restricted  to  a  family  of  continuous 
spectra  as  shown  in  Fig.  2.1.  Angular  frequency,  w,  is  a 
function  of  wavenumber,  k,  but  k  is  a  free  parameter. 

Requiring  the  phase  velocity,  w/k,  of  the  radiation  to  be 
equal  to  the  beam  velocity  fixes  the  frequency  and  wavenumber 
for  a  given  TMgn  mode.  The  resulting  discrete  frequencies  are 


WAVENUMBER  [INVERSE  CM) 


Fie;.  2.1.  TM0  dispersion  curves  for  a  12.5  millimeter 
diameter  wave^Oide  with  a  two  millimeter  thick  annular  liner 
of  boron  nitride  (dielectric  constants  .  2)  .  The  beam  line 
(dotted  line)  must  lie  between  the  vacuum  and  dielectric  liftht 
lines  (so lid  lines). 


indicated  by  the  points  of  intersection  of  the  beam  line  with 

the  dispersion  curves  in  Fig.  2.1.  Phase  velocity  synchronism 

allows  the  beam  to  interact  with  Ez  ,  rather  than  having  its 

effect  averaged  to  zero  by  rapid  fluctuations  in  the  beam's 

rest  frame.  The  synchronous  condition  can  be  derived  by 

conserving  energy  and  momentum  for  forward  emission  of  a 

photon  by  an  electron  which  loses  a  small  fraction  of  its 
1  3 

momentum.  J 

2.2  DISPERSION  RELATIONS 

Dispersion  relations  for  uniform  dielectric-lined  guides 

have  been  derived  elsewhere  for  the  single  slab,  double 

slab  and  cylindrical  ’  D  geometries.  A  short  analysis  is 

given  here  which  accounts  for  the  three  geometries 

simultaneously  and  serves  to  illustrate  their  similarities  and 

differences.  The  geometries  are  shown  in  Fig.  2.2. 

Incorporating  immediately  the  choice  of  modes  (TMQn)  »  two 

of  the  boundary  conditions  (Ez  symmetric  at  the  origin  and 

vanishing  at  x=b),  and  the  harmonic  temporal  and  longitudinal 

i.  (  wt — k  z } 

spatial  dependence  of  the  fields,  e  ,  the  solution  to 

Maxwell’s  equations  for  the  x  or  r  dependence  of  the  electric 
fields  is  given  by 


i 
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Fig.  2.2a. 

Single  slab  geometry. 


ELLCTKON  BEAM 


a  <  x  <  b 


E 


z 


E 

x 


1(A)sn(qx  ,qb) 

( B) Cs( ipx ,0)  0  <  x  <  a 

( ik/q) (A)cs( qx  ,qb)  a  <  x  <  b 
(-k/p) ( B)Sn( ipx ,0)  0  <  x  <  a. 


(2.1a) 


(2.1b) 


where  x,  for  the  single  and  double  slab  geometries  and  r,  for 

the  cylindrical  geometry,  are  represented  by  x  in  these 

general  expressions.  The  perpendicular  wavenumbers,  q  and  p, 
2  o  ?  2  222? 

are  given  by  q  =ew  /c  -k  and  p  =k  -w  /c  ,  where  e  is  the 
dielectric  constant  of  the  liner.  Wave  numbers,  p  and  q,  are 
real  for  the  slow  wave  region,  c//T<w/k<c,  within  which  the 
beam  and  wave  may  interact. 

The  generalized  trigonometric  functions  sn ,  cs,  and  Sn 
and  Cs  are  to  be  identified  as  follows.  For  the  cylinder  sn 
and  cs  are  the  small  radial  sine  and  cosine  functions. 


sn(u,v)  =  (Jq(v)Y0(u)  -  Yq(v)Jq(u)  )/(2/irv)  (2.2a) 

c  s  (  u  ,  v )  =  (Y0(v)J1(u)  -  J0(v)Y1  (u)  )/(2/irv)  .  (2.2b) 


Sn  and  Cs  are  the  large  radial  and  cosine  functions  given  by 

Sn  (  u  ,  v )  =  (J1(v)Y1(u)  -  Y1  (v)  J1  (u)  )/(2/irv)  (2.3a) 

Cs(u,v)  =  (J1(v)Y0(u)  -  Y1(v)J0(u) ) /(2/wv)  .  (2.3b) 

Useful  properties  of  these  functions  are  given  in  Appendix  1. 
For  the  double  slab  and  single  slab  geometries, 


10 


sn(u,v)  =  Sn(u,v)  =  sin(v-u) 


(2.4a) 
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cs(u,v)  =  Cs(u,v)  =  -cos(v-u)  (2.4b) 

Incorporating  the  other  two  boundary  conditions  (E  and 
e(x)Ex  continuous  at  x=a)  gives 

(A)sn(qa,qb)  =  (B)Cs(ipa,0)  (2.5a) 

e ( Ai/q) cs( qa ,qb)  =  - ( B/ p) Sn( ipa , 0 ) .  (2.5b) 

The  dispersion  relation  is  the  quotient  of  the  above  two 
equations . 

( q) tn( qa ,qb) (Tn( ipa ,0)/i)  =  tp,  (2.6) 

where  tn=-sn/cs  and  Tn=Sn/Cs.  By  substituting  function 
definitions  appropriate  to  each  geometry,  the  familiar  forms 
of  the  dispersion  relations  and  field  expressions  can  be 
recovered . 

In  general,  it  is  suggested  that  any  dielectric-lined 
guide,  uniform  in  two  orthonormal  coordinates,  will  satisfy 
equations  2.1  and  2.6.  The  functions  sn ,  cs,  Sn ,  and  Cs  will 
be  constructed  from  the  independent  solutions  of  the  Helmholtz 
equation,  in  the  appropriate  set  of  orthonormal  coordinates, 
and  the  derivatives  of  those  solutions.  Hereafter  the  names 
sn ,  cs,  Sn ,  and  Cs  will  normally  be  reserved  for  the  radial 
trigonometric  functions  (Equations  2.2  and  2.3). 

-11- 
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2.3  LIMITING  CASES 


The  slow  wave  requirement  can  be  expressed  as  cose>1//T, 
where  e  is  the  Cerenkov  angle  defined  by  cose=1/e/T.  The  sine 
of  the  angle  of  incidence  for  radiation  in  the  liner  impinging 
upon  the  dielectric-vacuum,  interface  is  equal  to  cose.  Thus 
the  slow  wave  requirement  is  equivalent  to  requiring  total 
internal  reflection  at  the  d iel ectr ic- vacuum  interface.  The 
fields  in  the  beam  region  must  be  evanescent,  as  can  be  seen 
by  evaluating  the  field  expressions.  If  the  fields  fall  off 
sharply  enough,  it  can  be  anticipated  that  the  single  slab 
expressions  will  describe  the  system.  The  fields  fall  off 
fastest  for  grazing  incidence,  where  cose=1  or  B*1//T,  that 
is,  near  the  Cerenkov  threshold.  Fig.  2.1  indicates  that  a 
beam  velocity  near  the  Cerenkov  threshold  limit  intersects  the 
dispersion  curve  at  large  w  and  that  w*ck//T. 

The  single  slab  dispersion  relation  is 


(q)tan(qd)  =  ep 


(2.7) 


where  d=b-a.  In  the  single  slab  geometry  the  liner  and 
conductor  are  located  infinitely  far  out  on  the  positive  x 
axis  so  that  Tn( ipa ,0)= tanhC pa)  =  1  . 

The  double  slab  dispersion  relation  is 


( q) tan( qd ) tanh( pa)  =  ep, 


(2.8) 


In  the  near  -  threshold  limit,  tanh(pa)“tanh(wa/e-1)=1 ,  giving 
( q) tan( qd ) = e p ,  the  single  slab  dispersion  relation.  The 
cylinder  dispersion  relation  is, 


( q) tn( qa ,qb) I 1 ( pa) /IQ(pa)  =  ep 


(2.9) 


where  I1  and  IQ  are  modified  Bessel  functions  of  the  first 
kind.  Their  ratio  approaches  unity  for  large  arguments. 
Using  ct=1/tn,  the  general  dispersion  relation  (Eqn.  2.6) 
becomes  in  this  limit 


ct(qa,qb)  =  q/ep. 


(2.10) 


The  right  hand  side  of  Eqn.  2.10  approaches  zero  at  threshold, 
implying  that  q  is  approaching  one  of  the  roots  of  ct(qa,qb). 
Taylor  expanding  gives 


(q-9n) (d/dq) (ct( qa ,qb) )  =  q/ep, 

q-  qn 


(2.11) 


q  =  qn  1  - 


sn(qna,qnb)  N 

ep( bSn( qna  ,qnb)  -  asn( qna  ,qnb)  ) 


(2.12) 


where  qn  is  the  nth  root  of  ct(qa,qb).  The  latter  expression 
is  used  to  evaluate  growth  rates  along  the  dispersion  curves 
near  threshold.  For  the  single  and  double  slab  Eqn.  2.12 


becomes 


3 


q  =  ( (2n+1  )*/2d)  (1  -  1/Upd)  )  . 


(2.13) 
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Being  near  threshold  is  not  a  sufficient  condition  for 
the  cylindrical  expressions  to  approach  the  rectangular  ex¬ 
pressions.  However,  in  the  limit  of  high  mode  number,  n>>1, 
or  thin  dielectric,  a,b>>b-a,  ct(qa,qb)  evaluated  for  the 
cylinder  approaches  cot(qd).  °  Coupling  to  the  beam  in  the 
region  where  these  expressions  are  valid  is  difficult  for  lack 
of  fields  in  the  vicinity  of  the  beam,  as  the  growth  rate 
results  will  confirm. 

As  the  beam  velocity  approaches  the  speed  of  light,  p 
approaches  zero,  q  approaches  infinity,  and  the  small  radial 
tangent,  tn(qa,qb),  approaches  tan(qd).  The  dispersion 
relation  becomes 


B 


single  slab:  tan(qd)  =  0 
double  slab:  tan(qd)  =  e/qa 
cylinder:  tan(qd)  =  2e/qa 


(2.14a) 
(2. 14b) 
(2.14c) 


with  q=w/e -1  /c .  The  possibility  of  exploiting  the  wider  mode 

spacing  of  the  single  slab  geometry,  as  part  of  a  scheme  for 

1  ? 

an  infrared  Cerenkov  device,  has  been  examined.  J 


2.4  FREQUENCY  MEASUREMENTS 
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Quartz  (£=3.78),  stycast  (e=5),  and  boron  nitride  (e=4.2) 
have  served  as  the  liner  material.  An  electron  beam  with  a 
minimum  diameter  of  three  millimeters  is  provided  by  the  beam 
focusing  and  collimating  system  shown  in  Fig.  2.3.  The 
diverging  beam  is  collected  on  the  copper  guide  downstream 
from  the  liner  and  guide  field.  The  radiation  is  guided 
through  a  transmitting  teflon  vacuum  window  and  a  circular  to 
rectangular  mode  transducer  into  the  detection  hardware. 

Frequency  measurements  are  made  with  a  mechanically 
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driven,  free  space,  Fabry-Perot  interferometer  (Fig.  2.4). 

The  drive  speed  of  about  ten  microns  per  second,  and  the 
repetition  rate  for  beam  generation  of  about  one  Hz,  insure 
that  more  than  one  pulse  occurs  within  the  20  micron  minimum 
linewidth  of  the  interferometer. 

For  almost  every  liner  used,  the  fundamental  output 
frequency  is  near  the  frequency  at  which  the  TMQ1  mode  is 
synchronous  with  the  beam.  Figures  2.5,  2.6,  and  2.7  show  the 
frequencies  obtained  in  12.5  millimeter  diameter  waveguide.  A 
large  systematic  offset  is  noted  for  the  three  millimeter 
quartz  and  three  millimeter  stycast  liners.  A  surface  charge 
build-up  on  those  liners  is  suspected  of  decelerating  the  beam 
to  velocities  synchronous  with  the  observed  output. 

The  charge  per  unit  length,  q/1,  along  the  inner  surface 
of  the  liner  is  limited  by  the  dielectric  strength  of  the 
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Fig.  2.3.  Experimental  cylindrical  liner  geometry. 
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Fig.  2.5.  The  frequency  of  a  coupling  mode  is  plotted  vs.  the 
synchronous  beam  voltage  to  give  a  "tuning  curve".  The 
frequency  measurements  made  with  quartz  liners  of  one,  two, 
and  three  millimeter  thicknesses  are  displayed  along  with  the 
TMQ  tuning  curves  for  those  liners.  The  guide  diameter,  2b, 
is  i2.5  millimeters. 
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Fig.  2.6.  Data  and  TMQ.  tuning  curves  for  two  and  three 
millimeter  thick  boron  nitr ide  liners  in  12.5  millimeter 
guide . 
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liner.  Thus  qm./  1-VD/ (2  )  ln(  b/a)  ,  where  VD  is  the  breakdown 
voltage  given  by  the  dielectric  strength  times  the  thickness 
of  the  liner.  The  length,  1,  of  the  charged  surface  is  at 
least  b-a  due  to  surface  breakdown.  In  approaching  the 
charged  section,  an  electron  on  the  z  axis  gains  a  potential 
of  q/a.  For  three  millimeter  liners  b=2a,  and 
^max^^max^ 1  ^  ln(2  )  •  For  quartz,  stycast,  and  boron 

nitride,  V0  is  48,  36,  and  60  kV,  respectively,  and  the  beam 
can  lose  a  large  fraction  of  its  kinetic  energy. 

Quartz  liners  are  frequently  found  to  have  suffered 
damage  attributable  to  dielectric  breakdown.  The  use  of  a 
thin  conducting  film  to  bleed  off  charge  has  been  studied, 
and  other  less  drastic  methods,  such  as  a  split  liner,  may 
prove  helpful.  So  far,  no  coated  liner  has  produced 
rad iation . 

If  the  charge  is  distributed  over  a  longer  section  of 
liner,  greater  biasing  is  possible.  If  the  dielectric  surface 
has  little  tendency  to  hold  charge,  then  q  may  remain  much 
less  than  qmax.  For  two  millimeter  liners  b*3a/2,  and 

^max^^max^^^^B^1*^11^  ’  which  is  still  85%  of  the 

three  millimeter  value.  However,  the  average  gap  between  the 
beam  and  the  liner  is  greater  in  the  two  millimeter  geometry, 
and  the  biasing  is  less  noticeable. 

The  longitudinal  electric  field  associated  with  the  one 
millimeter  liner  geometry  is  too  weak  to  give  good  coupling 
for  moderate  beams.  At  the  higher  voltages  and  currents  which 
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become  necessary,  the  beam  generator's  performance  is  less 
consistent.  Thus  the  probability  of  recording  a  noise-free 
inter ferogram ,  requiring  approximately  1000  pulses,  is  low. 

Since  the  phase  velocity  of  any  of  the  waveguide  modes 
decreases  with  increasing  frequency,  a  favorable  scaling  law 
is  obtained  for  output  frequency  with  respect  to  beam  energy 
for  a  given  dielectric  liner.  The  range  of  frequencies 
obtainable  on  a  single  mode  has  reached  nearly  an  octave 
(Fig.  2.5). 

In  some  of  the  12.5  millimeter  diameter  waveguide 
experiments  and  in  a  9.5  millimeter  waveguide  experiment,  it 
has  been  possible  to  suppress  the  TMg^  mode  and  observe 
coupling  to  the  TMg^  mode  (Figures  2.7  and  2.8).  This  can 
occur  when  the  frequency  is  below  the  cutoff  for  the  TMQ1  mode 
of  the  empty  guide.  Frequencies  between  100  and  120  GHz  have 
been  attained.  The  power  appears  to  be  less  than  that 
obtainable  on  the  TMQ.j  mode,  but  still  significant.  Output  of 
still  lower  apparent  power  has  been  observed  through  a  200  Ghz 
filter  with  the  two  millimeter  quartz  liner  in  a  12.5 
millimeter  guide. 

All  the  shorter  wavelength  power  measurements  are  made 
with  1 N 5 3  diodes  in  a  Ka  band  mount.  This  system  is  difficult 
to  calibrate  when  it  is  far  out  of  band.  The  voltage  response 
generated  at  these  short  wavelengths  would  have  corresponded 
to  one  tenth  to  one  quarter  of  the  power  measured  in  the  8  mm 
region,  if  the  wavelengths  had  been  in  the  design  band  of  the 
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mount.  Hence,  at  the  shorter  wavelengths  the  output  power 
would,  at  a  minimum,  Jbe  in  the  one  kW  range. 
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CHAPTER  III 


3.1  POYNTING  FLUX,  STORED  ENERGY,  AND  THE  DISPERSION  FUNCTION 

The  dispersion  relation  can  be  expressed  as  D(w,k)=0, 
where  D(w,k)  is  identified  as  the  dispersion  function.  Since 
the  value  of  D  is  zero  everywhere  along  a  dispersion  curve, 
the  total  derivative  of  D  with  respect  to  a  displacement  along 
the  curve  is  also  zero.  In  terms  of  partial  derivatives, 

0  =  dD/dK  =  ( aD/aw) (dw/dk)  +  aD/ak.  (3.1) 

The  group  velocity,  dw/dk  can  then  be  expressed  as 

dw/dk  =  ( aD/ak)/( aD/aw) ,  (3.2) 

the  ratio  of  the  partial  derivatives  of  D. 

The  group  velocity  is  also  the  ratio  of  the  Poynting  flux 
to  the  stored  energy  per  unit  length,  giving 

dw/dk  =  P /( U/L )  =  ( aD/ak) /( 3D/aw) ,  (3-3) 

where  the  time  averaged  P  and  U  are  given  by 
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(3.4a) 

(3.4b) 


where  L  is  the  effective  length  of  the  region  of  interaction 
between  the  beam  and  lined  guide,  and  the  integrations  are 
carried  out  over  the  cross  sectional  area  of  the  guide. 
Rearranging  3-3  gives 


P/(3D/3k)  =  (U/L)/(3D/3w) , 


(3.5) 


a  ratio  which  can  still  be  a  function  of  w  and  k. 

Assuming  unity  magnetic  permeability  and  TMQn  mode 
symmetry  we  have 


Hz  =  Bz  =  Hx  =  Bx  =  0  (3-6) 

(w/ck)E  0  <  x  <  a 

Hv  =  B  =  (3.7) 

y  y  (we/ck)Ex  a  <  x  <  b 

where  x  and  y  are  the  general  transverse  coordinates,  perpen¬ 
dicular  and  parallel  to  the  liner  respectively.  The  P  and  U 
integrals  can  now  be  performed.  The  ratio, 

P/( 3D/ak) = (U/L )/( 3D/3w) ,  is  found  to  be 
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single  slab:  B2  ( w/ 1  6p)  (  A  y/ c  it  q)  ( 1  /si  nqd )  (3.8a) 

2 

double  slab:  B  (  w/8p)  (  Ay/ e  ir  q)  cosh(  pa)  sinh(  pa) /sinqd  (3.8b) 
2 

cylinder:  B  ( w/8p)  (  a/q)  IQ(  pa)  1 1  ( pa) /sn(  qa  ,qb),  (3.8c) 

where  B  is  the  electric  field  coefficient  defined  by  2.1  and 
Ay  is  the  spatial  extent  of  the  interaction  in  the  $ 
direction.  The  derivatives  are 

double  slab: 

3  D/3k=  kS{b-a( 1 -e )cos2-a( 1-q2/Ep2)sin2+(1/q)(1-q2/p2)sincos} 

(3.9a) 

3D/3w=(w/c2)S{eb-a(e-q2/ep2)sin2+(l/q)(e-q2/p2)sincos}  (3.9b) 

single  slab: 

3D/3k=kS{d+(1/q)(1-q2/p2)sincus)  (3.9c) 

3D/3w=(w/c2)S{ed+(1/q)(E-q2/p2)sincos}  (3.9d) 

sin=sin(q(b-a))=sinqd  cos=cos(  q(  b-a)  )  =  cosqd 

cyl inder : 

3D/3k=kSlb2/a-a(1-e)cs2-a(1+q2/ep2) sn2- (2/q) ( 1  +q2/p2 ) sncs) 

( 3 . 9e) 

3D/3w=  (w/c2)S{  cb2/a-a(  £+q2/e p2)sn2-(2/q)  (  E  +  q/p2 )  sncs)  (3 .90 
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sn=sn(qa,qb)  cs= cs( qa  ,  qb)  , 

where  S=e  p(  Ct(  ipa  ,0) /i) /( (  q)  sn(  qa  ,qb)  )  in  each  case,  and  a 
convenient  form  is  used  for  D, 

D(w,k)  =  - ( q) sn( qa  ,qb)  -  ( ep) cs( qa ,qb) ( i( Ct( ipa  ,0) )  .  (3.10) 

3.2  COLLECTIVE  VS.  SINGLE  ELECTRON  GROWTH  RATE 

The  results  of  the  previous  section  are  used  in  Section 

3.3  for  the  calculation  of  growth  rates.  In  this  section  beam 
space  charge  interactions  and  beam  energy  spread  are  discussed 
with  respect  to  the  growth  rate.  It  is  concluded  that  the 
present  experiment  may  be  characterized  by  a  cold,  collective 
growth  rate. 

Assuming  the  Cerenkov  radiation  due  to  the  beam  and  the 
axial  bunching  of  the  beam  interact  as  a  linear  feedback 
system,  the  radiation  field  will  grow  exponentially.  The 
growth  rate  can  be  recovered  from  the  time  averaged  stored 
field  energy,  U,  given  by  Eqn.  3.4c.  Then 

w''  =  ( dU/dt) /2U  ,  (3.11) 
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where  w' ’ ,  the  growth  rate,  is  the  imaginary  part  of  the 
frequency,  w=Wg+w' + iw' ' .  The  no-beam  solutions  are  now 
designated  wQ,  kg  where  D(Wg,kg)=0. 

The  growth  rate  is  actually  temporal  only  in  the  electron 
beam’s  rest  frame.  The  spatial  gain  in  the  lab  frame  is  given 
by  k''=w''/v.  In  the  beam  frame  the  Poynting  flux  is  uniform 
in  space,  resulting  in  no  net  radiative  power  flow  out  of  a 
longitudinal  section  of  guide.  The  increase  in  stored  energy, 
dU/dt,  is  then  due  entirely  to  the  energy  given  up  by  the 
electrons  to  the  fields,  -  ( 1  /2  )  Re  JdV  ( 3*ii )  .  The  modulated 
portion  of  the  current  density,  the  only  portion  which 
contributes  to  this  integral,  is  related  to  the  electric  field 
through  Vlasov's  equation  for  the  electron  phase  space 
distribution  function. 

In  the  single  electron  case  the  direct  interaction  of 

electrons  with  one  another  is  ignored.  The  spatial  modulation 

of  the  electron  distribution  enters  the  problem  only  through 

Vlasov's  equation.  A  positive  growth  rate  is  possible  for  an 

electromagnetic  wave  whose  phase  velocity  is  slightly  less 

than  the  beam  velocity.  A  cold  beam  analysis  has  been 
29 

performed.  The  dependence  of  the  growth  rate  on  the 
velocity  difference  for  a  cold  beam  is  a  well  known  result  and 
is  independent  of  the  wave  supporting  structure.  The 
dependence  of  the  growth  rate  on  the  detuning  parameter, 

( w-kv) (L/2v) ,  is  shown  in  Fig.  3.1. 
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Electrons  may  influence  one  another  through  the  space 

charge  interaction  (Poisson's  equation)  if  given  sufficient 

time.  Space  charge  oscillations  can  serve  to  enhance 

bunching.  The  transit  time  of  electrons  through  the 

interaction  region,  L/v,  must  exceed  the  time  required  for  a 

V2  2  2 

plasma  oscillation,  2*y  /w  ,  where  w  =4*ne  /m.  In  this 

P  P 

case  the  growth  rate  can  exceed  that  which  is  calculated  from 
Vlasov's  equation  alone. 

A  direct  calculation  of  dU/dt  is  now  less  convenient. 

Instead,  the  beam  can  be  treated  as  a  dielectric  medium,  and 

the  beam  dispersion  function  becomes  coupled  to  the  passive 

guide  dispersion  function.  A  new  dispersion  function  is 

calculated,  and  its  roots  can  have  imaginary  parts  for  phase 

velocities  less  than  the  beam  velocity. 

The  procedure  for  calculating  the  dispersion  function  is 

analogous  to  the  usual  derivation  of  electrostatic  plasma 
12  30 

oscillations.  ’  In  the  cold  beam  limit  the  Vlasov  equation 
reduces  to  the  simple  continuity  equation  (Eqn.  4.3).  The 
electron  density  and  velocity  perturbations  are  assumed  small, 
permitting  linearization.  The  guide  boundary  conditions  must 

still  be  satisfied.  The  dispersion  functions  for  the  single 
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slab,  double  slab,  and  cylinder  have  been  calculated. 

The  result  is  given  as 

D(w,k,wp2)  =  (q)tn(qa,qb)(Tn(  ipa/T^  , 0 ) / i )  -  p //F£,  (3.12) 


-  31 


where 


eb  =  "I  -  wp2/U3(w-kv)2) .  (3.13) 

In  the  absence  of  a  beam,  Eqn .  3.12  clearly  gives  the 
earlier  result  (Eqn.  2.6)  for  D(w,k)=D(w,k ,0) .  The  roots  of 
D(w,k,Wp)  are  shifted  from  the  roots  of  D(w,k,0)  along  the 
real  and  imaginary  axes. 

In  the  absence  of  a  dielectric,  Eqn.  3.12  for  D=0  reduces 
to 

(P'//Tpct(p'b/I^,0)  =  0,  (3.14) 

2  2  2  ? 

where  e=1  and  a=  b .  The  definition  p'  =w  /c-k  is  used  in 
anticipation  of  fast  waves.  This  gives  an  empty  guide 
dispersion  function, 

D(w,k,wp2)  =  (1  -  wp2/(Y3(w-kv)2))(w2/c2-k2-xn2/b2)  ,  (3.15) 

where  xn  is  the  nth  root  of  Ct(x,0), 

cylinder:  JQ(xn)  =  0  (3.16a) 

double  and  single  slabs:  cos(xn)  =  0.  (3.16b) 

The  roots  of  Eqn.  3.15  are  all  real,  implying  no  growth. 
For  the  single  slab,  b  is  taken  to  be  infinitely  large.  The 
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usual  dispersion  function  for  electrostatic  plasma 

2  "3  2 

oscillations,  1-wp  /(y  (w-kv)  ),  results,  implying  that  a  half 
space  bounded  by  an  infinite  conducting  sheet  is  not  a 
waveguide . 

The  no-guide  dispersion  relation  is  satisfied  by  two 
modes , 

w  =  kv  ±  wp2/y3/2 ,  (3-17) 

the  fast  (+)  and  slow  (-)  space  charge  modes.  Modes  of  this 
nature  exist  in  the  presence  of  the  guide  as  well.  Only  the 
slow  space  charge  (negative  energy)  mode  can  couple  energy 
from  the  beam  into  the  radiation  field.  The  frequency  spread 
due  to  the  beam's  velocity  spread,  kAv,  must  be  less  than 

2  7.  /2 

2 Wp  /y  in  order  to  resolve  the  slow  space  charge  mode  from 

the  fast  space  charge  mode.  This  criterion  can  also  be 
expressed  as  an  upper  limit  on  the  beam's  fractional  energy 
spread  , 

Ay/ y  <  (By)2(?wp2/wy3/2).  (3.18) 

For  the  typical  experimental  parameters,  y=1.2,  wp=1010  sec-1, 
and  w=2x1 01 1  sec”1 , 

Ay/y  <  .05  (3.19) 
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is  necessary.  The  real  part  of  the  beam- induced  frequency 

shift,  or  detuning,  is  negative  for  the  growing  solution  to 
2 

D(w,k,Wp  )  =  0 .  Thus  it  can  help  to  resolve  the  space  charge 

"3  1 

modes.  This  shift  is  of  the  order  of  the  growth  rate  and,  as 
will  be  seen,  can  be  a  significant  fraction  of  Wp .  A  negative 
shift  also  means  that  the  phase  velocity  of  the  wave  is 
slightly  less  than  the  beam  velocity. 

In  Section  3.3  the  single  electron  growth  rate  is 
compared  to  the  collective  growth  rate.  For  this  reason  the 
restriction  on  beam  energy  spread  pertaining  to  the  single 
electron  case  is  stated  below.  For  significant  growth  the 
detuning  parameter  should  have  the  value  which  gives  the  large 
maximum  in  Fig.  3.1.  Expressed  as  a  frequency,  the  optimum 
detuning  found  for  the  single  electron  analysis  is  -»v/L. 

This  detuning  must  be  resolved  with  respect  to  kAv  .  Thus 

A  y/y  <  (By>2Uv/wL)  (3-20a) 

or 

Ay/ y  <  (By)2(x/2L)  (3.20b) 

or 

Ay/y  <  .02  (3.20c) 

for  the  experimental  parameters  above. 

Energy  spread  is  unavoidable  when  producing  collimated 

_3 

charge  beams,  although  Ay/y=10  is  possible  with  existing 
technology.  A  particle  simulation  code  is  available  to 
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calculate  the  trajectories  of  charged  particles  in  the 

vicinity  of  potential  surfaces  and  magnetic  field  which  have 
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rectangular  or  cylindrical  symmetry  (see  App.  III).  The 
charges  are  moved  a  small  distance  in  the  direction  of 
£+(v/c)x§.  Then,  the  electric  potential  field  is  adjusted  to 
satisfy  Poisson's  equation  at  all  the  mesh  points  of  the 
simulation  for  the  new  charge  distribution.  Self  magnetic 
fields  are  added  to  the  external  field. 

The  code  was  run  for  a  simplified  model  of  the  beam 
focusing  system.  Mesh  points  were  .5  millimeters  apart.  The 
maximum  transverse  energy  is  acquired  by  electrons  at  the 
maximum  beam  radius  (Figures  3.2  and  3*3).  It  is 

Ay/y  =  .023.  (3.21) 

The  average  transverse  energy,  weighted  by  the  beam  density 
profile  given  by  the  simulation  (Fig.  3.10,  is 

&y/y  =  . 01 .  (3.22) 


Thus  the  cold  beam  criteria  for  the  collective  and  single 
electron  regimes  are  just  satisfied  by  the  simulated  beam. 

The  transit  time  of  the  beam  through  the  resonator  lies 
in  the  range 


0.6  nsec  <  L/v  <  1.0  nsec. 


(3.23) 
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FIELD  LINE 


cHu<ui,  oyaucu  uuuiai  r  aan  uar  ieitj  .  me  outermost  trajectories 
are  least  well  aligned  with  the  guiding  magnetic  field.  The  coordinates  are  given 
in  mesh  units.  A  mesh  unit  is  defined  to  be  .5  millimeters  for  this  run  of  the 
particle  trajectory  code. 


o 


.  1 - 

.0 

FEl 


T~ 

1.3 


GUN 


i  r 

2.0  2.  C  4  .0 

26  JRNUART  1 S82 


6-0  ?.'o 


Fig.  3.4.  Final  simulation  beam  current  density  as  a  function 
of  radius,  normalized  to  the  maximum  value. 
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The  time  required  for  a  plasma  oscillation  lies  in  the  range 

K 

0.2nsec  <  2irY^/^/wp  <  1.7  nsec.  (3.24) 

^  There  is  hope,  then,  that  the  space  charge  interaction  may 

contribute  to  the  growth  rate. 


3.3  GROWTH  RATES  FOR  THE  CYLINDER 

To  evaluate  the  collective  growth  rate,  the  complex 

solution  to  the  combined  dispersion  relation,  Eqn .  3.12,  must 

be  found.  At  present,  codes  are  available  which  find  the  real 

solutions  to  the  no-beam  dispersion  relation,  Eqn.  2.6, 

(Appendix  III).  They  are  not  capable  of  solving  Eqn.  3.12, 

but  could  possibly  be  adapted  to  do  so.  Equation  3.12  has 

been  solved  in  the  weak  beam  approximation  for  the  single 
31  12 

slabJ  and  double  slab.  In  this  approximation  the  detuning 
and  growth  rate  due  to  the  beam  are  assumed  to  be  small  with 

p 

respect  to  wQ.  Thus  D(w,k,Wp  )  may  be  Taylor  expanded  about 
wQ,  kg  with  respect  to  w  and  with  respect  to  beam  density, 


D(w,k,wp2)  =  DQ  +  ( w-Wg ) 3Dq/ 3 w  +  (Wp2)3DQ/3( wp2) ,  (3.25) 


D0  =  D^wo»kO'0^  * 
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The  dispersion  relation  becomes 


w—w 0  =  wp2(3D(w0,k0l0)/3(wp2))/OD(w0,k0,0)/aw)  .  (3.26) 

For  the  cylinder 

3  D/ 3 ( wp2 )  =  -  ( 1 /2 y3 ( w-kv) 2 )(eplg(pa)/qsn)( q2asn2/e-e  p2acs2 ) , 

(3.27) 


using  D( w ,k , wp  )= q/F^l ^ ( p/e^a) sn- e pig ( p/e^a) cs .  It  is  further 
assumed  that  the  roots  are  near  synchronism,  w«kv ,  k*kg. 
Already  knowing  3D/3w  (Eqn.  3.9f),  we  have 


(W-Wg)3= 


2  2  2_  2 ,  _ 2  2 
-w  c  a  q  sn  /c  -  ep  cs 

2Y3w  ( cb2/a)-a( e-q2/ep2)sn2-(2/q) (e-q2/p2 


q  s  w/c/c-1/32,  p  =  w/cBy. 


- (3.28) 

)  sncs 


Equation  3.28  has  three  roots  for  w-wQ:  w' ,  w’+iw’1,  and 
w'-iw' The  detuning,  w'  ,  and  the  growth  rate,  w' *  ,  are 
given  by 


w'  =  |(w-w0)3l1/3  (3.29a) 

w”  =  (/T/2)w  (3.29b) 

The  no-beam  dispersion  curves  of  Fig.  2,1  are  replotted 
in  Fig.  3.5  for  the  lowest  TMgn  modes,  along  with  w' '  (k).  •  The 
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FREQUENCY  (INVERSE  CM) 


0  10  20  30  40 


WAVENUMBER  (INVERSE  CM) 


Fig.  3.5.  The  two  lowest  TM_  no-beam  dispersion  curves  and 
associated  growth  rates.  Attach  wavenumber  the  beam  is 
chosen  to  have  a  velocity  synchronous  with  the  no-beam  mode 
and  an  electron  density  of  10  cm”  .  The  liner  is  two 
millimeter  thick  boron  nitride  in  12.5  millimeter  guide. 


-  in  - 


4 


synchronous  beam  velocity  is  assumed  at  each  wavenumber.  A 

significant  growth  rate  for  the  TM^  mode  is  indicated  and  has 

led  to  output.  The  frequencies  corresponding  to  TMQ2  and 

higher  mode  coupling  (Fig.  3.6)  lead  to  violation  of  the  cold 

beam  criteria  (Equations  3-18  and  3.20).  Including  warm  beam 

13  31 

effects  will  reduce  the  calculated  growth  rates.  ’ 

The  peak  gain  as  a  function  of  velocity  is  plotted  in 
Fig.  3-7  for  several  dielectric  constants.  Liner  dimensions 
and  electron  density  are  chosen  for  comparison  to  a  previous 
double  slab  calculation  (Fig.  3.8).  The  cylinder  exhibits  a 
sharper  turn-on  at  threshold  and  a  slower  fall-off  than  the 
double  slab.  The  optimum  dielectric  constant  drops  from  about 
four  to  about  three. 

The  double  slab  growth  rate  expression  is  given  for 
comparison  as  well. 

2  2  2  2  2  2 
_  -w  c  a((q  /E)sin  -Ep  cos  )  +  qsincos 

(w-w0)^=  — ^ - 2 - 2 - ^ - - ,(3.30) 

2y  w  (eb-a(e+q  /ep  )sin  +1/q(e+q  /p  )sincos 

where  sin=sin(qd)  and  cos=cos(qd).  Conventional  trigonometric 
functions  replace  the  radial  trigonometric  functions.  The 
cross  term  in  the  numerator  no  longer  vanishes,  as  in  the 
cylinder  expression,  and  the  cross  term  in  the  denominator 
changes  by  a  factor  of  -2.  The  eb^/a  term  in  the  denominator 
of  Eqn .  3.28  is  now  ea  as  a  result  of  the  difference  in 
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Fig.  3*6.  Higher  mode  collective  growth  rates  for  the 
parameters  of  the  previous  figure. 
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Iw-kv/l  (INVERSE  SEC) 


Fig.  3.7.  Growth  rate  as  a  function  of  velocity  for  the  TMQ1 
mode  of  a  2.5  millimeter  thick  liner  in  ten  millimeter  '  1 
diameter  guide;  nslO1  cm"  .  The  dielectric  constant  is 
allowed  to  vary.  A  value  of  three  is  about  optimum. 


normalization  between  the  radial  and  rectangular  trigonometric 
functions  (Eqn.  A1.4). 

Taking  the  single  slab  limit  of  Eqn.  3.30  gives 


3  wn2°2 

(U-V  W- 


2  2  2  2 
a(ep  c  o  3  -(q  /£)sin  )-qsincos 


2yJw  e d+(a/p^)(ep^cos^-(q^/e)sin^)-(1/q) (e+q2/p2)sincos 


(3.31 ) 


The  single  slab  dispersion  relation  gives 

2  2  2  2 
ep  cos  -(q  /£)sin  =0.  Thus 


2  2 

■3  W  C 

(  U.U  )  *  —  - £—  .. 

w0;  -  .  3 


p  p  p 

q^  w  w  1  eB  -1 

_  _  _j> _ # _ _ _ 

2yJw  sqd+E+q^/p^  2y^  £qd+Y^(£-1) 


(3.32) 


The  fall-off  of  the  growth  rate  at  threshold  and  at  high 
energy,  indicated  by  Figures  3*7  and  3.8  for  the  cylinder  and 
double  slab,  is  apparent  in  the  expression  for  the  single 
si  ab . 

Experimentally,  liner  thickness  is  varied  as  well  as 
dielectric  constant.  A  thicker  liner  always  results  in  a 
higher  calculated  growth  rate,  assuming  the  beam  current  can 
be  maintained  for  the  smaller  hole  diameter.  Furthermore,  a 
smaller  hole  sees  less  of  the  curvatu.  e  present  in  the  guide 
field  now  in  use.  The  radial  magnetic  field  at  the  ends  of 
the  liner  increases  with  radius  and  with  distance  from  the 
mid-plane  of  the  guide  field.  This  guide  field  curvature 
leaves  a  gap  between  the  beam  and  the  dielectric,  degrading 
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the  growth  rate.  The  gap  is  largest  at  the  mid-plane  and  is 

largest  overall  for  a  thin  liner. 

In  order  to  choose  a  liner  thickness  for  optimum  growth 

rate  with  a  particular  beam  generating  system,  the  growth  rate 

should  somehow  account  for  beam-dielectric  gaps  and  for  non- 

uniform  beam  cross  sections  (Fig.  3.4).  In  principle,  one 

could  derive  a  new  dispersion  relation  for  a  beam,  gap,  and 

lined  guide  by  matching  the  new  boundary  conditions  for  Ez  and 
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Er  at  the  beam-vacuum  and  vacuum-dielectric  interfaces. 
However,  to  account  for  beam-dielectric  gaps  and  nonuniform 
beams  in  a  more  convenient  manner,  an  effective  beam  parameter 
can  be  used  . 

An  effective  beam  density, 


rrb  . 

n»  =  nJ0 Ez2(r)rdr  /Ez2(0)irb2, 


(3.33) 


for  a  solid  uniform  beam  has  already  been  incorporated  into 
the  calculation  of  a  single  electron  growth  rate  for  the 
cylinder.  For  the  single  electron  case,  n'  is  an  exact 
effective  beam  density.  The  reason  for  this  is  that,  in  the 
Vlasov  analysis  for  a  weak  beam,  electrons  lose  energy  at  a 
rate  proportional  to  the  square  of  the  local  longitudinal 
electric  field  amplitude. 

In  the  collective  case  n'  serves  as  an  approximate 
effective  beam  density.  An  alternate  form,  more  useful  for 
the  collective  case,  is 


-  L~- 


n’  =  n/oEz2(r>rd^OEz2(r)rdr' 


(3-34) 


as  n'=n  when  the  beam  just  fills  the  hole. 

Without  presenting  any  of  the  details  of  the  single 

electron  calculation,  the  growth  rate  will  be  evaluated  in  the 

cylindrical  geometry  for  comparison  to  the  collective  growth 

rate.  A  general  single  electron  expression  resulting  from  the 
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cold-beam  Vlasov  analysis  7  is  directly  quoted  as 


w  lb  n'  F3  /sir?  n* 


(By)  V„  n  ark  n 


(3.35) 


Some  definitions  are  required.  Ib  is  the  beam  current,  and 

is  the  Alfven  current,  ec/rQ  =  17  kA,  where  rQ  is  eVmc  ,  the 

classical  radius  of  an  electron.  is  the  mode  volume,  taken 

m 

2 

to  be  wb  L  for  the  cylinder.  n  is  the  dimensionless  single 
electron  detuning  parameter,  ( w-kv) (L/2v) ,  introduced  in 
Fig.  3.1.  The  maximum  value  of  the  n-dependent  factor,  known 
as  the  lineshape  function,  is  seen  in  Fig.  3.1  to  be  about 
.55.  F  and  n'  are  the  geometry-dependent  factors. 

Using  Eqn .  3.33, 


n’/n  =  (rb2/b2)(IQ2(prb)  -  I^Cpr^)). 


(3.36) 
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F  is  the  ratio  of  the  mode  energy  density  on  axis  due  to  Ez , 

p 

|  E  z  ( 0 )  |  /8",  to  the  average  total  energy  density,  U/Vm.  For 
the  cylinder 


(b2/a)  (q2c2/g  w2)(sn2/I  2(pa)  ) 

F  =  — p - P - p - p - 1 - p — p - ,  (3.37) 

eb  /  a- a(  c  +  q^/e  )  sn  - (2/q)  ( e+q^/p^ )  sncs 

where  again  sn=sn(qa,qb)  and  cs= cs( qa ,qb) . 

The  resulting  velocity  dependence  of  the  growth  rate  with 
the  lineshape  function  held  at  .55  is  shown  in  Fig  3*9  for  the 
parameters  of  Figures  3.7  and  3.8.  The  growth  rates  are  low 
but  not  insignificant.  The  optimum  dielectric  constant  is 
about  seven  in  this  case. 

In  Fig.  3.10  the  liner  thickness  is  varied  while  holding 
the  beam  radius  and  guide  diameter  constant.  The  growth  rates 
fall  off  sharply  as  the  gap  between  the  beam  and  the  liner 
widens.  Use  of  the  single  electron  effective  beam  density 
(Eqn.  3. 3^)  for  the  collective  calculation  gives  similar 
behavior  (Fig.  3-11).  In  both  cases,  the  peak  growth  rate 
shifts  to  higher  velocities  as  the  gap  widens  and  milder  Ez 
evanescence  becomes  more  important. 
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Fig.  3.10.  Single  electron  growth  rate  for  different  lin 
thicknesses.  The  beam  radius  is  held  at  three  millimeter 
The  guide  diameter  is  12.5  millimeters. 
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Fig.  3.11.  Collective  growth  rate  (wf,=  3/2iw-kv,)  for 
different  liner  thicknesses.  The  beam  radius  is  held  at  three 
millimeters.  The  guide  diameter  is  12.5  millimeters.  An 
effective  beam  density  is  used  to  account  for  the  gap  between 
the  beam  and  the  dielectric. 


CHAPTER  IV 


4.1  FEEDBACK 


In  order  for  the  present  device  to  produce  high  power,  it 
must  operate  as  a  feedback  oscillator.  The  guide  field 
geometry  restricts  the  interaction  region  to  a  maximum  length 
of  about  15  cm.  This  severely  limits  the  power  available 
through  single  pass  amplification. 

Experimental  evidence  for  multiple  passes  through  the 
liner  by  the  radiation  helps  to  explain  the  observation  of 
large  output  signals.  Most  noticeable  is  the  increased  beam 
energy  necessary  for  oscillation  when  the  impedance  mismatch 
at  either  end  of  the  liner  section  is  made  less  drastic. 

Also,  output  is  preferred  at  a  set  of  discreet  beam 
velocities,  coupling  to  the  appropriate  TMQ1  synchronous 
vrequencies.  This  set  of  discreet  frequencies  must  be  those 
for  which  the  multiply  reflected  waves  traveling  downstream 
add  constructively. 

In  a  multiple  pass,  oscillator  system  the  round  trip 
growth  must  exceed  the  round  trip  loss.  The  beam  is  assumed 
not  to  interact  with  waves  traveling  u-istream.  Thus 


w'  'T  >  -  InR .  -  InR  , 
d  u 
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(4.1) 


V 


where  T  is  the  one-way  transit  time,  L/v,  and  Ru  and  Rd  are 
the  fractions  of  incident  power  left  in  the  original  mode 
after  reflection  upstream  and  downstream  respectively.  All 
losses  due  to  absorption  are  assumed  to  be  negligible. 

Ideally,  Ru  should  be  unity  if  all  the  power  is  to  be  coupled 
out  downstream.  The  fields  are  concentrated  in  the  liner  for 
w/k<c,  allowing  a  reflective  washer  to  achieve  high  Ry  at  the 
upstream  end  of  the  liner. 

A  scheme  for  adiabatically  compressing  an  annular  beam  in 
the  region  between  an  on-axis  reflector  and  the  liner  has  been 
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proposed.  Such  a  scheme  presents  difficulties  with  respect 
to  beam  dynamics  and  millimeter  wave  optics,  but  the  added 
efficiency  may  be  worth  the  effort  for  some  applications. 

Reflections  from  several  downstream  elements  are 
possible.  Tapers,  bends,  and  horns  should  reflect  a 
negligible  amount;  the  vacuum  window  and  the  circular  to 
rectangular  mode  transducer  may  contribute  a  bit  more.  The 
reflections  from  the  downstream  end  of  the  liner  section  and 
the  Fabry-Perot  inter ferometer  are  the  most  significant  and 
help  to  explain  some  experimental  observations. 

The  liner  geometry  is  shown  in  Fig.  2.3.  An  upper  limit 
to  the  reflection  obtained  at  the  interface  between  the  lined 
and  empty  guides  will  be  that  obtained  at  an  infinite 
dielectric-vacuum  interface.  The  angle  of  incidence  is  non¬ 
zero  for  TM  modes.  In  terms  of  the  longitudinal  wavenumbers, 
k  and  k’ ,  for  the  dielectric  and  the  vacuum, 


Rd  =  ((k/c-k* )/(k/e+k* ))2.  (4.2) 

The  same  expression  is  obtained  for  the  junction  of  a 
completely  filled  guide  and  an  empty  guide  with  k  and  k' 
evaluated  for  the  TMq^  modes. ^ 

The  reflection  at  the  end  of  the  lined  guide  is  actually 
lower  for  two  reasons:  (1)  because  not  all  of  the  electric 
field  is  in  the  liner,  the  effective  dielectric  constant  in 
the  lined  guide  is  less  than  e;  and  (2),  the  radial  non¬ 
uniformity  of  the  lined  guide  will  cause  higher  TMQn  modes  to 
be  launched  and  reflected.  These  reduce  the  power  reflected 
in  the  TMQ1  mode  in  a  way  not  described  by  Eqn .  4.2. 

The  data  in  Fig.  4.1  indicate  a  dependence  on  end  reflec¬ 
tivity  for  the  threshold  beam  voltage.  The  tapered  end  is 
less  reflective  than  the  square  end,  and  a  check  can  now  be 
made  as  to  whether  or  not  feedback  from  the  end  of  the  square 
end  liner  is  sufficient  to  explain  the  output  at  the  lower 
voltage. 

The  dispersion  relations  for  the  lined  and  empty  guides 

give  k  and  k' .  For  the  loss  at  the  end  of  the  liner,  -InR^, 

the  value  9.5  dB  is  obtained  at  50  GHz  and  10.8  dB  at  30  GHz. 

The  additional  loss  obtained  by  tapering  the  end  would  be  at 
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least  10  dB  for  a  final  liner  thickness  of  zero.  However, 
the  final  thickness  achieved  was  still  about  half  the  original 
thickness.  Thus  only  a  few  dB  of  additional  loss  are 
guaranteed  . 
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The  growth  rate  for  the  lower  voltage,  50  GHz,  output  is 

significantly  lower  than  the  growth  rate  for  the  higher 

voltage,  30  GHz  output,  if  the  gap  between  the  beam  and  the 

liner  is  accounted  for  by  an  effective  beam  density.  The  gap 

is  about  one  millimeter.  The  middle  curve  in  Fig.  3.11  gives 

1 0  -3 

these  growth  rates  assuming  n=10  cm  J  (e.g.  a  7.5  ampere 
beam  of  uniform  cross  section  with  velocity  .55c  and  radius 
three  millimeters)  . 

The  square  end  liner  experiment  employed  a  20  ampere  beam 
(Fig.  4.2)  with  L=15  cm  giving  w’'T=8.4  dB  at  50  GHz  and 
w' 'T=13-2  dB  at  30  GHz.  Thus  the  oscillation  criterion, 
Inequality  4.1,  is  almost  satisfied  by  the  square  end  liner  at 
50  GHz  and  well  satisfied  at  30  GHz.  The  approximations  made 
in  calculating  w’ •  may  have  led  to  an  artificially  low  result. 
The  detuning  and  the  plasma  frequency  which  are  assumed  to  be 
small,  are  each  about  4J  of  Wq. 

If  the  50  GHz  oscillation  is  assumed  to  take  place  just 
above  threshold,  then  a  few  dB  of  additional  loss  due  to 
tapering  will  quench  it,  but  can  still  permit  oscillation  in 
the  higher  growth  rate  region  near  30  GHz,  as  the  data 
indicate.  Removal  of  the  upstream  copper  reflector  will  also 
increase  losses  and,  indeed,  was  found  to  increase  the  voltage 
threshold . 

The  interferometer  (Fig  2.4)  can  be  a  source  of  feedback. 
The  reflectivity  of  the  perforated  brass  mirrors  is  well  above 
90%  at  35  GHz.  The  largest  uncontrolled  loss  of  signal  is  by 


scattering  out  of  the  open  cavity  formed  by  the  launching 
horn,  collimating  lens,  and  upstream  mirror.  Feedback  is 
normally  reduced  further  by  a  waveguide  attenuator,  intended 
primarily  to  protect  the  diode  detector  of  the  interferometer. 
Even  so,  a  reflected  signal  only  10  dB  down  is  possible. 

When  the  interferometer  is  tuned  to  the  output,  both  the 
transmission  through  to  the  detector  and  the  reflection  back 
to  the  liner  are  maximized.  Amplified  feedback  should,  in 
this  configuration,  result  in  line-narrowing.  The  inter- 
ferograms  in  Fig.  4.3  have  narrower  than  normal  linewidths. 
They  also  have  signals  at  mirror  separations  tuned  to 
harmonics  of  the  output,  which  are  simply  higher  order  modes 
of  the  interferometer.  The  presence  of  harmonics  is  enhanced 
for  a  sharply  tuned  driving  signal.  This  effect  is  important 
to  the  operation  of  frequency-doubled  and  frequency-tripled 
lasers  and  musical  instruments. 

Once  the  threshold  oscillation  condition,  Inequality  4.1, 
is  satisfied,  whatever  the  source  of  feedback,  the  beam  pulse 
duration  of  100  to  1000  transit  times  guarantees  the  survival 
of  only  self-reproducing  cavity  modes.  At  a  given  point  in 
the  waveguide,  the  signal  will  have  the  same  phase  after  each 
round  trip.  The  or  TMqj  fields  can  satisfy  this 

condition  only  at  a  discreet  set  of  frequencies  or  longi¬ 
tudinal  modes. 

Observation  of  this  effect  is  simplified  by  the  slow  rise 
time  of  the  single  stage  pulse  forming  network  currently  in 
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23  Nov. '81  4C 
Pep  Rate  3.5 
Drive  Speed  20 


30  Nov  *81  IIC 
Rep  Rate  5  Hz 
Drive  Speed  18 
Gain  5  mV 


Fig.  4.3*  Inter ferograms  for  two  millimeter  thick  boron 
nitride  indicating  some  control  over  the  signal  by  the 
interferometer . 
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use.  Output  is  most  often  generated  when  the  beam  is  a  its 
peak  voltage  as  in  Fig.  4.4a.  A  small  increase  in  the  peak 
voltage,  which  is  insufficient  to  reach  synchronism  with  the 
adjacent  longitudinal  mode,  often  leads  to  output  on  the 
original  mode.  However,  the  output  is  now  generated  when  the 
beam  voltage  is  at  its  original  peak  value,  earlier  or  later 
in  the  pulse  than  the  new  peak,  or  both  earlier  and  later. 
(Figures  4.4b,  4.4c,  and  4.4d). 

Fig.  4.5  shows  the  signal  behavior  as  the  beam  voltage  is 
swept  from  oscillation  threshold  to  about  twice  that  value. 
Only  three  values  of  the  beam  voltage  correspond  to  output. 
When  such  a  value  is  realized  for  only  a  short  time,  as  on  the 
rise  of  the  beam  pulse,  the  duration  of  the  output  is  shorter. 
The  synchronous  frequencies  corresponding  to  these  voltages 
are  spaced  a  few  GHz  apart  as  confirmed  by  inter ferograms . 
Normally,  when  Inequality  4.1  is  satisfied,  the  single  pass 
gain  is  high  enough  to  generate  output,  masking  the  multiple 
pass  behavior  somewhat. 

The  frequency  measurements  for  the  tapered  liner 
(Fig.  4.1)  show  little  variation  over  a  wide  voltage  range. 
Preference  for  a  particular  axial  mode  is  indicated  by  these 
measurements  and  by  the  output  pulse  locations,  with  respect 
to  the  peak  beam  voltage,  during  the  same  measurements. 
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50  kV/div,  20mV/div 

(a) 


50  kV/div,  lOOmV/div 

(b) 


Fig  4.4.  With  a  three  millimeter  thick  stycast  liner,  output 
occurs  when  the  beam  is  synchronous  with  one  of  the  axial 
modes  of  the  resonating  cavity,  not  necessarily  at  the  peak 
beam  voltage.  The  top  trace  shows  the  beam  voltage.  (One 
microsecond  per  division.)  (Taken  10-28-81.) 
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4.2  OUTPUT  POWER 
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Having  determined  that  the  feedback  and  growth  rate  are 

sufficient  for  oscillator  operation,  an  estimate  of  the  output 

power  is  necessary.  It  is  noted  that,  for  single  pass 

w*  »t 

amplification,  the  output  power  should  scale  as  e  .  Since 
an  input  signal  is  not  provided,  the  coefficient  of  this 
exponential  depends  on  the  degree  of  pre-modulation  undergone 
by  the  beam  during  acceleration  and  focusing,  and  is  difficult 
to  predict.  Multiple  passes  lead  to  much  larger  power.  The 
restrictions  on  exponential  growth  lead  to  an  estimate  of  the 
achievable  power. 

The  linearized  analysis  which  gives  the  collective 
dispersion  function,  Eqn .  3*12,  assumed  an  electron  density  of 
the  form  nQ+n ^  wt_kz\  n1<<nQ.12  As  the  radiation  field 
amplitude  grows,  so  does  the  bunched  portion  of  the  beam,  n1# 
Eventually,  ng=0,  contrary  to  the  assumption,  and  the 
dispersion  function  is  no  longer  correct. 

The  continuity  equation, 

3n/3t  +  v*(nv)  =  0,  (4.3) 


becomes 


(-iw  +  ikvQ  +  ikv1)n1  =  0. 


(4.4) 
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where  v= Vq+v 1 e i ^ wt_k .  Equation  4.4  is  satisfied  when  kv ^  is 
equal  to  the  original  detuning,  w-kv^. 

The  force  equation,  a(  rmv) /3t=-eEz ,  which  assumes  vQ  and 
v1  parallel  to  the  strong  guiding  magnetic  field,  now  gives 
the  electric  field  in  terms  of  known  quantities. 

^saturating  =  1Wt) t3  < w-kvQ ) v, t 
or 

lESatl  =  (m/e) ( w-kvQ )2/k  , 

where  Eqn .  3.28  is  used  to  evaluate  w-kv^.  When  the  longi¬ 
tudinal  electric  field  component  reaches  this  value,  the 
collective  gain  mechanism  begins  to  saturate. 

The  entire  beam  cross  section  will  see  a  field  at  least 
as  large  as  that  given  by  Eqn.  4.5b  when  the  field  at  r=0 
reaches  that  value.  With  this  convention,  the  saturating 
power  is  given  by  multiplying  Eqn.  3.8c  by  Eqn.  3-9e.  The 
coefficient,  B,  in  Eqn.  3.8c  is  set  equal  to  IE„  .|.  Fig.  4.6 
shows  the  behavior  of  the  saturating  power  along  a  dispersion 
curve.  Figures  4.7  and  4.8  show  the  voltage  and  current 
dependence.  If  only  the  beam  cross-sectional  r.m.s.  average 
Ez  is  required  to  reach  l^satl»  the  power  at  saturation  will 
be  less,  as  shown  in  Figures  4.9,  4.10,  and  4.11.  A  beam 
radius  of  six  millimeters  corresponds  to  a  one  millimeter  gap. 

At  a  fixed  voltage  the  growth  rate  is  proportional  to 

1  /3 

Ijj  »  regardless  of  the  choice  of  saturation  criterion. 


(4.5a) 

(4.5b) 
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Fig.  4.6.  The  saturating  power  and  efficiency,  defined  as  the 
saturating  power  divided  by  the  beam  power,  for  a  two 
millimeter  thick  boron  nitride  liner  in  12.5  millimeter  guide. 
The  beam  current  is  held  at  ten  amperes. 
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Fig.  4.7.  The  saturating  power  and  efficiency  from  th 
previous  figure  as  a  function  of  beam  voltage. 
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Fig.  4.8.  Saturating  power  and  efficiency  vs.  current  for  a 
two  millimeter  thick  boron  nitride  liner  in  12.5  millimeter 
guide  and  a  constant  125  keV  beam. 
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Fig.  4.9.  Saturating  power  and  efficiency  for  a  two 
millimeter  boron  nitride  liner.  An  effective  beam  density  is 
used  to  account  for  a  beam  diameter  of  six  millimeters.  The 
gap  between  the  beam  and  the  dielectric  is  about  one 
mill imeter . 
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Fig.  4.10.  The  saturating  power  and  efficiency  from  the 
previous  figure  as  a  function  of  beam  voltage. 
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Fig.  4.11.  Saturating  power  and  efficiency  vs.  current  for  a 
two  millimeter  boron  nitride  liner  in  12.5  millimeter  guide 
and  a  constant  125  keV,  six  millimeter  diameter  beam. 
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The  power  and  efficiency  are  proportional  to  and  Ife  3 , 

respectively . 

4.3  POWER  MEASUREMENTS 

Power  measurements  are  made  with  calibrated  lengths  of  a 
lossy  material  (Emerson  &  Cuming  eccosorb  AN72),  conventional 
Ka  band  (21-42  GHz)  variable  waveguide  attenuators,  and 
square-law  detector  diodes  in  Ka  band  mounts.  Radiation  is 
guided  from  the  circular  to  rectangular  transducer  through  the 
attenuation  to  the  diode.  As  the  power  increases,  attenuation 
is  added,  avoiding  nonlinear  diode  response. 

At  least  a  factor  of  two  uncertainty  is  present  in  these 
measurements.  The  waveguide  attenuators,  used  in  measurements 
with  quartz,  and  the  diode  mounts  are  designed  for  the  TE1Q 
mode  of  the  rectangular  guide.  The  actual  mode  structure  is 
probably  TM^.  A  fixed  amount  of  eccosorb  in  a  fixed  position 
was  used  in  the  quartz  measurements. 

In  measurements  with  boron  nitride  the  amount  of  eccosorb 
was  changed  between  readings;  the  variable  attenuators  were 
not  used.  This  technique  should  be  less  sensitive  to  mode 
structure.  However,  the  standing  wave  pattern  in  the 
detection  hardware  is  sensitive  to  the  position  of  the 
eccosorb.  Thus  measurements  are  uncertain  by  a  factor  of  the 
standing  wave  ratio. 


The  general  behavior  of  the  power  as  a  function  of  beam 


1 


K 

s 


current  for  the  two  millimeter  quartz  liner  is  shown  in 
Fig.  4.12.  A  calibration  was  obtained  with  the  same  liner 
geometry.  A  12  ampere,  115  kilovolt  beam  generated 
30  kilowatts  of  radiation  at  50  Ghz. 

Power  vs.  current  data  for  two  millimeter  boron  nitride 
are  given  in  Fig.  4.13.  The  frequency  is  21  GHz,  and  the 
maximum  power  is  36  kilowatts.  Powers  up  to  100  kilowatts 
(10%  efficiency)  were  obtained  for  three  millimeter  thick 
boron  nitride  in  the  12.5  millimeter  diameter  waveguide 
(Fig.  4.14).  Power  is  often  sufficient  to  cause  atmospheric 
breakdown  in  X-band  waveguide. 

The  kilowatt  level  power  measurements  correspond  well 
with  estimates  of  the  saturating  power.  Signal  levels  two  and 
three  orders  of  magnitude  lower  are  also  observed  and  may 
correspond  to  single  pass  operation.  Different  mechanisms  may 
dominate  at  different  times  during  the  output  pulse. 
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Fig.  4.12.  Saturating  power  and  measurements  (•)  for  two 
millimeter  quartz  with  a  150  keV  beam.  The  saturating  powe 
is  evaluated  assuming  a  six  millimeter  beam  diameter. 
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Fig.  4.13.  Saturating  power  and  measurements  (+)  for  two 
millimeter  boron  nitride  with  a  150  keV  beam.  The  saturating 
power  is  evaluated  assuming  a  six  millimeter  beam  diameter. 
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Fig.  4.14.  Saturating  power  and  measurements  (+)  for  three 
millimeter  boron  nitride  with  a  114  keV  beam.  The  saturating 
power  is  evaluated  assuming  a  six  millimeter  beam  diameter. 


CHAPTER  V 


5.1  CONCLUSION 

Knowledge  of  the  cold,  collective  growth  rate  for  the 
cylinder  is  useful  to  the  present  experimental  effort.  The 
onset  of  multiple  pass  operation  and  subsequent  saturating 
power  can  be  estimated. 

Planned  for  the  future  are  experiments  which  will  employ 
a  colder,  more  dense  (more  collective)  beam.  The  beam  cross- 
sectional  profile  will  be  more  uniform  in  the  2  direction. 
Modifications  to  the  pulse  forming  network  can  make  the  beam 
energy  more  uniform  in  time  and  provide  information  about  the 
possible  DC  operation  of  Cerenkov  masers.  With  these 
improvements  the  calculated  growth  rate  will  be  subject  to 
less  uncertainty. 

The  growth  rate  related  measurements  will  be  more 
accurate  as  well.  A  low  power,  swept  source  is  available  for 
passive  resonator  feedback  measurements,  and  methods  for 
controlling  the  feedback  externally  are  under  development. 
Better  power  diagnostics  and,  eventually,  beam  calorimetry 
will  be  implemented. 

All  of  this  information  can  be  used  to  optimize  the 
length  of  a  liner  for  a  given  cross-sectional,  beam-liner 


geometry  and,  with  a  little  more  effort,  to  design  a  beam- 
liner  geometry  of  variable  cross  section.  The  resulting  large 
energy  extraction  combined,  possibly,  with  a  beam  recovery 
system  should  lead  to  an  efficient,  high  power,  tunable, 
millimeter  source. 


APPENDIX  I 


A 1 . 1  RADIAL  TRIGONOMETRIC  FUNCTIONS 

Definitions : 


sn(u,v)  =  (JQ(v)Y0(u)  -  Y0(v)J0(u)  )/(2/irv)  (Al.la) 

cs(  u ,  v )  =  (YQ  ( v)  J  1  ( u)  -  JQ  (  v>  Y 1  (u)  )/(2/irv)  (Al.lb) 

Sn(u,v)  =  (J 1 ( v)Y1 (u)  -  Y1(v)J1(u))/(2/*v)  (Al.lc) 

Cs(  u  ,  v)  =  (J1(v)Yq(u)  -  Y^vJJqCu)  )/(2/irv)  (AT. Id) 


Sn(u,0)  =  lim(v-->0)  Sn(u,v) 

=  J^u)  lim(v-->0)(-Y1(v)/(2/  v)  )  =  J^u)  (A1.2a) 

Cs(u,0)  =  lim(v-->0)  Cs(u,v) 

=  JQ(u)  lim(v-->0)(-Y1(v)/(2/  v) )  =  JQ( u) . ( A1 . 2b) 

Identities: 

sn(u,u)  =  Sn(u,u)  =  0  (A1.3a) 

cs(u,u)  =  Cs(u,u)  =  /v/u  (  A 1 . 3b) 

sn( u ,v) Sn( u  ,v)  +  cs( u,v)Cs( u,v)  =  v/u.  (A1.H) 
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Derivatives : 


(d/du) sn( u,v)  r  cs(u,v)  (A1.5a) 

( d/du) cs( u ,v)  =  Sn(u,v)  (A1.5b) 

(d/du)Sn( u,v)  =  Cs( u ,v)-Sn( u ,v) /u  (A1.5c) 

( d/du) Cs( u  ,v)  =  -Sn(u,v)  (A1.5d) 

( 1 /(2/irqb)  )  (d/dq)  ( (2/irqb)  sn(  qa  ,qb) )  =  (a)cs-(b)Cs  (A1.6a) 


(1 /(2/wqb) )  (d/dq)  ( (2/nqb)  cs(  qa  ,qb)  )  =  -  (  a)  sn+ ( b)  Sn-cs/q .  (  A1 . 6b) 


APPENDIX  II 


A2.1  ANISOTROPIC  DIELECTRIC  LINERS 

Dielectric  tubes  and  sheets  are  manufactured  which  have  a 
lower  dielectric  constant  along  the  short  (50  axis  than  that 
along  the  other  two  axes  (Union  Carbide  "Boralloy”  boron 
nitride).  Such  a  material  suggests  the  possibility  of 
operating  at  a  low  phase  velocity  and  high  frequency  with 
milder  evanescence  of  the  longitudinal  electric  field. 
Preliminary  analysis  indicates  that  some  advantage  may  be 
gained  . 

The  wave  equation  for  Ez  becomes  more  complicated,  since 
the  substitution  v»£=Hirne  is  no  longer  possible.  The  fields 
are  now  given  as 


E 


z 


(  A)  sn(  q'  x  ,q  ’  b)  a  <  x  <  b 
( B)Cs( ipx ,0)  0  <  x  <  a 


E 


x 


(  ikq’ /qx^ ) ( A) cs( q* x ,q ' b)  a  <  x  <  b 
(-k/p)  (B)Sn(  ipx  ,0)  0  <  x  <  a, 


( A2. la) 


( A2. 1b) 


where  q'^  =  c  w^/c^-(e  /e  )k^,  q  ^  =  e  w^/c^-k^ ,  and  p^  =  k^-w^/c^. 

Z  2  a  XX 

The  dispersion  relation  is  now 
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( A2.2) 


(Qx2/Q' )tn(qa,qb)(Tn( ipa,0)/i)  =  exp. 

Equations  A2.1  and  A2.2  reduce  to  to  Equations  2.1  and  2.6 
when  ex=  e z . 

Fig.  A2.1  shows  dispersion  curves  for  isotropic  liners 
with  dielectric  constants  of  3.4  and  5.12  and  the  curve  for  an 
anisotropic  liner  with  ex=3.4  and  ez=5.12.  At  small 
wavenumbers  the  mode  of  the  anisotropic  liner  has  a  phase 
velocity  nearly  as  slow  as  the  phase  velocity  of  the  mode  of 
the  e=5.12  isotropic  liner.  At  large  wavenumbers  it 
approaches  the  e=3.4  phase  velocity.  The  group  velocity  of 
the  anisotropic  mode,  however,  is  quite  a  bit  larger  than  that 
of  either  of  the  isotropic  modes  in  the  range  of  wavenumbers 
shown.  The  experimental  wavenumbers  fall  in  this  range.  The 
growth  rate  of  the  collective  Cerenkov  instability  is 
connected  to  the  group  velocity  as  described  in  Section  3.1, 
and  may  also  be  larger  for  the  anisotropic  case. 


WAVENUMBER  (INVERSE  CM) 


Fig.  A2.1  Dispersion  curves  for  the  TMni  modes  of  two 
isotropic  liners  and  an  anisotropic  linei1  with  e  =  E  =  3.H  and 
e  =5.12.  The  liner  thickness  is  three  millimeters, rand  the 
guide  diameter  is  12.5  millimeters. 
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APPENDIX  III 


A3.1  SOFTWARE 

The  following  programs  and  catalogs  are  stored  in  the 

catalog  *D 1 1 1 40 : FELC AT  and  on  tape  Q1172. 

2SLABCAT:  Catalog  containing  double  slab  dispersion  program 
and  related  software. 

ANIS.DIS:  CYLDIS  modified  to  account  for  a  radial  dielectric 

constant  different  from  that  along  the  other  axes  of 
the  liner. 

COLGNEFF:  Program  to  calculate  the  cylinder  growth  rate  using 
an  effective  beam  density  (Eqn.  3-34)  to  account  for 
a  gap  between  the  liner  and  a  three  millimeter 
radius  beam.  The  guide  dispersion  curve  is  read 
from  files  CYLTMOn.  Liner  parameters  are  read  from 
CYLTMCO.  When  the  program  reaches  the  end  of  the 
dispersion  file,  it  attempts  to  continue  the 
calculation  by  using  an  approximation  (Eqn.  2.12). 
The  roots  of  ct(qa,qb),  necessary  for  the 


approximation,  are  read  from  CT. ZEROS.  Growth  rates 
are  output  to  COLGAINn. 

COLLGAIN:  Operates  exactly  as  COLGNEFF  assuming  no  gap. 

CTZEROER:  Calculates  roots  of  ct(qa,qb)  and  stores  them  in 
CT. ZEROS  (compiled). 

CUTOFFS:  Calculates  cutoff  frequencies  for  the  TEQn  and  TMQn 

modes  of  the  cylinder.  The  liner  parameters  and 
frequencies  are  stored  in  CYLTECO  or  CYLTMCO 
(compiled)  . 

CYLDIS:  Calculates  dispersion  curves  for  the  TEQn  and  TMQn 

modes  of  the  cylinder.  The  curves  are  stored  one 
point  (k  cm  ,  w  cm  )  to  a  line  in  terminal  format 
files  CYLTEOn  and  CYLTMOn,  n=1,...,9.  Liner 
parameters  and  cutoff  frequencies  are  read  from 
CYLTECO  or  CYLTMCO  (bound). 

IMSL8:  The  public  fortran  Bessel  routines  used  by  the 

various  cylinder  programs  were  recovered  from 
IMSLLIB**# : EDITI0N8 ,  appended  together,  compiled 
under  XF84,  and  saved  in  FELCAT  under  this  name. 

The  routines  are  accesse  by  basic7  programs  by  the 


statement  LIBRARY  "IMSL8" .  Fortran  programs  can  use 
the  statement  $LIBRARY  "IMSLLIB*## : IMSL8" . 

PLOT:  Calls  GRAPHLIB***  routines  to  perform  various 

plotting  tasks.  For  example,  empty  guide  dispersion 
curves,  beam  lines,  and  axes  appropriate  for 
dispersion  curves  can  be  drawn  with  simple  commands 
after  the  "dispersion"  or  "d"  command. 

POWER:  Calculates  growth  rate  using  effective  beam  density 

and  uses  it  to  calculate  the  saturating  power.  The 
saturating  electric  field  can  be  equated  to  the 
field  at  r=0  or  the  r.m.s.  average  field.  Either 
current  or  voltage  can  be  varied.  Input  files  are 
CYLTMCO  and  CYLTMOn.  POWERDAT  is  the  output  file. 

RAMANCAT:  Catalog  containing  programs  for  converting 
dispersion  curves  to  tuning  curves  when  the 
intersecting  beam  line  is  shifted  by  the  cyclotron 
frequency  or  the  pump  wavenumber. 

SCODELIB:  Catalog  containing  source  code  for  compiled 
programs . 

SPGAIN:  Calculates  single  electron  growth  rate.  Reads  liner 

parameters  from  CYLTMCO,  dispersion  curves  from 
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CYLTMOn,  and  attempts  to  approximate  using  CT. ZEROS. 
Growth  rates  are  output  to  SPGAINn. 

TUNECONV:  Converts  a  dispersion  curve  to  a  tuning  curve 
(  voltage , frequency)  . 

The  Herrmannsfeldt  electron  trajectory  program  and 
related  programs  are  stored  in  the  first  six  positions  on  tape 
Q-773  (format  Geos,  seven  track)  at  Kiewit  Computation  Center, 
Dartmouth  College. 
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Cerenkov  radiation  in  dielectric-lined  waveguides 
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Using  a  dielectric-lined  waveguide  in  combination  with  a  relativistic  electron  beam,  we  are  able  to 
produce  1  kW  of  Cerenkov  microwave  radiation.  We  have  observed  the  dependence  of  the  output 
microwave  frequency  on  the  beam  voltage.  In  addition,  we  have  detected  microwave  radiation  in 
excess  of  150  GHz,  corresponding  to  higher-order  mode  interaction. 

PACS  numbers:  32.30.Bv 

The  need  for  a  high-power  source  of  short-wavelength  of  using  modest  electron  beams  ( 1  SO  kV,  20  A)  in  studying 
microwave  radiation  for  use  in  communications,  plasma  the  Cerenkov  radiation  produced  in  dielectric-lined  wave¬ 
heating,  etc.,  has  motivated  much  experimental  work.  Re-  guides  and  in  scaling  such  a  device  up  to  higher  frequency 
cently,  work  has  been  directed  toward  making  the  cyclotron  operation. 

maser  (gyrotron)1  and  the  free-electron  laser2  into  viable  Cerenkov  radiation  is  produced  when  a  charged  parti- 

short-wavelength  microwave  generators.  The  Cerenkov  ra-  cle  is  directed  through  or  in  close  proximity  to  a  dielectric 
diation  experiments  presented  in  this  letter  represent  a  material,  at  a  speed  in  excess  of  the  speed  of  light  in  the 

slightly  different  approach  to  the  problem  of  producing  dielectric.  The  resulting  “shock”  radiation  is  generally  very 

short-wavelength  microwaves  at  significant  power  levels.  weak  and  is  emitted  in  a  continuous  spectrum.  However,  if  a 
Ginzburg  first  suggested  using  the  Cerenkov  effect  to  resonant  structure  such  as  a  dielectric-lined  waveguide  is 
produce  microwave  radiation  in  1947, 3  and  Abele4  pointed  used  in  combination  with  a  relativistic  electron  beam,  the 
out  that  Cerenkov  radiation  produced  in  a  waveguide  would  resulting  Cerenkov  instability  bunches  the  electron  beam  in 

lead  to  a  discrete  spectrum.  Since  then  several  experiments  the  waveguide  and  makes  it  possible  to  obtain  high-power 

have  substantiated  Ginzburg’s  predictions.3'7  These  early  microwave  radiation  at  discrete  frequencies, 
experiments  produced  only  small  amounts  of  power  and  The  frequency  of  the  Cerenkov  radiation  produced  in  a 

generally  used  prebunched  electron  beams.  In  1977,  Walsh,  dielectric-lined  waveguide  is  approximately  the  synchro- 
Marshall,  and  Schlesinger  used  an  intense  relativistic  elec-  nous  frequency  of  the  electron  beam  and  the  partially  filled 

tron  beam  (0.3  MeV,  10  kA)  in  combination  with  a  dielectric-  waveguide  modes,  as  shown  in  Fig.  1 .  T ransverse  magnetic 

lined  waveguide  and  produced  approximately  1  MW  of  mi-  waveguide  modes  with  aximuthal  symmetry,  TMoN  modes, 

crowave  radiation  at  60  GHz.8  The  present  experiments  will  be  the  modes  excited  by  the  axially  directed  electron 

have  been  performed  in  an  effort  to  investigate  the  viability  beam.  It  can  be  seen  that  by  varying  the  dimensions  of  the 
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FIG.  1.  Plot  of  the  dispersion  relation  of  a  dielectric-lined  waveguide  for 
TMon  waveguide  modes,  showing  the  synchronous  frequencies  of  the 
waveguide  modes  and  a  150-keV  electron  beam. 

dielectric  liner  or  by  changing  the  electron  velocity  the  syn¬ 
chronous  frequency  may  be  tuned  to  both  higher  and  lower 
values.  In  addition,  if  one  is  able  to  couple  to  higher-order 
waveguide  modes,  one  is  able  to  obtain  extremely  high  fre¬ 
quencies.  The  gain  obtainable  on  the  higher  branches  will, 
however,  usually  be  smaller  than  that  on  the  fundamental 
branch. 

A  schematic  diagram  of  the  electron  beam  generator 
and  resonator  structure  is  shown  in  Fig.  2.  The  cathode  is 
either  a  dispensor  or  a  barium-oxide-cathode,  which  pro¬ 
vides  a  solid,  cylindrical  electron  beam  of  up  to  30  A.  A  pulse 
from  a  pulse  forming  network  is  fed  into  a  high-voltage  pulse 
transformer  with  bifilar  secondary  windings.  The  output 
voltage  pulse  may  be  varied  between  0  and  -300  kV  and  is 
between  4  and  6  fi s  in  width.  The  machine  is  fired  1  to  7  times 
per  second  and  has  the  ability  to  operate  at  200  pps.  The 
electron  beam  is  accelerated  from  the  cathode  through  a  hole 
in  the  molybdenum  anode  and  into  the  quartz-lined 
(e  =  3.78)  cylindrical  waveguide.  The  beam  is  held  on  the 
axis  of  the  waveguide  by  an  axial  magnetic  field  of  between 
2.5  and  8  kG.  As  the  electron  beam  emerges  from  the  dielec¬ 
tric-lined  waveguide  it  follows  diverging  magnetic  field  lines 
until  it  is  deposited  on  the  waveguide  walls. 

Microwaves  produced  in  the  resonator  are  fed  out  the 
end  of  the  machine  where  power  and  frequency  measure¬ 
ments  are  carried  out.  Power  measurements  are  made  using 
crystal  detectors  and  attenuators  in  suitable  waveguide 


FIG.  2.  Schematic  diagram  of  the  experimental  apparatus  for  the  Cerenkov 
microwave  generation  experiments. 


FIG.  3.  Oscilloscope  traces  of  the  electron  beam  cunrent  and  voltage  and  the 
resultant  Cerenkov  microwave  pulse. 

mounts,  and  frequency  measurements  are  made  by  using  a 
series  of  high  pass  filters  or  by  using  a  microwave  Fabry- 
Perot  interferometer. 

Data  from  a  typical  shot  are  shown  in  Fig.  3.  The  micro- 
wave  pulses  usually  occur  on  or  near  the  peak  of  the  high- 
voltage  pulse  after  a  certain  voltage  threshold  has  been 
reached.  The  Cerenkov  threshold  for  quartz  (/?  2e  >  1, 
e  =  3.78)  is  about  85  kV,  but  we  do  not  observe  microwave 
output  until  voltages  of  about  120  kV  have  been  reached. 
The  value  of  the  guiding  magnetic  field  has  no  effect  on  the 
microwave  output  other  than  the  fact  that  it  must  be  present 
to  enable  the  electrons  to  pass  through  the  resonator.  Elec¬ 
tron  beam  currents  as  low  as  500  mA  are  sufficient  for  mi¬ 
crowave  generation,  while  higher  currents  serve  to  enhance 
the  power  and  pulse  width  of  the  microwave  signal.  For  a 
typical  150-kV,  10-A  electro”  beam,  we  detect  peak  micro- 
wave  powers  of  1  kW. 

The  inset  in  Fig.  4  shows  an  interferogram  obtained 
using  a  Fabry-Perot  interferometer  consisting  of  two 
screens,  one  fixed  and  one  movable.  The  interferogram  is  the 
result  of  hundreds  of  successive  shots  while  the  moving 
screen  is  slowly  swept  through  a  distance  of  3-5  cm.  The 
distance  between  successive  maxima  in  the  interferogram  is 
a  direct  measurement  of  one-half  the  free-space  wavelength 
of  the  radiation. 

A  series  of  such  interferograms  was  taken  for  numerous 
different  voltages.  A  plot  of  the  output  frequency  versus 
beam  voltage  is  shown  in  Fig.  4.  Included  in  Fig.  4  is  a  theo¬ 
retical  plot  of  the  synchronous  frequency  of  the  beam-wave¬ 
guide  system  as  a  function  of  beam  voltage.  We  note  that  the 


FIG.  4.  Plot  of  the  output  microwave  frequency,  obtained  using  the  Fabry- 
Perot  interferometer  [a  typical  interferogram  it  shown  in  the  inset),  as  a 
function  of  beam  voltage.  The  theoretical  plot  of  the  synchronous  frequen¬ 
cies  for  the  TMo,  mode  is  shown  for  comparison. 
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*  data  plot  has  the  same  general  shape  as  the  theoretical  one, 
but  that  it  is  displaced  to  higher  voltages.  We  find  that  by 
dividing  the  observed  voltage  by  a  factor  of  1 .3,  the  data  will 
fit  on  the  theoretical  curve.  As  of  yet,  we  have  found  no 
direct  reason  for  this  discrepancy. 

-  Using  high-frequency  filters,  we  have  been  able  to  de- 

IP  tect  radiation  in  excess  of  150  GHz  with  the  present  resona¬ 

tor  structure.  These  high-frequency  signals  are  probably  due 
to  the  interaction  of  the  electron  beam  with  higher-order 
waveguide  modes,  since  the  gain,  although  smaller,  is  not 
negligible.  The  power  of  the  150-GHz  signals  is  1  W  at  the 
diode  junction  and  may  be  much  higher  if  one  takes  into 
H  account  the  sensitivity  of  the  detector  diode  [  1 N53]  at  1 50 
GHz. 

In  conclusion,  the  above  experiments  have  demonstrat¬ 
ed  the  ability  to  use  the  Cerenkov  interaction  produced  when 
a  relativistic  electron  beam  is  accelerated  through  a  dielec¬ 
tric-lined  waveguide,  as  a  tunable  source  of  high-power, 
short-wavelength  radiation.  We  have  observed  the  basic  de¬ 


pendence  of  the  output  frequency  on  electron  beam  energy 
and  have  detected  microwave  signals  in  excess  of  1 50  GHz  at 
significant  power  levels. 

We  would  like  to  acknowledge  the  support  of  the  U.S. 
Army.  Grant  No.  DAAG-29-79-C0203. 


1 J.  L.  Hirshfield,  “Gyrotrons,”  in  Infrared  and  Millimeter  Wave  ,  edited  by 
Kenneth  J.  Button  (Academic,  New  York,  1979),  Vol.  1. 

2P.  Sprangle,  R.  A.  Smith,  and  V.  L.  Granatstein,  “Free  Electron  Lasers 
and  Stimulated  Scattering  from  Relativistic  Electron  Beams,”  in  Ref.  1. 
'V.  L.  Ginzburg,  Dok.  Akad.  Nauk  SSSR  56, 253  (1947). 

4M.  Abele,  Nouvo  Cimento  9,  207  (1952). 

SM.  Danos,  S.  Geschwind,  H.  Lashinsky,  and  A.  van  Trier,  Phys.  Rev.  92, 
828(1953). 

*P.  D.  Coleman  and  C.  Enderby,  J.  Appl.  Phys.  31, 1695  (1960). 

7R.  Ulrich,  Z,  Phys.  199,  171  (1967). 

* J .  E.  Walsh,  T.  C.  Marshall,  and  S.  P.  Schlesinger,  Phys.  Fluids  20, 709 
(1977). 


a 


009 


Appl.  Phys.  Lstt  M  (0),  15  April  19B1 


603 


FILMED 


